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The  monograph  discusses  the  mathematical 
principles  of  certain  allowed  variants  of  iner¬ 
tial  guidance  of  the  flight  of  ballistic  missiles t 
i.e.,  the  control  without  the  use  of  any  external 
information  (radio  signals t  radiation  of  stars 
and  others). 

It  is  assumed  that  the  providing  of  the 
assigned  flight  range  of  the  missile  is  produced 
as  a  result  of  the  well-timed  switching  off  of 
its  engine  by  a  signal  entering  from  the  computer. 
Fed  to  the  input  of  the  latter  are  current  read¬ 
ings  of  sensitive  elements  of  the  system  of  iner¬ 
tial  guidance  3  which  measure  the  apparent  accel¬ 
eration  of  the  missile  or  integrals  from  the 
apparent  time  acceleration.  The  control  should 
be  such  that  the  deviation  of  the  actual  motion 
of  the  missile  from  the  rated  does  not  have  an 
effect  on  the  range  of  its  flight.  For  this  mea¬ 
suring  instruments  of  the  system  of  inertial 
navigation  are  placed  aboard  the  rocket  in  a 
quite  definite  manner 3  and ,  specificallyt  the 
direction  of  their  axes  of  sensitivity  in  a  num¬ 
ber  of  cases  is  stabilised  by  means  of  gyroscopes . 
The  setting  of  the  necessary  directions  of  axes 
of  sensitivity  and  the  search  of  relations  which 
should  satisfy  current  readings  of  sensing  ele¬ 
ments  for  the  formation  in  the  computer  of  the 
command  for  the  switching  off  of  the  engine ,  and 
this  is  the  main  content  of  the  book.  Questions 
of  the  inertial  guidance  of  motion  of  the  missile 
in  q  lateral  direction  are  also  considered. 

The  monograph  is  intended  for  specialists  in 
the  field  of  the  theory  of  control  processes.  It 
can  prove  to  be  useful  in  the  investigation  of 
new  problems  of  this  discipline  and  also  as  a 
training  manual. 
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THE  AUTHOR'S  COMMENT 


This  book  came  about  as  a  result  of  thoughts  of  many  years  over 
the  general  problems  of  Inertial  guidance  systems  of  moving  objects. 
Come  of  them  are  presented  in  my  other  monograph,  published  In  19^3 
under  the  name  of  "Mechanics  of  Gyroscopic  Systems,"  In  which,  in 
accordance  with  Its  name,  considerable  space  has  been  given  also  to 
the  theory  of  sensing  elements  of  systems  of  inertial  guidance  —  a 
gyroscope  in  a  gimbal  suspension,  gyroscopic  stabilizer  e.nd  gyroscopic 
integrator  of  the  apparent  unit  accelerations.  In  this  book,  on  the 
contrary,  studied  basically  is  the  possibility  of  the  solution  of 
problems  of  inertial  guidance  of  ballistic  missiles  by  means  of 
different  combinations  of  gyroscopic  and  other  mechanized  sensing 
elements  and  computers.  Of  course,  questions  examined  here  exhaust 
by  no  means  the  whole  diversity  of  ideas  and  possibilities  of  auto¬ 
nomous  control  of  ballistic  missiles.  In  recent  years  in  the  theory 
of  specific  inertial  systems  great  successes  have  been  achieved. 

They  have  not  received  as  yet  the  proper  reflection  in  monographic 
literature.  It  can  be  hoped  that  for  the  assimilation  and  develop¬ 
ment  of  new  concepts  of  inertial  control  the  proposed  book  proves  to 
be  useful. 

The  author  thanks  A.  S.  Kachanov,  D.  M.  Klimov,  D.  F.  Klim, 

I.  S.  Kovner  and  M.  Ye.  Temchenko  for  their  valuable  advice  in  the 
editing  and  preparation  of  the  manuscript  for  press. 


INTRODUCTION 

Given  in  this  small  monograph  are  mathematical  bases  of  one  of 
the  possible  aspects  of  the  theory  of  flight  control  of  ballistic 
missiles  without  any  external  information  (radio  signals,  radiation 
of  the  stars  and  others).  The  initial  data  for  such  control  are 
readings  of  instruments  located  aboard  the  missile.  They  record 
both  the  orientation  of  the  missile  with  respect  to  fixed  directions 
(at.  fixed  stars)  and  the  difference  in  its  motion  from  free  flight 
without  resistance  of  the  medium.  The  effect  of  such  instruments 
is  based  on  the  phenomenon  of  inertia.  Because  of  this  the  mentione 
guidance  is  called  inertial.  Besides  Instruments  which  incorporate 
laws  of  mechanics,  in  systems  of  inertial  guidance  the  application 
of  instruments  founded  upon  other  physical  principles  (spins  of 
elementary  particles,  standing  waves  of  coherent  radiation  of  lasers 
and  others)  Is  possible. 

Creation  of  the  variant  examined  in  the  monograph  of  the  theory 
of  inertial  guidance  is  based  on  the  assumption  that  aboard  the 
missile  there  Is  stabilization  (specifically,  gyroscopic)  of  direc¬ 
tions  of  axes  of  sensitivity  of  meters  of  apparent  accelerations  — 
newt  unmeters  or  Instruments  directly  measuring  the  time  integral 
from  current  values  of  apparent  acceleration.  The  latter  are  called 
below  Integrators  of  accelerations;  they  can  also  be  named  impulso- 
meters  [Translator’s  note:  This  term  is  not  verified].  The  selec¬ 
tion  of  these  sensing  elements  was  determined  by  the  many  years  of 
interest  of  the  author  on  gyroscopic  instruments  and  other  mechanics 
measuring  devices. 
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The  instruments  have  a  different  kind  of  error.  The  effect  of 
the  latter  on  the  accuracy  of  flight  of  the  missile  requires  a 
special  Investigation,  which  goes  beyond  the  limits  of  this  monograph. 
Here  it  is  assumed  that  instruments  of  systems  of  inertial  guidance 
are  ideal,  i.e.,  they  function  as  though  without  errors.  The  same 
is  referred  to  servodrives  and  reading  devices'  and  also  any,  which 
accomplish,  specifically,  arithmetical  operations  and  the  integration 
of  current  values  of  different  parameters. 

Airborne  instruments  and  devices  should  not  be  excessively 
complex  and  bulky.  Therefore,  keeping  in  mind  the  use  of  simplified 
schemes,  with  the  solution  of  problems  of  inertial  guidance  certain 
fundamental  errors  can  be  allowed.  The  latter,  however,  should  not 
involve  great  deviations  in  the  missile  from  its  target.  The  selec¬ 
tion  of  an  appropriate  compromise  (at  which  deviations  are  small, 
and  the  system  of  inertial  control  is  not  too  complex)  constitutes 
one  of  the  basic  tasks  of  the  designer  of  a  specific  system  of 
inertial  guidance. 

Independently  of  the  form  of  the  guidance  system,  in  its  design 
one  should  proceed  from  the  fact  that  the  same  target  -  the  end  of 
the  free- flight  section  —  can  be  hit  by  the  ballistic  missile,  with¬ 
out  necessarily  moving  along  the  programmed  trajectory  but  according 
to  an  Infinite  set  of  other  adjacent  trajectories.  Therefore,  for 
an  accurate  hitting  of  the  target  it  is  not  necessary  that  at  the 
instant  of  termination  of  the  active  flight  section  of  the  missile 
its  coordinates  and  component  velocities  In  a  certain  system  of 
coordinates,  fixed  relative  to  the  earth,  would  be  equal  to  the 
calculated  provided  by  the  program  of  control.  The  last  observation 
is  very  important  since  the  providing  of  the  Just  mentioned  equalities 
places  before  the  system  of  control  of  flight  of  the  missile  very 
difficult,  not  always  feasible  tasks.  The  attaining,  a  rather 
accurate  determination  of  current  values  of  the  coordinates  them¬ 
selves  of  the  moving  missile  and  projections  of  its  velocity  is 
considerably  simpler.  It  is  possible,  therefore,  to  interrupt  the 
powered-flight  section  of  the  missile  at  exactly  that  instant  when 
the  totality  of  deviations  in  its  moving  coordinates  and  projections 
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of  velocity  from  the  appropriate  calculation  values  provides  sub¬ 
sequent  motion  along  one  of  the  trajectories  leading  to  the  target. 

The  determination  of  relationships  which  should  be  satisfied  by  the 
mentioned  deviations  at  instant  of  the  switching  off  of  the  engine, 
constitutes  the  main  task  of  the  stated,  theory  of  inertial  guidance 
of  flight  of  ballistic  missiles.  In  this  case  without  apparent 
deterioration  of  the  accuracy  of  hitting  of  the  missile  on  the  target, 
considerable  simplifications  of  guidance  systems  can  be  produced 
because  of  the  separate  control  of  the  range  and  lateral  motion  of 
the  missile  and  also  due  to  the  selection  of  specific  directions  of 
axes  of  sensitivity  of  inertial  instruments  and  the  forming  of 
signals  from  them. 

With  gyroscopic  stabilization  of  axes  of  sensitivity  of  newton- 
meters  and  integrators  of  accelerations,  it  is  natural  to  construct 
calculated  equations  and  equations  of  inertial  guidance  in  the 
system  of  coordinates  not  taking  part  In  the  rotation  of  the  earth. 

Let  us  r.ote,  however,  that  range  of  flight  of  the  missile  is  deter¬ 
mined  in  this  case  not  only  by  values  of  its  coordinates  and  com¬ 
ponents  of  velocities  relative  to  such  a  system  at  the  instant  of 
termination  of  the  power-flight  section,  but  also  of  duration  of  the 
latter. 

* 

In  the  computer  of  the  given  system  of  inertial  guidance  oil 
flight  range  according  to  current  readings  of  newtonmeters  and 
integrators  of  accelerations,  there  is  generated  a  certain  alternating 
magnitude,  uniquely  connected  at  each  instant  with  the  magnitude 
of  range  of  the  missile.  Understood  by  the  latter  is  that  magnitude 
of  range  which  is  obtained  if  the  thrust  of  the  engine  at  the  instant 
of  time  suddenly  becomes  zero.  To  provide  accurate  hitting  of  the 
missile  at  the  target,  the  switching  off  of  the  engine  should  be 
produced  at  the  instant  of  time  of  achievement  by  the  mentioned 
magnitude  of  assigned  value.  Since  this  magnitude  is  determined  by 
the  actual  course  of  the  change  in  coordinates  and  projections  of 
the  velocity  of  the  missile,  then  it  is  called  functional.  As  it 
will  be  shown  below,  it  is  possible  to  propose  different  forms  of 
functionals.  With  their  construction  in  the  computer,  subsequently. 


FTD-MT-24-291-70 


x 


there  will  not  be  taken  into  account  the  squares  and  products  of 
deviations  in  the  actual  coordinates  and  projections  of  the  velocity 
of  the  missile  at  the  instant  of  end  qf  operation  of  the  engine  from 
the  rated.  The  same  refers  to  deviations  of  moving  coordinates  and 
projections  of  the  velocity  of  the  rocket  and  also  to  the  time 
interval  between  the  necessary  and  calculated  instants  of  the  switch¬ 
ing  off  of  the  engine.  The  assignment  of  the  control  of  range  is 
thus  solved  in  a  linear  approximation.  This  should  not  lead  to 
great  errors  in  flight  at  the  assigned  range  of  the  missile  with 
well-controlled  thrust  of  the  engine  and  with  the  possibility  of 
the  switching  off  of  the  latter  after  the  feed  of  the  appropriate 
command  with  a  minimum  delayed  pulse.  In  the  case  when  control  of 
the  thrust  is  hampered,  subsequent  refinements  of  equations  deter¬ 
mining  the  termination  of  the  power- flight  section  are  necessary. 

These  problems  and  also  the  problem  of  inertial  control  without  the 
switching  off  of  the  engine  are  not  considered  in  this  monograph. 

Results  of  investigations,  given  in  this  monograph  indicate 
that  in  the  flight  of  a  missile  over  great  ranges  it  is  expedient 
to  have  on  the  gyroscopic  stabilizer  two  integrators  of  acceleration. 
Axes  of  sensitivity  of  these  integrators  shoula  :.e  definitely  oriented 
in  two  different  directions.  Fed  into  the  comparer,  which  determines 
the  instant  of  the  switching  off  of  the  engine,  are  the  current  value 
of  one  integrator  and  the  result  of  complementary  integration  of 
readings  of  the  other.  In  principle,  it  is  possible  to  bypass  one 
integrator  of  acceleration  if  in  the  process  of  flight  of  the  missile 
oroperly  both  parameters  of  the  Integrator  and  direction  of  Its  axis 
of  sensitivity  are  changed. 

The  use  of  two  stabilized  integrators  of  accelerations  or  one 
but  with  alternating  parameters  allows  constructing  the  so-called 
complete  functional  most  accurately  resolving  the  problem  of  the 
control  of  range  of  the  missile.  Subsequent  development  of  the 
theory  leads,  specifically,  to  the  basis  or  the  arrangement  of 
Integrators  of  accelerations  directly  aboard  the  missile  under  the 
condition  of  introduction  into  the  guidance  system  of  the  so-called 
standard  integrator  of  accelerations  with  accurate  execution  of  its 
commands . 

.  xi 
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The  acceleration  of  force  of  the  earth' 3  gravity  depends  on  the 
distance  between  the  missile  and  the  center  of  the  earth  and  also 
on  geographical  coordinates  of  the  missile.  Consequently,  forces  of 
gravity  acting  on  the  missile  in  its  actual  and  calculated  motions 
are  net  equal.  This  fact  can  condition  the  error  in  the  determining 
of  the  necessary  instant  of  switching  off  of  the  engine  even  during 
accurate  operation  of  the  integrators  of  accelerations  and  of  gyro¬ 
scopes  of  the  system  of  inertial  control.  For  the  complete  elimina¬ 
tion  of  this  error  the  airborne  computer  should  be  supplemented  by 
a  certain  comparatively  complex  integrating  equipment.  However,  a 
considerable  decrease  in  this  error  can  be  achieved  also  by  simpler 
means  -  because  of  a  small  change  in  the  orientation  of  the  integrator 
of  acceleration. 

Deviation  of  the  rocket  from  the  target  in  a  lateral  direction 
can  be  eliminated  with  the  help  of  introduction  into  the  system  of 
inertial  guidance  of  the  so-called  lateral  integrator,  which  controls 
the  motion  of  the  rocket  in  the  direction  of  the  normal  to  assigned 
programmed  plane  of  its  flight.  The  distinction  of  the  actual 
magnitude  of  duration  of  the  free-flight  section  from  Its  calculated 
value  leads  to  additional  lateral  deviations  of  the  missile  because 
of  the  rotation  of  the  earth.  However,  in  principle  it  is  possible 
to  avoid  errors  of  such  a  kind  if  we  use  current  values  of  integrators 
of  acceleration  of  system  of  range  control. 

Content  of  book  Is  the  following.  Chapter  I  gives  a  solution 
to  the  problem  of  inertial  guidance  of  ballistic  missiles  in  a 
simplified  formulation  for  the  purpose  of  explaining  basic  Ideas  and 
methods  examined  in  detail  in  subsequent  chapters.  Here  the  earth 
Is  taken  as  being  flat  and  not  rotating.  The  effect  of  the  atmosphere 
is  not  considered  and  the  force  of  gravity  is  considered  constant  in 
magnitude  and  direction.  It  is  natural  that  under  such  simplifying 
assumptions  the  motion  of  the  missile  in  the  free-flight  section  of 
its  flight  is  completely  described  by  known  equations  of  theoretical 
mechanics  about  the  motion  of  a  material  particle  in  a  vacuum  in  a 
uniform  field  of  the  force  of  gravity.  Further  it  is  explained 
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which  relationship  should  be  satisfied  by  small  changes  in  parameters 
of  the  end  of  the  power-flight  section,  i.e.,  coordinates  of  the 
missile  and  projections  of  its  velocity  at  the  instant  of  the  switch¬ 
ing  off  of  the  engine  in  order  that  it  hits  the  assigned  target  with 
an  error  of  not  more  than  the  second  order  of  smallness.  This 
relationship  serves  as  the  basic  point  for  the  formation  of  a  certain 
variable,  called  the  ballistic  function.  In  the  examined  case  the 
value  of  this  function  at  any  fixed  instant  of  time  is  the  approximate 
expression  for  the  error  in  the  range  of  the  rocket,  which  will  arise 
if  the  engine  of  the  missile  is  turned  off  at  precisely  this  instant. 
It  is  natural  that  for  the  hitting  of  the  missile  on  the  target  the 
engine  should  be  turned  off  at  the  instant  of  the  passage  of  current 
values  of  this  function  through  zero. 

To  construct  a  ballistic  function  aboard  the  rocket  in  the  form 
of  a  certain  electrical  or  mechanical  magnitude  with  the  variant  of 
the  inertial  system  selected  in  the  monograph  the  presence  of  special 
sensing  elements  is  proposed  -  newtonmeter  and  a  computer,  which 
contains  in  Its  composition  Integrating,  multiplying  and  summing 
elements.  In  the  example  of  the  simplest  meter  of  apparent  accelera¬ 
tion  Is  established  the  connection  between  the  magnitude  being  mea¬ 
sured,  the  real  acceleration  of  the  missile  and  acceleration  of  the 
force  of  gravity.  Concepts  of  apparent  velocity  and  apparent  .travel 
of  the  missile  are  introduced  and  it  is  indicated,  specifically,  that 
for  the  direct  measurement  of  the  apparent  velocity  it  is  necessary 
that  the  axis  of  sensitivity  of  the  integrator  of  accelerations  be 
stabilized.  It  turns  out  that  the  ballistic  function,  which  is 
referred  to  the  given  case  of  motion  of  the  missile  in  a  uniform 
field  of  gravity  without  allowing  for  resistance  on  the  side  of  the 
atmosphere,  can  be  constructed  aboard  the  missile,  using  only  the 
Integrator  of  acceleration  with  subsequent  integration  of  its  current 
readings  by  means  of  the  computer.  At  the  end  of  the  chapter  an 
analysis  of  this  question  is  given  both  from  an  analytical  and 
geometrical  point  of  view.  It  is  indicated  that  the  flight  range  of 
the  missile  will  be  changed  by  a  magnitude  of  the  second  order  of 
smallness  if  the  vector  of  its  actual  velocity  at  the  end  of  the 


power-flight  section  appears  the  same  as  that  in  the  calculated  case, 
and  the  position  of  the  missile  will  be  somewhat  displaced  in  a  cer¬ 
tain  definite  direction,  or,  on  the  contrary,  with  a  fixed  position 
of  the  missile  the  difference  between  the  real  and  calculated  vectors 
of  its  velocity  at  the  same  instant  of  time  proves  to  be  definitely 
directed.  In  the  examined  case  the  mentioned  directions  coincide 
with  each  other.  They,  specifically,  are  parallel  to  the  vector  of 
velocity  of  the  missile  at  instant  of  its  hitting  of  the  target  with 
calculated  motion.  In  general  when  the  curvilinearity  of  the  form 
of  the  earth,  the  effect  of  the  atmosphere  on  the  descending  branch 
of  the  free-flight  section  of  the  trajectory  and  heterogeneity  of 
the  field  of  the  earth's  gravity  are  taken  into  account,  such 
directions  also  exist,  but  they  are  not  parallel.  The  perpendiculars 
to  them,  the  so-called  A-  and  u-directions  can  serve,  as  it  is 
indicated  in  more  detail  in  Chapter  III,  for  orientation  of  axes  of 
sensitivity  of  integrators  of  accelerations  of  the  appropriate  system 
of  autonomous  control  of  the  flight  range  of  the  missile. 

Chapter  II  gives  ballistic  functions  of  various  kind  for  the 
general  case  of  flight  of  the  missile,  i.e.,  without  additional 
simplifying  assumptions,  which  took  place  in  Chapter  I.  Here  it  is 
impossible  to  express  the  magnitude  of  range  of  the  missile  in  the 
form  of  a  definite  equation,  which  would  contain  parameters  of  the 
end  of  the  power-flight  section,  for  example  in  the  starting  system 
of  coordinates.  Nevertheless,  correct  to  smallness  of  the  .second 
order  relative  to  differences  between  current  parameters  of  the  real 
motion  of  the  missile  and  parameters  of  the  end  of  the  power-flight 
section  of  its  calculated  motion,  it  is  possible  (with  the  help  of 
the  use  of  the  same  Instruments  as  in  Chapter  I)  to  determine  the 
error  in  the  range  which  would  take  place  with  the  switching  off  of 
the  engine  at  the  current  instant  of  time.  The  expression  for  such 
an  arbitrary  error  in  range  can  be  accepted  as  the  ballistic  func¬ 
tion.  Derivatives  of  the  magnitude  of  flight  range  of  the  missile 
entering  into  this ' function ,  according  to  parameters  of  the  end  of 
the  power- flight  section,  i.e.,  according  to  coordinates  and  pro¬ 
jections  of  velocity,  are  taken  at  calculated  values  of  the  latter. 
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Thus  they  are  constant  quantities,  which  can  be  determined  earlier 
for  each  assigned  case  of  flight  of  the  missile. 

With  the  inertial  control  of  range  the  current  value  of  the 
ballistic  function  should  be  determined  onboard  the  missile  itself. 
This  leads  to  the  necessity  of  continuous  determination  by  means  of 
the  computers  of  current  coordinates  of  the  missile  and  projections 
of  its  velocity  according  to  readings  of  newtonmeters  or,  in  other 
cases,  integrators  of  apparent  accelerations.  The  starting  system 
for  this  goal  proves  to  be  barely  adequate  because  of  the  necessity 
of  the  calculation  of  translational  and  coriolis  accelerations. 
Therefore,  it  is  more  expedient  to  pass  from  coordinates  and  pro¬ 
jections  of  velocity  of  the  missile  in  the  starting  system  of 
coordinates  to  appropriate  magnitudes  referring  to  the  nonrotational 
system  with  the  origin  at  the  center  of  the  earth.  As  a  result, 
after  the  rejection  of  terms  of  the  second  order  of  smallness  and 
terms  dependent  on  the  motion  of  the  missile  in  a  direction  per¬ 
pendicular  to  the  plane  of  flight,  the  so-called  initial  ballistic 
function  is  formed,  and  they  are  linear  with  respect  to  the  moving 
coordinates  and  projections  of  velocity  of  the  missile  and  time.  The 
expression  for  the  initial  ballistic  function  in  a  nonrotating  system 
of  coordinates  contains  a  derivative  of  the  magnitude  of  flight  range 
of  the  rocket  with  respect  to  duration  of  the  power-flight  section. 
However,  this  derivative  can  be  excluded  from  the  analysis  by  means 
of  the  use  of  some  auxiliary  relation.  The  latter  follows  from  the 
equality  of  the  actual  range  to  its  calculated  value,  if  coordinates 
and  projections  of  velocity  of  the  end  of  the  power-flight  section 
accurately  coincide  with  current  values  of  coordinates  and  projections 
of  velocity  of  any  particles  of  the  calculated  free-flight  section. 
Finally  the  basic  ballistic  equation  is  obtained  for  deterr'  on 
of  the  necessary  current  time  of  the  switching  off  of  the  engine  in 
flight  of  the  missile  at  the  assigned  range.  Determination  of 
current  coordinate's  and  projections  of  velocity  of  the  missiles, 
which  enter  the  left  part  of  this  equation,  requires  continuous 
integration  aboard  the  missile  of  appropriate  nonlinear  differential 
equations,  which  connect  the  second  derivatives  of  coordinates  with 


FTD-M7- 


yil-29i-70 


xv 


appropriate  projections  of  the  apparent  acceleration  and  accelera¬ 
tion  of  the  force  of  gravity.  However,  we  can  substantially  simplify 
such  a  problem,  if  into  the  basic  ballistic  equation  we  substitute 
expressions  of  coordinates  and  projection  velocities  obtained  as  a 
result  of  the  integration  of  approximate  linear  differential  equations 
in  the  so-called  isochronal  variations  of  coordinates  of  the  missile. 
The  latter  are  differences  between  the  acutal  and  calculated  values 
of  appropriate  coordinates  of  the  missile,  which  refer  to  the  same 
instant  of  time.  With  the  proper  arbitrary  extension  of  the  cal¬ 
culated  power-flight  section  of  flight  of  the  missile  during  the 
calculated  instant  of  the  switching  off  of  the  engine  isochronal 
variations  of  coordinates,  projections  of  velocity  and  projections 
of  apparent  acceleration  can  be  considered  small  magnitudes  during 
the  whole  interval  of  time  which  corresponds  to  the  real  power-flight 
section. 

Differential  equations  for  isochronal  variations  of  coordinates 
are  derived  In  Chapter  III,  and  their  approximate  integration,  allow¬ 
ing  r or  terms  appearing  as  a  result  of  the  heterogeneity  of  the  field 
of  the  earth's  gravity,  is  developed  in  Chapter  IV.  In  Chapter  III 
in  differential  equations  for  isochronal  variations  of  coordinates, 
the  mentioned  terms  are  completely  droped,  whereupon  variations  of 
projections  of  velocity  become  equal  to  variations  of  projections 
of  the  apparent  velocity  of  the  missile,  and  variations  of  the 
coordinates  themselves  -  to  variations  of  projections  of  the  apparent 
path.  This  allows  converting  the  basic  ballistic  equation  to  one  of 
the  forms  allowing  the  construction  of  current  values  of  its  left  side 
directly  onboard  the  missile  according  to  current  readings  of  integra¬ 
tors  of  accelerations. 

Different  orientation  of  axes  of  sensitivity  of  the  integrators 
of  accelerations  are  possible.  Depending  on  this,  the  computing 
part  of  the  system  of  inertial  control  of  the  range  and  the  number 
of  elements  entering  it  are  changed.  For  example,  with  orientation 
of  axes  of  sensitivity  according  to  invariable  X-  and  w-directions , 
which  was  already  mentioned  above,  it  is  possible  to  manage  with  only 


F?;>MT-2«-?91-70 


xvi 


one  element  of  repeated  integration  of  accelerations  unlike  the  case 
of  using  two  such  elements  necessary  in  the  arrangement  of  axes  of 
sensitivity  in  parallel  to  the  axes  of  the  nonrotating  system  of 
coordinates . 

It  proves  to  be  possible  to  construct  a  ballistic  equation  even 
with  the  help  of  a  single  meter  of  acceleration,  if  we  according  to 
the  earlier  assigned  law  change  the  orientation  of  its  axis.  In  the 
computer  there  should  be  provided  in  such  a  case  the  integration  of 
current  readings  of  this  meter,  preliminarily  multiplied  by  the 
coefficient,  whici  also  is  changed  with  the  course  of  time. 

The  problem  indicated  above  on  the  creation  of  a  ballistic 
equation  can  be  solved  approximately  also  with  the  help  of  one 
Integrator  of  acceleration,  as  this  is  indicated  at  the  end  of  the 
chapter  In  the  example  of  construction  of  the  system  of  inertial 
control  of  range  with  accurate  control  of  the  direction  of  thrust 
of  the  engine  according  to  indications  of  standard  Integrator  of 
accelerations.  In  this  case  the  axes  of  sensitivity  of  the  standard 
integrator  and  Integrator  of  the  system  of  control  of  the  range  are 
located  perpendicular  to  each  other.  Finally,  if  we  consider  the 
deviation  In  the  force  of  thrust  of  the  engine  from  the  longitudinal 
axis  of  the  missile  to  be  insignificant,  then  for  solution  of  the 
problem  of  control  of  range  one  can  use  current  readings  of  the 
integrator,  the  housing  of  which  is  directly  fastened  to  the  side  of 
the  missile,  and  the  axis  of  sensitivity  is  parallel  to  its  longi¬ 
tudinal  axis  (so-called  longitudinal  integrator). 

In  Chapter  IV  there  is  formed  the  integration  of  the  totality 
of  differential  equations  for  isochronal  variations  of  coordinates 
of  the  missile  with  the  approximate  calculation  in  them  of  terms 
which  characterize  changes  in  projections  of  the  force  of  gravity 
during  motion  of  the  missile  according  to  the  law  different  from 
the  calculated.  Corresponding  to  this,  the  basic  ballistic  equation 
derived  out  in  Chapter  II  is  converted  to  the  form  similar  to  that 
examined  in  Chapter  III.  Because  of  a  specific  method  of  integration 
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of  the  mentioned  totality  of  equations,  the  problem  on  inertial  con¬ 
trol  of  the  range  is  solved  by  the  same  means  as  it  was  earlier,  i.e., 
with  the  help  of  current  readings  of  integrators  of  accelerations  and 
computers,  which  contain  only  integrating,  multiplying  and  summing 
elements.  This  method  is  based  on  the  preliminary  replacement  of 
equations  of  the  mentioned  system  by  equivalent  integral  differential 
and  cumulative  relations  with  the  subsequent  introduction  into  the 
terms  being  integrated  instead  of  the  sought  functions,  which  repre¬ 
sent  current  values  of  coordinates  of  the  missiles  approximating 
their  polynomials  of  the  second  or  third  degree.  Coefficients  of 
approximating  polynomials  of  the  second  degree  are  selected  so  that 
the  polynomials  would  satisfy  the  same  initial  conditions  as  the 
sought  functions,  and,  furthermore,  coincide  with  the  latter  for  the 
Instant  of  termination  of  the  power-flight  section.  With  a  more 
accurate  solution  to  the  problem,  one  should  approximate  the  isochronal 
variations  of  coordinates  by  means  of  polynomials  of  the  third  degree. 
In  this  case  time  derivatives  are  equated  also  to  appropriate  pro¬ 
jections  of  velocity  of  the  missile  at  the  instant  of  termination  of 
the  power-flight  section.  As  a  result  for  the  determination  of 
variations  of  coordinates  and  projections  of  velocity  of  the  missile, 
algebraic  equations  with  coefficients  dependent  on  time  are  obtained. 
The  latter,  with  sufficient  accuracy,  can  be  replaced  by  constants 
equal  to  values  of  these  variations  at  the  Instant  of  termination  of 
the  power-flight  section  of  flight  of  the  missile.  Finally  the 
desired  variations  are  expressed  in  terms  of  readings  of  integrators 
of  accelerations  and  integrals  of  their  time  readings.  Thereby,  the 
subsequent  creation  of  different  forms  of  ballistic  functions  similar 
to  those  given  in  Chapter  III,  fundamental  has  no  basic  difficulties 
and  is  reduced  basically  to  certain  changes  in  coefficients  and 
directions  of  axes  of  sensitivity  of  Integra*. ors  of  accelerations. 

Finally,  Chapter  V  of  the  book  is  devoted  to  the  general 
problems  of  Inertial  control  of  the  lateral  motion  of  ballistic 
missiles.  The  removal  of  lateral  deviation  in  the  rocket  from  the 
target  requires  different  methods  than  the  providing  of  the  assigned 
range  of  its  flight,  where  the  required  accuracy  can  be  obtained 
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because  of  the  well-timed  switching  off  of  the  engine.  The  duration 
of  the  power-flight  section  with  flight  at  the  same  range  depends  for 
rockets  of  the  same  design  on  a  number  of  random  facts  and,  con¬ 
sequently,  is  itself  a  random  magnitude,  the  mean  value  of  which  is 
near  to  the  magnitude  of  duration  of  the  power-flight  section  of  the 
calculated  motion.  Therefore  it  is  sufficient  that  the  control  of 
lateral  motion  of  the  missile  would  be  especially  accurate  only  during 
the  interval  of  time  in  which  switching  off  of  the  engine  can  occur. 

The  current  removal  of  the  missile  from  the  programmed  plane  is 
connected  with  projections  on  the  normal  to  this  plane  of  its  apparent 
acceleration  and  acceleration  of  the  force  of  gravity  by  means  of  a 
differential  equation  of  the  second  order.  To  control  the  lateral 
motion  of  the  missile,  it  is  necessary  to  know  current  values  of 
this  removal  and  its  time  derivative.  They  can  be  obtained  by  means 
of  the  so-called  lateral  integrator  of  accelerations  and  simplest 
computer.  As  a  first  approximation  can  in  the  mentioned  differential 
equation  drop  terms  containing  the  projection  of  the  force  cf 
gravity,  whereupon  it  is  Immediately  integrated.  The  lateral  depar¬ 
ture  of  the  rocket  proves  to  be  equal  to  the  projection  of  its 
apparent  path  on  the  normal  to  the  programmed  plane  of  flight,  and 
the  value  of  velocity  of  this  departure  —  respectively,  to  the  pro¬ 
jection  of  the  apparent  velocity.  With  the  obtaining  of  more  accurate 
expressions  for  the  mentioned  sought  magnitudes,  the  approximation 
method  of  integration  of  differential  equations  given  in  Chapter  IV 
is  used. 

In  Chapter  V  it  is  also  indicated  how  one  can  use  readings  of 
the  instruments  of  inertial  control  of  the  range  in  order  to  remove 
the  additional  lateral  deviation  in  the  rocket  from  the  target  because 
of  noncoincidence  of  its  actual  and  calculated  motions  ir.  the  pro¬ 
jection  on  the  programmed  plane. 

In  the  appendix  of  the  book  an  analytical  derivation  of  cosines 
of  angles  between  axes  of  system  of  coordinates  with  a  finite  turn 
is  given.  The  table  of  cosines  is  used  in  Chapter  II. 
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C  H  A  P  i’  E  R  I 


THEORY  OF  INERTIAL  CONTROL  OF  FLIGHT  RANGE  OF  THE 
BALLISTIC  MISSILE  IN  SIMPLIFIED  FORMULATION 


S  1 .  Equation  Expressing  Flight  Range  of  the  Missile 
In  Terms  of  the  Cooralnate  ana  Projection  of 
Its  Velocity  at  the  End  of  the 
Power-Flight  Section 

The  general  formulation  of  the  problem  of  control  of  flight 
range  of  the  ballistic  missile  is  most  convenient  to  explain  i'  *he- 
following  simplest  example. 

Let  us  assume  that  the  earth  is  flat  and  does  not  rotate,  the 
atmosphere  surrounding  the  earth  is  absent, *and  the  acceleratior.  of 
the  force  of  the  earth1 3  gravity  /,  coinciding  in  this  case  with 
gravity  g,  is  constant  in  magnitude  and  direction.  Equations  of  the 
motion  of  the  missile1  on  the  free-flight  section  of  its  trajectory 
under  such  simplifying  assumptions  are  completely  integrated,  and 
thereby  the  whole  solution  to  the  problem  can  be  conducted  in  general 
form  up  to  the  end. 

Let  the  power-flight  section  of  flight  of  the  missile  proceed 
from  the  origin  0  of  the  fixed  system  of  coordinates  xy ,  the  x  axis 
of  which  is  horizontal,  and  the  y  axis  is  directed  along  the  vertical 
upward  (Fig.  1).  The  initial  velocity  of  the  missile  is  considered 
equal  to  zero. 

‘Here  and  further  the  term  "motion  of  the  missile"  is  understood 
as  motion  of  its  center  of  mass. 


i 
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Pig.  1. 


The  calculated,  or  programmed  motiqp  of  the  missiles  is  that  one 
at  which  its  coordinates  are  changed  in  accurate  conformity  with  their 
precalculated  values. 

Let  us  designate  coordinates  of  the  real  motion  of  the  missile 

at  instant  a,  termination  of  the  power- flight  section  of  its  flight, 

i.e.,  at  the  instant  of  complete  switching  off  of  the  engine,  by  x 

and  y,  and  the  same  coordinates  in  calculated  motion  —  by  x*  and  y* 

(Fig.  1).  Let  us  designate  by  u  and  »  projections  on  the  axis  x 

x  y 

and  y  of  the  velocity  of  the  rocket  at  the  instant  of  the  termination 

of  the  power-flight  section  of  its  real  motion  and  respectively  by 

v*  and  v*  —  their  calculated  values.  The  magnitude  of  duration  of 
x  y 

the  power-flight  section  of  the  actual  flight  of  the  missile  o  in 
genera i  is  distinguished  from  its  calculated  value  a*. 

The  integration  of  equations  of  motion  of  the  missile  in  the 
pi-'v..,.  Might  section  of  motion 

^-o,  (l.i.D 

immediately  leads  to  the  following  expressions  for  current  projec¬ 
tions  of  velocity  v  ft)  and  v  ft): 
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MO 


rfy(l) 

~3T~ 


*»  — fO 


and  its  moving  coordinates  x(t)  and  y(t): 

y(0  — y  +  V— 


(1.1.2) 


(1.1.3) 


Time  t  is  read  off  here  anew,  i.e.,  from  the  instant  of  termination 
of  the  power-flight  section  and,  consequently,  from  the  beginning  of 
the  power-flight  section  of  motion  of  the  missile.  Arbitrary  con¬ 
stants  with  integration  are  selected  so  that  there  would  take  place 
the  apparent  equalities 

MO)—;**,  ,M°)  ->».  .*(0)-*,  y(0 )-y.  (1.1.0) 

Considering  that  the  free-flight  section  is  finished  at  instant 
t  =  t  at  point  C  —  intersection  of  the  trajectory  of  the  missile 
with  the  x  axis,  it  is  easy  to  determine  the  connection  between  the 
range  of  its  flight  l  and  appropriate  parameters  x,  y,  v  and  vt  - 

x  y 

end  of  the  power-flight  section.  Setting,  first  of  all,  in  the 
second  equation  (1.1.3)  t  =  r  and  y(x)  =  0,  we  obtain  the  following 
quadratic  equation  for  determination  of  the  duration  of  the  free- 
flight  section  t: 

ft*— 2»ft  ~2y  —  0.  (1.1.5) 

Hence 

T-7-(«v  +  K*T+2fy).  (1.1.6) 

since  one  takes  the  positive  root  of  equation  (1.1.5).  Having 
substituted  now  the  found  value  t  into  the  first  equation  (1.1.3) > 
we  obtain  the  desired  expression  for  the  range  of  the  rocket  as 
functions  of  parameters  x,  y ,  v£  and  v —  end  of  the  power-flight 
section,  namely: 

/ -«(t) -*  +  -!£.(*,  + J/p»+2fy).  (1.1.7) 
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§  2.  Condition  of  Invariability  of  Flight  Range  with 
Little  Distlnction~ln  Coordinates  and  Projections 
of  Velocity  of  the  Missile  from  Their 
Calculated  Values  at  the  Instant 
of  Termination  of  the 
Power-Plight  ' 

Section 


Let  us  assume  that  the  actual  values  of  parameters  of  the  power- 

flight  section  of  flight  of  the  missile  i,  y,  v  and  y  are  distin- 

x  y 

guished  from  their  calculated  values  by  small  magnitude 


r— x\  Ay  —  y  ~  y*.  Ac,  «»  n*,Ap,  =  wv—v\. 


(1.2.1) 


The  flight  range  of  the  rocket  Z,  in  accordance  with  that  given  in 

§  1,  and,  specifically,  according  to  equation  (1.1.7),  is  a  function 

of  magnitudes  x,  y,  y  and  y  ,  i.e., 

x  y 

/  —  /  (x,  y,  P,).  (1.2.2) 

In  accordance  with  equality  (1.2.1),  we  have 

!-/(**+ Ax,  jr*+Ajr,  v\  +  A*„  v\  +  A»f).  (1.2.3) 


Expanding  now  the  right  side  of  equation  (1.2.3)  into  Taylor  series 
for  the  function  of  many  variables,  we  obtain  correct  to  smallness 
of  the  first  order  relative  to  magnitude  Ax,  Ay,  Ay^  and  Ay^ 

f  - 1  <«",  F ,  £  * #r)  +  ^  Ax  +  Ay  +  Ar,  +  *L  Arr  ( 1 . 2 .  H ) 


Partial  derivatives  of  function  l(x,  y,  v  ,  ,  )  with  respect  to 

x  y 

variables  x,  y,  y  and  y  are  themselves  functions  of  the  same 

x  y 

variables.  In  this  case  they  can  be  presented  in  evident  form  as 
a  result  of  the  differentiation  of  function  (1.1.7)  according  to 
the  appropriate  variable,  namely: 


(1.2.5) 


According  to  rules  of  the  composition  of  Taylor  series,  in 

equation  (1.2.4)  arguments  of  aforementioned  partial  derivatives 

should  be  assumed  correspondingly  equal  to  magnitudes  x*t  y*t  v* 

and  v *.  Thus  they  are  for  the  selected  calculated  motion  of  the 
d 

missile  certain  constants.  Taking  into  account,  furthermore,  that 
expression 

(1.2.6) 

is  the  calculation  flight  range  of  the  missile,  let  us  transform 
equation  (1.2.4)  co  the  form 


**  -  **  +t«+  VW+W)  *»■+ 

+^(,+F5feK 


where  the  difference 


a/-/— r 


(1.2.7) 


a. 2. 8) 


Is  the  change  in  flight  range  of  the  missile  because  of  noncoincidence 

of  magnitudes  of  parameters  of  the  end  of  the  power- flight  section 

x,  y*  and  y  with  their  calculated  values  x*.  y*.  v*  and  v*. 
a  x  y  *  a  *  x  y 

The  expression  for  AZ,  allowing  for  equation  (1.1.6),  can  be 
represented  now  in  the  form 


where 


hi  -  A*  +  *Ay  +  +  ***,). 


a. 2. 9) 


Let  us  note,  that  coefficient  k  has  a  simple  geometric  meaning.  It 
Is  easy  to  show,  using  equations  (1.1.2)  and  (1.1.6),  that  it  is 
equal  to  the  magnitude  of  the  angular  coefficient  of  the  normal  v  to 
the  calculated  trajectory  of  motion  of  the  missile  at  point  C, 
its  Intersection  with  the  x  axis  (Fig.  1).  Thus, 


(1.2.11) 


where  a  -  angle  between  the  mentioned  direction  v  and  x  axis  (and 
also  between  the  vertical  line  and  vector  of  velocity  of  the  missile 
at  the  Instant  of  Its  incidence  on  the  earth). 

Prom  equation  (1.2.9)  it  follows  that  correct  to  smallness  of 
the  first  order  Inclusively  the  range  of  the  missile  remains  the 
very  same  if  changes  in  parameters  Ax,  Ay,  Lvg  and  -  end  of  Its 
power-flight  section  in  comparison  with  their  calculated  values,  will 
be  subordinate  to  the  condition 

+  +  t'(A».  +  *A»,)-0.  (1.2.12; 

It  is  important  that  this  condition  is  the  linear  relation  with  con¬ 
stant  coefficients  known  earlier  for  the  selected  flight  of  the 
missile . 

Tht  duration  o'"  the  power-flight  section  of  flight  o,  in  con¬ 
trast  to  parameters  cf  its  end  x,  y ,  vg  and  ,  does  not  play  any 
role  in  the  determination  of  the  magnitude  of  range  l.  Therefore, 
there  is  no  meaning  to  strive  with  execution  of  the  specific  program 
of  motion  of  the  missile  for  the  rated  value  of  duration  o*  of  the 
power-flight  section.  On  the  contrary,  by  slightly  lengthening  or 
shortening  the  power-flight  section.  It  is  possible  to  obtain  a 
fulfilling  of  Just  the  given  condition  (1.2.12)  and,  consequently, 
provide  the  calculated  range  of  the  missile  at  not  exactly  an  accurate 
sequence  of  the  actual  coordinates  to  their  programmed  values. 

5  3 •  Error  in  Flight  Range  of  the  Missile  Presented 
in  the  Form  of  a  Function  of  the  Duration 
of  the  Power-Flight  Section 

Let  us  assume  that,  as  also  in  5  1,  x(t)  and  y(t)  are  functions 

which  are  the  actual  (not  calculated)  change  in  time  of  coordinates 

of  the  missile  on  the  power-flight  section,  a  v  ft)  and  v  (tj,  re- 

x  y 

spectively,  are  projections  of  its  velocity,  and  the  Instant  t  »  0 
at  this  time  corresponds  tc  the  beginning  of  the  power-fiight  section. 


fc 


Let  us  substitute  into  equation  (1.2.9)  for  AZ  quantities  Ax, 

Ay,  Au  and  Ay  respectively  by  differences 
x  y 

*(0— *\  y(0— y\  *,(')—«£.  *»(0— «v  (1.3.1) 

As  a  result  let  us  obtain  the  time  function 

* 

.(*)-*(<)-* I  + 

*Jl.  (1.3.2) 

which  we  call  ballistic  function. 

At ' the  instant  of  the  switching  off  of  the  engine  t  »  a*  this 
function,  according  to  equation  (1.2.9),  turns  into  A l  and,  con¬ 
sequently,  determines  the  error  in  flight  range  of  the  rocket  because 
of  an  accurate  execution  of  the  program  of  the  power-flight  section 
but  accurate  observance  of  its  duration.  However,  if  we  interrupt 
the  power-flight  section  at  that  instant  when  the  ballistic  function 
e (t)  turns  into  zero,  then  the  error  into  the  range  proves  to  be  a 
magnitude  of  the  second  order  of  smallness,  and  the  desired  accuracy 
of  flight  of  the  missile  will  be  achieved.  Hence  it  follows  that 
the  problem  of  inertial  control  of  the  flight  range  of  the  missile 
can  be  reduced  to  the  construction  of  current  values  of  the  ballistic 
function  directly  aboard  the  missile  and  to  the  switching  off  of  its 
engine  with  the  passage  of  function  z(t)  through  zero.  For  this 
purpose,  besides  the  summing  and  multiplying  devices,  the  presence 
aboard  the  missile  of  special  instruments  which  measure  the  apparent 
acceleration  with  subsequent  integration  of  their  current  readings, 
or  instruments  directly  recording  the  time  integral  from  the  apparent 
accelerations  -  integrators  of  accelerations  is  necessary.  Meters 
of  apparent  accelerations  will  subsequently  be  called  newtonmetere . 

§  *4.  Connection  Between  the  Acceleration  of  the  Missile 
in  Its  Motion  Relative  to  the  Fixed  System  of 
Coordinates  and  Reading  of 
Newtonmeter  5et  on  ft 


In  the  system  of  control  of  the  flight  range  of  the  ballistic 


missile  when  using  newtonmeter  one  should  keep  in  mind  that  the  latter 
can  measure  only  the  apparent  and  not  the  real  acceleration  of  that 
place  of  the  missile  where  they  are  located.  The  apparent  accelera¬ 
tion  is  usually  called  the  difference  Detween  the  acceleration  of 
any  points  in  the  fixed  system  of  coordinates  and  the  acceleration 
of  the  force  of  gravity.  Specifically,  standard  uniaxial  newton- 
meters  should  measure  the  projection  of  the  apparent  acceleration 
on  its  axis  of  sensitivity  v,  i.e.,  magnitude 

(1.4.1) 

where  —  projection  on  axis  v  of  the  real  acceleration  of  the 
newtonmeter  relative  to  the  fixed  system  of  coordinates  and  — 
projection  on  the  same  axis  of  acceleration  of  the  force  of  the 
earth's  gravity. 

Let  us  explain  equation  (1.4.1)  in  the  example  of  the  newton¬ 
meter,  the  sensing  element  of  which  is  a  small  weight  of  mass  m 
attached  to  the  end  of  the  spring  of  rigidity  c  (Fig.  2).  The  other 
end  of  the  spring  is  sealed  in  the  housing  of  the  newtonmeter.  The 
small  weight  can  be  moved  within  the  housing  on  straight  line  v, 
which  is  the  axis  of  sensitivity  of  such  a  meter  of  apparent  accelera¬ 
tion. 


Fig.  2. 


The  equation  of  motion  of  the  Small  weight  relative  to  the 
housing  of  the  newtonmeter,  if  we  disregard  the  mass  of  the  spring 
and  friction  of  small  weight  about  the  internal  cylindrical  surface 
of  the  housing,  can  be  presented  in  the  form 


—  m-gjr  +  iM*:  ""  mfi  +  c6. 


(1.4.2) 


Here  6  —  movement  of  the  small  weight  from  that  position  E*at  which 
the  spring  has  not  been  stretched;  w*  -  projection  on  direction  v  of 
translational  acceleration,  i.e.,  acceleration  with  respect  to  the 
fixed  system  of  coordinates  of  that  place  cf  the  housing  where  at  the 
giver,  instant  the  small  weight  is  located  (projection  of  the  coriolis 
acceleration  to  direction  v  is  equal  to  zero). 

With  translational  movement  of  the  housing  of  the  newtonmeter 


(1.4.3) 


where  uv  -  acceleration  of  the  place  of  fastening  of  the  spring  to 
the  housing. 

If,  however,  the  housing  of  the  newtonmeter  has,  furthermore, 
angular  motions,  then  equality  (1.4.3)  should  be  considered  as 
approximate.  However,  because  of  the  comparatively  small  dimensions 
of  newtonmeters  the  difference  between  and  is  important.1 


Let  us  assume  that  the  frequency  p  of  natural  oscillations  of 
the  small  weight,  determined  cy  equation 


V±- 


*It  is  possible  to  show  that  this  difference  is  equal  to  the 

product  u>q c  where  -  projection  cf  angular  velocity  of  the  nousing 

cf  the  newtonmeter  on  the  plane  perpendicular  to  the  axis  ■>  and  p  — 
distance  between  the  small  weight  and  point  of  attachment  of  the 
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is  sufficiently  great.  Then  the  amplitude  of  these  oscillations, 
conditioned  basically  by  the  change  in  acceleration  of  the  housing 
of,  the  newtonmeter  will  be  small.  In  this  case  the  term  -  m 
in  equation  (1.4.2)  can  be  dropped,  and  allowing  for  equality  (1.4.3) 
we  can  obtain  relation 

+  (1.4.5) 

The  member  c6  of  equation  (1.4.5)  is  the  elastic  force  of  the 
spring.  Deformation  6  (see  Pig.  2)  will  be  considered  positive  if 
the  small  weight  is  displaced  in  a  negative  direction  of  axis  v. 

When  6*0  the  spring  is  not  stretched,  and  the  small  weight  is 
located  at  position  E  (see  Pig.  2).  The  observable  shift  in  the 
small  weight  from  position  E  can  be  graduated  so  that  it  would 
directly  measure  a  certain  magnitude  connected  to  the  deformation 
of  the  spring  by  relation 


(1.4.6) 

Replacing  here  the  magnitude  of  deformation  of  the  spring  6  by  its 
expression,  by  the  following  from  equation  (1.4.5),  we  will  arrive 
at  equation  (1.4.1).  Thus,  the  shift  in  the  small  weight  is  pro¬ 
portional  to  the  projection  of  its  apparent  acceleration  on  the 
direction  v,  i.e.,  on  the  axis  of  sensitivity  of  the  newtonmeter. 

The  current  reading  of  the  integrator  of  the  apparent  accelera¬ 
tions  or  simply  integrator  of  accelerations  (it  can  also  be  called 
impuleome ter)  is  a  time  integral  from  the  projection  of  the  apparent 
acceleration  on  its  axis  of  sensitivity,  i.e.,  magnitude 

i 

(1.4.7) 


If  the  axis  of  sensitivity  of  the  integrator  of  acceleration 
retains  a  fixed  direction  with  respect  to  the  fixed  system  of 
coordinates,  then 
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(1.H.8) 


where  v^(t)  -  projection  on  the  directi-on  v  of  velocity  of  the  hous¬ 
ing  of  the  integrator  in  the  same  fixed  system.  Therefore,  sub¬ 
stituting  into  equation  (1.4.7)  the  expression  for  av  from  equation 
(1.4.1)  and  using  equality  (1.4.8),  we  obtain  the  relation 


V.(0  -  B,<0-  0.(0) -  j  /.(/)<*/,  (1.4.9) 

which  connects  the  projection  of  velocity  v^ft)  with  the  reading  of 
the  integrator  of  acceleration  V^(t) .  The  latter  is  called  in  this 
case  the  projection  of  the  apparent  velocity  on  direction  v. 

Function  (tj  which  stands  under  the  sign  of  the  integral  in 
the  right  side  of  relation  (1.4.9)  when  the  earth  is  not  proposed  to 
be  flat,  is  changed  with  the  course  of  time  because  of  a  change  in 
the  position  of  the  missile  relative  to  the  earth. 

With  alternating  orientation  of  the  axis  of  sensitivity  v,  for 
example,  when  the  integrator  of  acceleration  is  located  directly 
aboard  the  missile,  equations  (1.4.8)  and  (1.4.9),  of  course,  are 
not  applicable,  and  the  reading  of  the  integrator  of  acceleration  is 
no  longer  equivalent  to  the  corresponding  projection  of  the  apparent 
velocity . 


5  5 .  Construction  of  a  Ballistic  Function  by  Means  of 
.  Two  Integrators,  of  Accelerations  and  a  Computer 


In  this  section  let  us  examine  the  use  in  the  system  of  inertial 
control  of  the  flight  range  of  two  Integrators  of  accelerations. 

Let  us  assume  that  the  integrators  are  stabilized  so  that  the  axis 
of  sensitivity  of  one  of  them  during  the  whole  power-flight  section 
of  the  missile  would  remain  parallel  to  the  horizontal  x  axis,  and 
the  axis  of  sensitivity  of  the  other  would  be  directed  parallel  to 
the  vertical  y  axis. 
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In  the  examined  case,  i.e.,  under  the  assumption  of  a  nonrotat¬ 
ing  flat  earth,  projections  of  accelerations  of  the  force  of  gravity 
on  the  x  and  y  axes  are  expressed  by  equations 


/.-0.  (1.5.1) 

Consequently,  in  accordance  with  formula  (1.4.1),  projections  of  the 
apparent  acceleration  on  the  same  axis  are  the  quantities 


where 


«,-«M0.  «*-«MO+r.  (1.5.2) 

(1.5.3) 


are  corresponding  projections  of  the  real  acceleration  of  the  missile 
in  the  fixed  system  of  coordinates  xy .  Here  x(t)  and  y(t),  Just  as 
earlier,  are  current  coordinates  of  the  missile  on  the  power-flight 
section  of  its  flight. 

The  current  readings  of  the  integrators  with  axes  of  sensitivity 
parallel  to  axes  of  coordinates  x  and  y,  according  to  equation 
(1.4.7),  are  quantities 


V. 


(1.5.4) 


which  are  projections  of  the  apparent  velocity  of  the  missile, 

respectively,  on  the  axis  of  the  fixed  orientation  *  and  y.  In 

virtue  of  equations  (1.5.2)  and  (1.5.3)  they  are  connected  with 

projections  v  (t)  and  v  (t)  of  the  real  velocity  of  the  missile 
*  x  y 

relative  to  the  system  of  coordinates  xy  by  relation 

F.(0  -  MO.  F,(0  -  MO +(*■  (1.5.5) 

Into  the  latter  constants  of  integration  are  absent,  since  projec¬ 
tions  of  the  real  velocity  of  the  missile  at  the  initial  instant 


of  the  power-flight  section,  i.e.,  when  t  •  0,  are  equal  to  zero. 

Coordinates  of  the  missile  x(t )  and  y(t)  at  this  Instant  are  also 

equal  to  zero.  Therefore,  as  a  result  of  the  integration  with 

respect  to  time  of  current  indications  V  ft)  anc  V  (t)  of  corre- 

x  y 

sponding  integrators  of  accelerations  will  be  reduced  on  the  ba3is 
of  equations  (1.5*3)  and  (1.5*5)  to  the  following  relation: 

J.(0-*(0+4f.  (1*5*6 

In  them  time  functions 

<  t 

*.<*) -  j  5,(1) -  j  V9{t)it  (1.5*7; 

can  be  called  projections  of  the  apparent  path  of  the  missile. 

By  means  of  relations  (1.5*5)  and  (1.5*6)  one  can  express! 

coordinates  x(t)  and  y(t )  in  terms  of  functions  5  (t)  and  S 

x  y 

and  projections  vx(t)  and  v^(t)  of  its  real  velocity  -  in  terms  of 
current  readings  of  integrators  of  accelerations  V(t)  and  V,  (tv. 

If  further  we  substitute  these  results  into  equation  (1.3.?)  for 
the  ballistic  function  c(t) ,  then  after  the  simplest  conversions  we 
will  obtain  for  it  the  following  expression: 

•(I)  -  5.{I) +*5,(0 + f  IM0+*V,(0I  - 

-«*-V-t>;+^-*J kgt-k%,  U.5.8) 

where  the  coefficient  k,  as  before,  is  determined  by  equation 
(1.2.10). 

For  the  construction  aboard  the  missile  of  a  current  value  of 
the  ballistic  function  c(t)  in  the  form  of  an  alternating  mechanical 
or  electrical  magnitude  tn  the  examined  case  the  presence  of  two 
integrators  of  accelerations,  clocks  and  computer  is  necessary.  The 
composition  of  the  latter,  specifically,  should  include  two  elements 
of  additional  integration  of  readings  of  integrators  of  accelerations 
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§  Control  of  Flight  Range  of  the  Missile  by  Means 
of  a  Single  Integrator  of  Accelerations  with 
Inclined  Orientation  of  Its  Axis 

An  attentive  examination  of  expression  (1.5.8)  for  the  ballisti 
function  c(t)  leads  to  the  conclusion  that  sums 

+  0.  J.(i)  +  W,<  Or  (1.6.1) 

entering  into  its  composition  can  be  formed  aboard  the  missile  by 
means  of  only  the  integrator  of  accelerations  with  subsequent  addi¬ 
tional  time  integration  of  its  instant  readings.  The  axis  of 
sensitivity  v  cf  such  an  integrator  should  be  inclined  during  the 
whole  time  of  the  power-flight  section  of  flight  of  the  missile  to 
the  horizontal  axis  x  at  a  constant  angle  a,  the  magnitude  of  which 
is  connected  witn  coefficient  k  by  relation  (1.2.11).  Thus  (see  §  2 
of  this  chapter),  the  direction  of  the  axis  of  sensitivity  of  the 
mentioned  integrator  is  perpendicular  to  the  vector  of  velocity  of 
•no  missile  at  the  instant  of  reaching  by  it  of  the  target  C  with 
calculated  motion  (Fig.  1). 

Let  us  note,  first  of  all  (Fig.  3),  that  for  projections  on 
the  direction  v  of  the  acceleration  of  the  force  of  gravity  and  real 
acceleration  of  the  missile  we  have  equations 

U  —  —  w,(/)eo»a  +  w,(*)*in<*.  (1.6.2) 

Consequently,  according  to  equation  (l.ft.l),  we  have 

«,  -  »,(!)— /.<*)  -  »•(!) «»«  +  to (0  +  g]  (1.6.3) 

or,  taking  into  account  equality  (1.5.2), 

•,(0-*,(0«m«+*»(0»1“«'  Cl. 6. ft) 

The  last  relation  indicates  that  the  apparent  acceleration  is  a 
vector  wh'  ?e  components  are  parallel  to  axes  of  coordinates  x  and  y , 
respectively  equal  to  quantities  a  (t-1  and  a  ( • 

*  v 


Integrating  expression  (1.6.4)  for  the  projection  of  the 
apparent  acceleration  on  the  direction  v  with  respect  to  time,  and 
taking  into  account  equations  (1.4.7)  and  (1.5.4),  we  obtain  equality 

MO™  Ka(l)eos<*-f-Ft(t)»lna,  (1.6.5) 

which  indicates  that  quantity  V^ft)  is  in  this  case  a  projection  on 
the  direction  v  of  a  certain  vector  with  components  Vx(t)  and 
along  axes  x  and  y . 

Function  V^(t)  is  the  current  reading  of  the  integrator  of 
acceleration  whose  axis  of  sensitivity  has  the  direction  v.  Integrat 
ing  this  time  function  and  taking  into  account  equation  (1.5.7),  we 
obtain  a  new  relation 


5,(1)  -  5,  (t)  eo«  «  +  S¥  (<)  sin  «, 


in  which  quantity 

5,(0- J  MO* 


(1.6.6) 


(1.6.7) 


can  be  called  the  projection  of  the  vector  of  apparent  path  of  the 
missile  on  direction  v. 

Relations  (1.6.5)  and  (1.6.6)  can  be  given  with  the  help  of 
equation  (1.2.11)  the  following  form: 
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and 


f.(0  -  cos  air,  (o + tr,(oi 
■*•(')  -  «"« +  kSv(l) ]. 


(1.6.8) 

(1.6.9) 


From  the  obtained  relations  It  follows  that  sums  (1.6.1)  and,  con¬ 
sequently,  and  expression  (1.5.8)  for  the  ballistic  function  t(t ) 
can  be  constructed  aboard  the  missile  by  means  of  the  use  of  current 
readings  of  the  single  integrator  of  acceleration  and  of  its  time 
integral.  Really,  on  the  basis  of  relations  (1.6.8)  and  (1.6.9), 
expression  (1.5.6)  is  converted  to  the  form 

•(0  -  “  l*  (0  +  t-f.  (i)i -(*•+**’>- 

—  t* (fit  4- —  t'kft— k .  (1.6.10) 

As  was  already  mentioned  above,  the  duration  of  the  power-flight 
section  is  determined  by  equation  eft)  ■  0.  In  accordance  with 
equations  (1.6.10)  and  (1.2.11),  this  equation,  which  will  sub¬ 
sequently  be  called  ballistic ,  is  reduced  to  the  following,  more 
convenient  form: 

F*(0  +  4r5.(0*htos*+K  +  rt,'na  + 

+4-[**eos«  +  (y*  +  -^-)9ln«].-  (1.6.11) 


§  7 .  Conversion  of  the  Ballistic  Equation  by  Using 
a  New  Rectangular  System  of  Coordinates 

The  ballistic  equation  (1.6.11)  acquires  a  more  transparent  form, 
if  we  introduce  a  new  fixed  system  of  coordinates  Cn  with  the  same 
origin  as  that  for  the  system  xy .  Axis  £  of  this  system  of  coordi¬ 
nates  is  directed  parallel  to  straight  line  v,  i.e.,  at  the  same 
angxe  a  to  the  horizon  (to  the  x  axis)  at  which  the  axis  of  sensitivity 
of  the  integrator  of  accelerations  should  be  located  (Fig.  *0.  Then 
axis  n  will  be  the  antiparallel  to  the  vector  of  velocity  of  the 
missile  at  the  instant  of  its  hitting  of  the  target  with  calculated 
motion  (see  Fig.  1). 


Coordinates  of  the  end  of  the  calculated  power-flight  section 
of  the  trajectory  of  the  missile  in  the  new  system  of  coordinates 
are  designated  by  and  n*.  They  are  connected  with  coordinates 
x *  and  y*  of  the  same  points  in  the  xy  system  by  equations  of  con¬ 
version 

T  —  *‘»ln>+Yeo*«.  (1.7.1) 

Let  us  designate  projections  of  velocity  of  the  missile  on  the  axis 

5  and  n  with  its  calculated  motion  at  the  instant  of  the  switching 

off  of  the  engine  by  v *  and  v*.  They  are  connected  with  projections 

s  n 

of  the  same  velocity  on  the  axes  x  and  y,  with  quantities  v*  and  v*3 

*  d 

by  similar  equations 

+  *••!»•,  »a->  —  »attna4-»^eoea.  (1.7.,) 

For  moving  coordinates  of  the  missile  in  the  system  £n,  which  refers 
to  its  real  motion,  let  us  introduce  designations  Z(t )  and  pi 't;,  for 
projections  of  its  velocity  on  these  axes  -  v^(t)  and  and, 

finally,  for  projections  of  acceleration  -  w^(t)  and  w^(t). 

Since  the  axis  f  and  direction  v  are_parallel ,  then  the  pro¬ 
jection  of  acceleration  of  the  missile  on  the  direction  v  is  the 
equatl on 

•»<<)- (1.7. 
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and,  consequently,  according  to  equations  (1.4.1)  and  (1.6.2),  we 
obtain 

(1.7.4) 

At  the  initial  instant  (t  =  0)  of  the  power-flight  section, 

coordinates  of  the  missile  Z,(t)  and  n (t)  and  also  projections  of 

its  velocity  vF(t)  and  v  (t)  are  equal  to  zero.  Therefore,  succes- 
s  n  2-2 

sively  integrating  expression  (1.7.4)  for  d  E,(t)/dt  with  respect  to 

time,  we  obtain  equations 


IW -*.(*>— J (1.7.5) 


where  V^(t)  and  Sv(t),  as  before,  are  projections  of  the  apparent 
velocity  of  the  missile  and  its  apparent  path  on  direction  v  or, 
which  is  the  same,  on  axis  £.  They  are  connected  with  the  projec¬ 
tion  of  the  apparent  acceleration  i^(t)  by  relations  (1.4.7)  and 
(1.6.7). 

When  U3ing  equations  (1.7*5) »  (1.7.2)  and  (1.7.1)  the  ballistic 
equation  (1.6.11)  for  determination  of  the  instant  of  the  switching 
off  of  the  engine  is  reduced  to  the  following  simple  fo  m:1 

M0  +  7-t(0-«*+7*r.  (1,7.6) 


Not  entering  into  equation  (1.7.6)  is  either  the  coordinate, of 
the  missile  r, (t)  or  projection  of  its  velocity  on  axis  n.  This 

means  that  the  small  distinction  in  the  mentioned  magnitudes  at  the 
end  of  the  power- flight  section  from  their  calculated  values  n*  and 


‘Equation  (1.7.6)  can  be  obtained  directly  from  the  ballistic 
equation  z(t)  *  0,  in  which  function  c(t)  is  taken  in  the  initial 
form  (1.3.2).  In  order  to  be  convinced  of  this,  it  is  enough  to  sub¬ 
stitute  in  expression  (1.3.2)  coefficient  k  by  its  value  following 
from  the  relation  (1.2.11)  and  to  use  equations  (1.7.1)  and  (1.7.2) 
and  also  those  similar  to  them  for  quantities  Z(t)  and  v^(t). 
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correct  to  smallness  of  the  second  order  should  not  affect  the 
flight  range  of  the  missile.  In  essence  this  is  explained  by  the 
fact  that  for  the  control  of  the  flight  range  of  the  missile  it  is 
enough  to  use  only  the  integrator  of  accelerations  whose  axis  of 
sensitivity  is  parallel  to  the  axis  £  or,  which  is  the  same,  as 
direction  v. 

§  3.  Geometric  Examination  of  the  Condition  of 
Hitting  of  the  Missile  on  the 
Assigned  Target 

From  a  geometric  point  of  view  the  ballistic  equation  (1.7.6)  of 
the  previous  section  can  be  given  the  following  interpretation.  Let 
us  assume  that  the  motion  of  the  missile  i3  such  (Fig.  5)  that  at 
the  instant  of  the  switching  off  of  the  engine  it  proves  to  be  at 
point  A1,  shifted  with  respect  to  points  A  -  the  end  of  the  cal¬ 
culated  trajectory  of  the  power-flight  section  on  a  small  segment 
parallel  of  axis  n.  If  its  velocity  at  this  instant  appears  the 
same  as  that  at  the  end  of  the  power-flight  section  of  the  calculated 
motion,  then  the  trajectory  of  the  missile  in  the  free-fllght  section 
will  have  the  same  form  as  that  in  the  calculated  case.  In  order  to 
obtain  this  trajectory,  it  is  enough  to  shift  by  the  magnitude  of 
segment  AA '  the  whole  calculated  trajectory  forward  in  the  directicr. 
of  axis  n  or,  which  is  the  same,  in  the  direction  of  the  tangent  to 
the  trajectory  at  point  C  of  its  intersection  with  axis  x  (see  Fig.  5). 
It  is  obvious  that  the  point  of  intersection  of  the  "shifted"  tra¬ 
jectory  with  the  x  axis  will  be  remote  from  points  C,  i.e.,  from  the 
calculated  points  of  fall  of  the  rocket,  on  a  smallness  of  the 
second  order. 


Let  us  assume  that,  on  the  contrary,  at  the  end  c.f  the  power- 
flight  section  the  rocket  arrives  accurately  at  the  calculated  point 
tut  at  a  velocity  somewhat  distinguished  from  the  calculated  in 
magnitude  and  in  direction  (Fig.  6  and  Fig.  7).  The  form  of  the 
trajectory  proves  to  be  different,  and  in  general  the  rocket  will 
no  longer  hit  the  assigned  target. 


Nevertheless ,  as  follows  from  the  ballistic  equation  in  the 
form  (1.7.6),  if  the  velocity  of  the  missile  will  be  different  from 
the  calculated  by  a  small  vector  parallel  of  axis  n  (Pig»  8),  then 


the  deviation  of  the  rocket  from  the  target  will  also  be  a  magnitude 
of  the  second  order  of  smallness.  Thus,  in  the  case  of  a  nonrotating 
flat  earth  deprived  of  an  atmosphere,  it  is  possible  to  show  such 
coinciding  directions  of  small  vectors  of  the  change  in  velocity  or 

rocket  at  the  end  of  the  power-flight  section,  in 
which  the  change  in  flight  range  of  the  rocket  has 
*  the  second  order  cf  smallness.  Similar  directions 
car.  be  shown  in  general  if  one  considers  the  rotation  of  the  earth 
and  its  noiiopher-lcal  state  in  the  presence  of  an  atmosphere.  However, 
they  prove  tc  he  different,  and  for  the  solution  of  the  problem;  on 


*  ■"  -•  *  J  4.  »  *  4 
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Pig.  8. 


selection  of  the  proper  instant  of  the  switching  off  of  the  engine 
of  one  integrator  of  accelerations  with  the  axis  of  sensitivity  of 
constant  direction,  in  general,  they  are  insufficient.  The  mentioned 
problem  (al3o  correct  to  smallness  of  the  second  order)  Is  most 
simply  solved  with  the  help  of  two  integrators  of  accelerations  (see 
§  3  cf  Chapter  III)  with  additional  integration  of  readings  of  ,ne 
of  them. 


CHAPTER 


-  r  • 
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THEORY  OP  INERTIAL  CONTROL  OP  PLIGHT  RANGE  OP 
BALLISTIC  MISSILES  IN  THE 
GENERAL  FORMULATION 


§  1.  Expression  of  the  Error  in  Flight  Range  of  the 
Missile  with  Small  Changes  in  Parameters  of  the. 

End  of  the  Power-Flight  Section. 

Ballistic  Function 

The  theory  of  inertial  control  becomes  incomparably  more  com¬ 
plex  if  with  the  motion  of  the  missile,  in  contrast  to  the  simplify¬ 
ing  assumptions  accepted  in  Chapter  I,  we  take  account  of  the  change 
Ln  acceleration  of  the  force  of  gravity  both  in  magnitude  and  in 
direction,  consider  the  earth  no  longer  to  be  fixed,  and  in  the  cal¬ 
culation  of  the  free-f light  section  take  account  of  the  resistance 
of  the  atmosphere.  The  basic  difficulty  consists  here  in  the  selec¬ 
tion  of  a  certain  rather  simple  function  of  parameters  measured 
aboard  the  missile  by  Inertial  sensing  elements.  The  function 
should  be  such  that  with  the  achievement  by  it  of  the  earlier 
assigned  value,  it  was  possible  to  produce  a  switching  off  of  the 
engine,  having  provided  the  calculated  range  of  the  missile.  An 
example  of  a  similar  kind  was  given  in  Chapter  I.  This  function, 
cailr’1  nuhr  oqurntly ,  Just  a?  1  r;  Chapter  I,  ballistic  (sometimes  it 
is  called  cuyitrolliny  and  also  controlling  functional) ,  is  con¬ 
structed  aboard  the  missile  by  means  of  the  computer,  which  uses 
current  readings  of  the  integrators  of  accelerations.  The  magnitude 
of  the  ballistic  function  should  be  directly  connected  with  the  error 
in  vhe  flight  range  of  the  rocket  which^would  occur  if  the  switching 
off  of  th<  engine  occurred  it  the  current  instant  of  time.  The 
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system  of  inertial  control  in  turn  should  give  a  signal  for  cessation 
of  the  operation  of  the  engine  upon  achievement  by  the  ballistic 
function  of  a  value  which  corresponds  to  the  turning  of  the  mentioned 
error  into  zero. 

One  of  the  ballistic  functions  is  the  error  itself  in  deter¬ 
mining  the  flight  range  of  the  missile,  expressed  in  terms  of  its 
moving  coordinates  x(t)3  y(t)t  z(t)  and  projections  v  (t),  v  (t), 

*  v 

v  (t)  of  velocity  relative  to  the  so-called  starting  system  of 
coordinates  xyz  connected  with  the  rotating  earth.  In  this  case  it 
is  proposed,  cf  course,  that  at  precisely  the  instant  t  complete 
switching  off  of  the  engine  occurs. 

It  is  convenient  to  put  the  beginning  of  starting  system  at  the 
center  of  the  earth,  having  directed  the  y  axis  along  the  radius  cf 
the  earth  through  the  points  of  start  0  and  the  x  axis  so  that  the 
coordinate  plane  xy,  called  the  program,  contains  the  assigned  points 
of  fail  of  the  missile  C.  If  we  direct  the  x  axis  to  the  side  of 
points  Ct  then  direction  of  the  a  axis  is  thereby  completely  deter¬ 
mined  (Fig.  9).  It  Is  useful  to  note,  that  the  trajectory  of  the 
missile  does  not  lie  in  plane  xy ,  passing,  nevertheless,  in  the  cal¬ 
culated  case  through  the  origin  of  coordinates  0  and  through  point 

For  the  magnitude  of  the  flight  range  of  the  missile  l  we  can 
take  the  shortest  geographical  distance  between  the  point  of  start  '• 
and  actual  point  of  fall  of  the  missile  D,  and  for  its  lateral  devia¬ 
tion  b  -  distance  of  points  D  to  the  programmed  plane,  i.e.,  plane 
(Fig.  9). 

hot  us  introduce  brief  designations 

*(«) 

**(«)-*«..  v* (<*)-»*. 

for  values  of  coordinates  of  the  missile  and  projections  of  its 
velocity  relative  to  the  starting  system  of  coordinates  xyz  at  the 
instant  of  the  switching  off  of  the  engine  t  -  o.  In  this  rase  it 


*(<*)-*. 


Fig.  9. 


is  considered  that  the  instant  t  *  0  is  referred  to  the  beginning 

of  motion  of  the  missile  from  the  start.  With  fixed  values  x,  y ,  z 

and  v  ,  v  ,  v „  range  l  depends  on  the  width  of  the  place  of  the 
x  y  -t 

start  and  on  the  location  of  the  plane  xy  relative  to  the  point  of 

the  compass.  The  position  at  the  instant  of  termination  of  the 
active  section  of  the  celestial  bodies  -  moon  and  sun  relative  to 
the-  earth,  practically  does  not  govern  the  flight  range  of  the 
missile.  A  certain  small  effect  on  the  magnitude  of  the  range  of 
the  missile  proves  to  be  the  state  of  the  atmosphere,  especially 
near  the  points  of  ics  fail,  which  subsequently,  however,  is  not 
considered. 


Having  in  mind  any  specific  position  oi'  start  of  the  missile 
anti  its  assigned  point  of  fall,  we  take  tha  .  range  l  is  only  a  func¬ 
tion  of  coordinates  xt  y ,  s  and  projections  of  velocity  v  ,  v  ,  v 

x  y  z 

of  the  missile  at  the  end  of  the  power-flight  section  of  flight  and 
dot.;  not  depend  on  the  duration  of  the  latter.  Thus,  we  will  Cott¬ 
le  i  or  v  that- 

y,  X,Vx,  Vy,  »,).  (2.1.2) 


The  construction  of  functions  l(x,  y ,  zs  v£>  v  ,  on)  in  general 
in  the  form  of  an  explicit  function  of  their  variables  is  impossible, 
since ,  in  contrast;  to  that  given  in  Chapter  I,  the  eou'aion  >f  motion 
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of  the  missile  in  the  free-flight  section  in  quadratures  are  not 
integrated.  Nevertheless,  mathematical  machines  with  great  accuracy 
allow  comparatively  rapidly  calculating  the  magnitude  of  range  of 
the  missile  l  according  to  data  of  its  coordinates  x,  y,  z  and  pro¬ 
jections  of  velocity  v  ,  v  ,  y  ,  which  corresponds  to  the  end  of  the 

*  y  — 

power-flight  section  in  the  starting  system  of  coordinates  xyz. 

Let  us  designate  by  x*,  y*,  z*  and  v y*,  v*  those  values  of 

x  y  z 

coordinates  and  of  projections  of  velocity  of  the  missile  in  the 
starting  system  of  coordinates  xyz ,  which  corresponds  to  the  end  of 
the  power-flight  section  in  the  calculated  motion  and  by  l*  the  cal¬ 
culated  range  of  its  flight.  In  accordance  with  equation  (2.1.2), 
we  have 

t  - »\ *\  (2.i.3) 


Let  us  form  expression 

Al  mm  1 1* (I),  »(<).*(*),  »«(*).»»(').  *«<0I  — 

■  —  f  (**,  jr*,  *\  »«,  *>i).  (2.1.4) 

It  can  be  accepted  as  the  ballistic  function  mentioned  in  the 
beginning  of  this  section.  Actually,  if  we  turn  off  the  engine  at 
any  arbitrary  instant  t ,  then  value  of  the  selected  function  corre¬ 
sponding  to  this  Instant  is  directly  the  error  in  the  flight  range 
of  the  missile.  If,  however,  the  engine  is  turned  off  at  that 
instant  t  *  a  at  which  this  function  turns  into  zero,  then,  naturally, 
there  will  not  be  an  error  in  the  flight  range. 

because  of  the  absence  of  the  equation  which  expresses  function 

(.2.1.2)  in  an  explicit  form  in  terms  of  its  arguments,  the  practical 

use  of  the  difference  (2.1.4)  as  a  ballistic  function  is  difficult. 

Considerably  simpler  is  another  ballistic  function,  which  is  obtained 

from  the  expansion  of  function  l  into  Taylor  series  about  calculated 

values  of  its  arguments,  i.e.,  about  the  totality  of  quantities  x*, 

y *,  z*,  v*  oA,  and  v*.  According  tc  expression  (2.1.4),  retaining 
xyz 

only  smallness  of  the  first  order,  we  have 


(2.1.5) 


m  -  (*<<)-*' >\1L  + 1* (/)-„•)  £  +  *.  + 

+  mo -*:i  ^ + mo—  »:i  £. . 


Here  the  partial  derivatives  called  ballistic 

coefficients,  are  taken  at  values  of  their  arguments,  respectively 
equal  to  x*,  y*,  z *,  v*,  v*  and  v *,  and,  consequently,  for  the 
assigned  flight  of  the  rocket  are  constant  quantities.  The  ballistic 
efficients  can  be  calculated,  specifically,  on  high-speed  computers. 

It  is  obvious  that  with  a  small  deviation  in  the  motion  of  the 
missile  from  calculated  as  a  ballistic  function  one  can  take  expres¬ 
sion 


4  (I)  -  (x  </)  -  x‘l  £  +  ||f  (»)  -  y‘|  £  +  U  (0  -  «*1  £  + 

+IM0  -  *1  £-+  K(0  -  *!  5T+  MO -*.*]£.  (2.1.6) 

i 

• 

nctually,  with  the  switching  off  of  the  engine  at  the  instant  when 
this  expression  becomes  equal  to  zero,  the  difference  between  the 
real  and  calculated  flight  range  of  the  missile  proves  to  be  a 
magnitude  of  the  second  order  of  smallness.1  In  this  case  deviations 
in  the  parameters  of  the  actual  motion  of  the  missile  from  the  cal¬ 
culated  &»  the  end  of  the  power-flight  section  of  Its  flight  are 
taken  smallness  of  the  first  order. 

§  2.  Equations  Connecting  Coordinates  of  the  Missile 
in  Starting  and  Nonrotating  System  of  Coordinates 

To  o u  1 problems,  vf  inertial  guidance-,'  the  starting  system  of 
'ordinates  represents  certain  inconveniences.  Specifically,  the 
determination  of  coordinates  and  components  of  velocities  of  the 

’If  we  do  not  take  into  account,  of  course  other  facts  disturb¬ 
ing  the  accuracy  of  flight  of  the  missile,  specifically,  the  effect 
- f  atmospheric  conditions  at  the  end  of  the  power-flight  section. 
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missile  in  this  system  according  to  readings  of  newtonmeters  or 
integrators  of  accelerations  is  complicated  by  the  necessity  of 
calculation  of  translational  and  coriolis  forces  of  inertia.  Let 
us  introduce,  therefore,  the  nonrotating  system  of  coordinates  CnC 
with  the  same  beginning  at  the  center  of  the  earth  as  that  for  the 
starting  system  xyz ,  which  rotates  together  with  the  earth.  Let  us 
assume  that  at  the  instant  of  the  beginning  of  motion  of  the  missile 
t  =  0  the  axes  of  system  of  coordinates  and  xyz  respectively 

coincide . 

Let  us  designate  by  U. ,  U  and  Ur  projections  of  angular 
velocity  of  the  earth  on  corresponding  axes  of  the  system  of  co¬ 
ordinates  SnC  and  by  11,11  and  U  —  on  axes  of  the  system  xyz. 

x  y  z 

For  an  arbitrary  instant  of  time  the  following  equalities  are 

valid 


Ut~Ua»lU,  £/,  =  £/,=.  mil. 

Ui-Ut-nU, 

where  l,  n  -  direction  cosines  of  the  earth's  axis  and,  con- 

..eurent,  iy ,  vector  U  in  the  system  of  coordinates  CnC  or,  whi  ur.  is 
me  same,  in  the  system  of  coordinates  xyz. 

After  the  time  t  after  the  start  of  the  missile,  trie  system  cf 
coordinates  xyz  will  be  turned  counterclockwise  at  angle 

?  mm  Ut 

relative  to  the  system  SnS,  if  we  observe  rcmtioi.  or.  the  silo  of 
the  positive  direction  of  the  vector  of  angular  ciocit.v  :•?  the 
earth  V  (i.e.,  on  the  side  of  the  star  Polaris.). 

Together  with  the  current  coordinates  of  the  missile  x(t),  yd) 
and  z(t)  in  the  starting  system  of  coordinates  xyz ,  let  us  introduce 
its  coordinates  Ut)t  r\(t)t  t,(t)  in  the  system  Axes  of  the 

ratter,  as  follows  from  the  aforementioned,  do  not  change  their 
rier.talion  relative  to  directions  at  fixed  stars. 


axes  of  system  of 


The  table  of  cosines  of  angles  between  the 
coordinates  and  xyz1  is  the  following: 


*  (I  — cos  9)/*+ 

+  C06<P 

y  (i  —  cos  V)lm— 
—  a  sin  9 

*  (1  —  eo*9)/n+ 

+  *»sin  9 


n 

(1  —  eos»)«/4- 
4nsin<p 

(1  —  cos  9)  «*  + 
+  cos  9 

(1  —  cos  9)  mn — 
—  /  sin  9 


(1  —  cos«p)n/— J 
—  mslaqi, 

(i  —  eos«p)nm+ 
+  (  sin  9, 

(1  —  cos  9)  a*  4* 

+  «*f. 


(2.2.3) 


According  to  table  (2.2. 3) ,  one  can  determine  coordinates  of 
the  missile  in  the  system  £ns,  i.e.,.  quantities  t )t  r\(t),  and  z,<t) } 
if  direction  cosines  l,  m,  n,  angle  9  and  coordinates  x(t),  y(t), 
z(t)  are  known.  The  appropriate  equations  have  the  form 

E  (I)  “  l(i  —  cos  9)  l*  +  COS9I  *  (I)  4-  1(1  —  cos <f)lm  — 

—  n  aln^l  if  ( I )  +•  f(l  —  eo«9)  In  4.  m  sin 9)  t  (0. 

»l  (/)  *  [(1  —  C0S9)  ml  -f-  n  sln<pl  x  {<)  4.  |(1  —  eo99)  m*  4- 
4-COS9I  j r(0  +  f(l  —  COS9)  mn  —  l  sin <p1  *(<). 

5(0  -*  1(1  —  COS9)  nl  —  m  sin  9!  x  (t)  +  1(1  —  00*9)  mn  + 

4-  /sin9ly  (0  +  ((1  —  C0S9)  «*  +C0S9I  i(0-  (2.2.10 


The  reverse  conversion,  which  expresses  coordinates  x(t).  y(t ) 
and  zltj  in  terms  of  Z(t) ,  n (t)  ana  Z(t) ,  corresponds  to  the  finite 
turn  of  the  system  of  coordinates  £n5  relative  to  the  system  xyz  at 
tii*  same  angle  <fc  but  in  the  opposite  direction.  Specifically,  it 
is  obtained  from  the  previous  equations  (2. -.4)  if  we  replace  in 
them  angle  0  by  and  exchange  places  of  quantities  C  and  x,  n  and 
y ,  and  also  (,  and  z.  As  a  result  we  obtain  equations 

X  (f)  =*  |(1  —  ccct)  I*  -f  1 1  (/)  ;  ((1  —  co«<f)  ml-\ 

+  n  sin?)  r\  (0-1- 1(1  —  C099)  »/-'/»  sin9l  5(0. 
y  ( 1 )  -  ((1—  C0S9)  Im  —  n  sin  9I  l  (0  +  1(1  —  cos*)  m*4* 

4-COS9I  n  (0  4  1(1  —  cosy)  nm  4  l  sin 9 1  5  (0. 
s  (0  “  HI  —  eosf)  (n  4  m  sin 9I  5  (0  4  K*  —  cos f)mn— 

~/sin9l»)(0  4  K1  —  cos  9)  n*  4COS9I  5(0.  (2.2.5) 


u;pe,  cp'-ciflcaliy ,  toe  Appendix  on  page  ?8. 


which,  of  course,  can  be  obtained  directly,  using  table  (2.2.3). 


Projections  of  the  Velocity  of  the  Missile 
Relative  to  the  Nonrotating 
System  of  Coordinates 


Earlier  there  were  introduced  designations  v£(t),  and 

v  (t)  for  current  values  of  projections  on  axes  x,  y  and  z  of  veloc- 
2 

ity  of  the  missile  relative  to  the  starting  system  of  coordinates 
xyz  rotating  together  with  the  earth.  Projections  of  velocity  of 
the  missile  in  nonrotating  system  of  coordinates  £n?  on  the  axis  of 
this  system  are  designated,  respectively,  by  u^(t) ,  nn(t)  and 
It  is  obvious  that 


MO— (2.3.1) 

in  the  same  way  as 

*.<«)-  *.«»>— Tp-  (2.3.2) 

Quantities  vx(t)t  and  vjt)  can  be  expressed  by  u^(t), 

u^(t)  and  inversely.  According  to  equations  of  kinematics  relative 
to  the  motion  of  the  points,  we  have 


■„(*)  -  <0 + 

«,(*)  -M0  +  *.»(*>  -  (0. 

MO-MO  +  tfdrW— 0*«O.  (2.3.3) 


where  ux(t),  u  (t),  ug(t)  -  projections  on  the  axes  x,  y  and  z  of 
the  velocity  of  the  missile  In  the  rotating  system  of  coordinates 
Snc. 

In  equations  (2.3.3),  in  accordance  with  equalities  (2.2.1), 
let  us  replace  V  ,  U  and  Uz ,  and  coordinates  x(t),  y(t)  and  z(t) 
will  be  expressed  by  means  of  relations  (2.2.5)  in  terms  of  co¬ 
ordinates  Z(t),  r\(t)  and  c (t).  If  further  projections  ux(t) , 


Th 

th 

re 
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and  u  (t-  are  expressed,  using  table  (2.2.3),  in  terms  of  projections 
and  u  (t),  then  we  will  arrive  at  the  following  desired 

equations 

n*(«)  =*  1(1  —  cos  9)  I*  +  cos^l  u*  (i)  + 

+  I(i  —  cos  9)  im  -f-  n  sin  9 1  u,  (i)  + 

+  1(1  —  cos  9)  In  —  m  jin  f  I  (»)  —  U%  {1(1  —  cos  9)ni  -f 
+  m  siof  i  {  (1)  +  ((1  —  co*<p)  m/i  —  i  sin  9I  n  (*)  + 

+  1(1  —  cos  9)  n*  +  cos9l  £  (i)}  +  Ui  1 1(1  —  cos  9)  Im  — 

—  n  sin  9  I  5  (i)  +  ((1  —  cos  9)  m*  +  cos  9  I  t)  (1)  + 

+  1(1  -  oo«f)  mn  +  i  da  9I  8  (0). 

*/t) »  |(1  —  cosf)  bn  —  n  sic  9  1 14  (I)  + 

+  1(1  -  cos?)  m*  +  cos  9  1  u,  (1)  + 

+  1(1  —  coif)  ns  +  i  sin  9 !  u«  «- 

—  {((1  —  cos 9)  P  +  cos 9]  l  (I)  + 

+  1(1  —  cos  9)  ml  +  n  sin 9!  n  (0  +  - 

+  1(1  —  cost)  ml  —  m  sin  9!  {  (/)}  +  ! 

+  Vl(l(l—eo*9)/n+mslaflJ(0+((l—co6f)mn— 

—  i  sU  »I  tj  (1)  +  ((1  —  cos 9)  n*  +  cos 9 1  C  to), 

/»,<*)  =  1(1  —  cos?)  In  +  m  S109]  «n  (<)  + 

+  1(1  —  cca 9)  mn  —  l  sin  9!  u,  (0  +1(1  —  coef)  n *  + 

+  cos 9 1  Ji*(i)  —  I/t  { 1(1  —  COS9)  Im  —  it  sin 9I  i  (<)  + 

+  1(1  —  cos 9)  m*  +  C0S9I  i\  (i)  + 

+  ((1  —  cos  9)  mn  +  i  3U9I  C  (t)}  + 

+tf»  (1(1  —  c<»9)  I*  +  eo»9l  t  {<)  + 

+  1(1  —  cos 4)  ml  +  n  si  a  9]  f|  (<)  + 

+  1(1  —  cos  4)  ni  —  m  SID9I  $  (<)}.  (2.3.14) 

The  last  equations  are  somewhat  simplified  after  the  replacement  in 
them  of  projections  of  angular  velocity  of  th<  earth  U and  U 
respectively  by  products  W,  mV  and  nU  in  accordance  with  the  same 
equalities  (2.2.1).  As  a  result  we  obtain 

*.(<)  -  ((1  —  e»9)  P  +  eoef]  H  (0  + 

+  1(1  —  COS9)  mi  +  n  sJd9F  •*,  (I)  + 

+  1(1  —  cos 9)  nl  —  m  sin  9!  (f)— 

—  U  (cos 9  (m$  (<)  —  irr)  (<)|  — 

-  i  sin  9  1/5  (<)  +  mn  (0  +  n  C  (i) I  +  l  (/)  s/n  9). 


*VW“  f(l  —  eos  9)  Im  —  n  sin  ip  >  at  (/)  + 

4- 1(1  —  eos^p)  m*  4-  eos<pl  u,  ( t )  4- 
4- 1(1  —  coflf)  mn  4-  /  sin  9 1  (<)— 

—  U  (cos 9  In£  (/)  —  /£  (<)I  —  m  sin<p[J!  (t)  4-  m*j  (<)  + 

+  "#(<)!  +*I  (0  s»n<p}, 

v/t)=-  ((1  —  cos  9)  /«  +  m  sin  9)  m  (<)  4- 
+  1(1  —  eosf)  mn  —  l  sin 9]  u*  ( t )  4- 
4.  ((1  —  C0S9)  n*  4-  eos 9I  U*  (0  — 

—  U  {cos  9  I/ij  (!)  —  m£  (01  — 

-  n  sin  9  US  (0  +  wt|  (0  4-  «6  (01  +  t  (0  (2.3.5* 


It  is  obvious  that  by  similar  means  it  is  possible  to  arrive  at 
equations  which  express  quantities  ,  u ^  and  u ^  in  terms  of  v x, 

v  ana  v  .  We  have 

y  2 


«H(*)-  1(1  —  eos 9)  i*  4-  cos 9]  uM  ( t )  4- 
4-  1(1  —  cos 9)  tm  —  n  8(09)  vv  ( t )  4- 

4.  ((1  —  COS9)  In  4-  m  sin9l  v,  (0  4- 

4-  V  {C039  fmz  (0  —  »y  (01  4- 1  s*n«P  U*  (0  + 

4-  my  (t)  +  ns  tyf—  x  (t)  sin  9}, 

«,(/)=*  1(1  —  C0S9)  ml  4-  n  sin9)  v«  (0  + 

4-  1(1  —  C0S9)  m*  4-  COS9I  vv  (/)  4- 

4-  1(1  —  CO69)  mn—  l  sin9l  v,  ( t )  4-  U  {cos 9^/11  (/)  — 

—  U  (01  4-  m  sin  9  [lx  (t)  f  my  (/)  4-  nz  (/>!  -  y  (0aii»9), 

u c(f)=*  1(1  —  cos9>  nl  —  m  3109]  vx  ( t )  4- 
4-  ((1  —  COS9)  nm  4-  /  5109!  vv  ( t )  4- 
4-  [(1  —  C0S9)  n*  4-  cos<pl  v»  (0  4* 

4-  U  {cos  9  Uy  (l)  —  mx  (01  + 

4-  n  sin  9  U*  (0  +  m»  (0  4-  nz  (01  —  z  (/)  *109}.  (l.  5. 0 


Error  in  Flight  Range  of  the  iiissile  as  • 
1  unction  of  Changes  'in  Parameters  of  the  Ena 
of  the  powered- F ! i ght  Sent! on  in  the 
" Nonr-utatlng  System  of  Coordinates. 
Initial  Ballistic  Function 


Wo  will  designate  coordinates  and  projections  of  velocities  o 
i h"  missile  in  the  nonrotating  system  of  coordinates  CnC  at  the 
instant  of  the  end  of  the  powered-flight  section  of  its  flight 
■  -- -  ? ,  respectively,  by  letters  Z ,  <) ,  C  ana  ur,  ur  and  u  (  .  They 


are  similar  to  designations  (2.1.1)  introduced  earlier  for  quantities 
xt  y,  2  and  v£l  vy,  yg,  which  refer  to  the  starting  system  xyz.  Thus, 

60*)- 5.  1(a)  -n.  C0*)«l. 

«»(a)«-u„  (2.4.1) 


Setting  in  equations  (2.2.5)  and  (2.3.4)  4s  »  Uo ,  we  obtain  expres¬ 
sions  of  quantities  x,  y ,  a  in  terms  of  £,  n,  z  and  o  and  also 
quantities  v  .  u  ,  u  in  terms  of  u,,  u_,  u_,  5,  n ,  Z  and  o.  This 

xyz  s  n  4 

allows  considering  magnitude  of  range  Z  as  a  function  of  values  of 
coordinates  of  the  missile  i,(t) ,  n (t)  and  z(t)  at  the  instant  t  *  o. 
projections  of  its  velocity  u^(t),  u^(t)  and  u.(t)  in  the  nonrotational 
system  of  coordinates  SnC  (at  this  instant  of  time)  and,  finally,  the 
duration  itsel^  of  the  powered-flight  section  o.  In  accordance  with 
formula  (2.1.2)  we  now  have 

(2.4.2) 

It  is  easy  to  explain  why  In  this  case  the  range  l  clearly  depends  on 
variable  o.  Actually,  the  position  of  the  rocket  relative  to  the 
earth  at  the  same  values  of  its  coordinates  in  the  ncnrotating  system 
substantially  depends  on  the  position  of  the  latter  with  respect 
to  starting  system  of  coordinates  xyz  connected  to  the  earth,  i.e., 
on  angle  =  fo .  The  same  refers  to  magnitudes  of  projections  of 
the  vector  of  velocity  of  the  missile  in  the  starting  system  at  the 
instant  of  termination  of  the  powered- f .tight  section. 

Calculated  values  of  parameters  of  the  end  of  the  powered- 
flight  section  in  the  starting  system  of  coordinates  xA,  y  * ,  z*f 
y<,  y *,  and  v *  correspond  to  calculation  values  of  parameters  £* , 
n* ,  ut,  a*,  and  u *  in  the  nonrotating  system  £n5.  They  are 

connected  with  each  other  by  equations  (2,2.5)  and  (2.3.4)  or  (2.2.4) 
and  (2.3.6  ,  in  which  it  follows  to  assume  4  =  yc»,  considering  cr* 
as  the  designation  of  the  calculated  duration  of  the  powered-flight 
sootier  ,'f  flight  of  tne  missile.  It  is  obvious  that  in  accordance 


/•  mg  •  •  •  •  •  *.  .  ,*•  •  »•  •  •  •  •. 

-  /(* .  y .  * .  Py.  9,)  -  HI ,  t| ,  C . »».  at;  « ), 


(2.4..  3) 


where  l*,  as  previously,  is  the  rated  value  of  the  range. 

The  difference  LI  between  the  magnitude  of  the  actual  range  of 
the  missile  l  and  its  calculated  value  l *  can  be  represented,  similar 
to  equation  (2.1.5),  in  the  form  of  the  expansion 

w  -  a  -  D -J+ (»1  - 11*)  (C  -  D +  («t  -  ^  + 

(2.4.<o 

where  terms  of  the  second  and  higher  order  relative  to  differences 
are  dropped 

I— 1\  n — n.  C — C*.  h — *1.  *% — ■«  *c — o — 9. 


In  contrast  to  the  expansion  (2.1.5) >  entering  into  equation 
(2.4.4)  is  a  term  of  the  first  order,  which  contains  by  a  factor 
the  difference  a  -  d*  between  the  actual  time  of  the  powered-flight 
section  of  flight  of  the  missile  a  and  its  calculated  value  a*.  The 
same  difference  is  contained  by  a  number  of  terms  of  expansion 
(2.4.4)  of  a  higher  order.  The  reason  for  this  fact  was  already 
explained  in  the  beginning  of  this  section. 


Derivatives 


91  9t  M  M 

1C’  "J* *  '  S3J’ 


are  themselves  functions 


of  variables  £,  n,  C,  u ^  and  a.  In  the  expansion  of  (2.4.<i) 

they  should  be  taken  at  calculated  values  of  enumerated  variables 
£*,  n*,  C*»  u|,  u*,  o#,  and,  consequently,  for  the  specific 

assigned  flight  of  the  missile  arc  constant  qu-.nti tics  -  ballistJ 
coefficient,  which  refer  to  the  nonrotating  system  of  coordinates. 

Let  us  show  that  they  are  all  expressed  in  terms  of  ballistic 
coefficients  of  the  starting  system  -Jj- ,  *|j-,  jj~,  and 

Actually,  according  to  rules  of  differentiation  of  complex  functions, 
specifically,  we  have: 
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(2.4.5) 


at  at  »x  .  at  a*  ,  at  it  ,  a'®*,  a/  ®»,  a/  *», 
1C  "  "ST  ^  ^  ®v  +  IT  ®r  +  %  aT  +  ®T +5^  af 


Analogous  expressions  take  place  for  derivatives 
and  for  the  partial  derivative  -Jj-  equation 


at  tai  at  at 

*T*  ic- 


a/  ,a/  a»  x,a/  ay  .at  'a*-,  a/  •*«  ,  at  ,  *L*», 
laT  *“-l*  a«  as  •  •  at  "ST"*’  5?jj  as  a«  a*,  as  ' 


(2.4.6) 


is  valid.  Entering  into  enumerated  expressions,  the  partial  deriva¬ 
tives  of  variables  x,  w,  z,  v  ,  v  v  with  respect  to  variables  £, 

x  y  z 

a ,  c,  u,,  o  ,  y  and  a  are  formed  by  means  of  the  application  of 

S  M  S 

appropriate  operations  of  differential  equations  of  the  type  (2.2.5) 
and  (2.3.5),  in  which  it  follows  preliminarily  to  set  t  *  o  and 
i,  =  Ua . 


Thus,  for  instance, 

■  j-  •  co»W+(l  co*  U9)  l*, 

A-  || 

«•  (1  —  cos  Ua)  Im  +  n  sin  Ua , 

^  »  — V  (1  —  /*)  sin  Ua, 

gjL  mmU  | — (1  — 7*)(sin Ua  4-  ( Im  sint/o  4-neos  t7a)jj  4. 
+  (to  sin  Ua  —  moot  £to)Cl. 

-gf  ■■  V  t — (1  — -  P)  ii(Sln  Ua'  +  (Im  sin  Us  4- »  cos  Ua)  a;  + 
+  (to  sin  Ua  —  moos  Ua)  u5J  4-  U*  ((mg  —  n»j)  sin  Ua  4- 
4- 1  (/£  4-  m  tj  4-  ng)  cos  Ua  —  \  cos  Us  j . 


If  now  in  equations  (2.4.5),  (2.4.6),  (2.4.7)  and  ones  similar 
to  them,  we  assume  the  values  of  all  variables  equal  to  the  cal¬ 
culated  values,  then  thu3  the  desired  coefficients  in  expression 
(2.4.4)  will  he  completely  determined. 


,jl-  wllJ  assume  that  as  a  resuit  of  the  small  deviation  in  the 
r-oai  me.  t ion  of  the  missile  from  the  calculated  in  the  expansion  of 
11 ),  the  terms  caused  by  the  so-called  lateral  deviation  in  the 

missile  (8  —  8*)  -Jg-  and  («c  —  “t)^»  are  small  in  comparison  with  other 

terms.  \s  a  result  the  difference  between  the  actual  and  calculated 
.Lght  range  of  the  missile  LI  can  be  represented  by  the  approximate 
e  i  nation 


<i  -v>^+ <»*  -  »;i  £ + 


g .  U  .  S ' 


91  31  91  dt  91  ,  , 

whore  -fir,  and  -55-.  as  was  already  snown  above,  are  number:- 

which'  are  completely  defined  by  the  selected  calculated  motion  cf 
the  missile. 


expression  lc.4.b)  corresponds  to  the  following  balii: 
motion,  which  we  call  the  initial  ballistic  function : 


•(0*  18(0— ri-Jf  +  !n(0— +  l»e(0  — »Jl^+ 


.l-re  r,  t  .1  .  n*(  t)  and  u  (t),  as  before ,  are  moving  coordinate’s 

*  0  n 

mil  project !  one  of  velocity  of  the  missile  in  the  nonrotatir.g  system 
•  f  coor-  ilnates  £,nb. 


if  the  engine  of  the  missiie  is  turned  off  at  any  .u.  tant  •  -  1 , 
tr.vn  the  bn.'*;'  of  equ^t  •  -  nr  (  ?  .  4  .  3  )  and  (P.4.8)  t!v  "t  1  u-  of  ■ 
ballistic  function  (2.4.9)  corresponding  ti>  this  instant  coiwect  t. 
smaianess  of  the  second  order  determines  the  error  in  tne  flight 
range  of  tne  missile.  To  avoid  this  error,  one  snould  turn  off  the 
.-agine  ...f  the  missile  at  the  instant  when  the  ballistic  function*' 

•  * :  ’urns  Into  aero.  Thus,  tne  necessary  instant  of  the  switching 

f f  f  • ne  engine  is  determined  by  the  root  f  equation 


(2.4.10) 


15(0-  r  I  -Jr  +  ln<0  -  »TI  ^  + 1 MO  -  «t*l  ^  + 

+  [MO  -  «;i  ^  +  (<  -V)  •£■  -  o. 

Equation  (2.4.10)  will  be  called  the  initial  ballistic  equation. 

The  basic  content  of  subsequent  sections  of  this  chapter  anvi 
also  Chapters  III  and  IV  consists  in  the  construction  of  a  series 
of  other  ballistic  functions,  with  the  help  of  which  it  is  possible 
as  simply  as  possible  to  solve  the  problem  about  the  inertial  control 
of  the  flight  range  of  the  ballistic  missiles.  They  all  are  obtained 
by  means  of  certain  conversions  from  function  (2.4.9)  and  different 
presentations  of  coordinates  and  projections  of  velocity  of  the 
missile  in  the  nonrotating  system  £n 5  in  terms  of  current  readings 
of  integrators  of  accelerations. 

§  5 .  Differential  Equations  Which  Determine  Current 
Coordinates  of  Motion  of  the  Missile 


The  most  natural,  but,  as  it  appears,  not  the  simplest,  method 
; f  construction  of  the  left  side  of  the  initial  ballistic  equation 
(1.4.10)  aboard  the  missile  Is  the  preliminary  obtaining  of  most 
moving  coordinates  t,(t)  and  n  (tJ  by  means  of  solution  of  the  dif¬ 
ferential  equations,  which  connect  these  coordinates  with  current 
readings  of  two  newtonmeters .  The  axes  of  sensitivity  of  the  latter 
should  maintain  the  fixed  directions,  respectively,  parallel  to  axes 
£  and  n  of  the  nonrotating  system  of  coordir  -tes  £nc.  Readings  of 
the  mentioned  newtonmeters  are  projections  c.it)  and  of  the 

apparent  acceleration  of  the  missile  on  the  axis  £  and  n.  According 
to  equation  (1.4.1)  of  Chapter  I,  in  this  case  we  have 


wr^ 


in  «=  *5(0- 
•M0  * 


MO  “M0 — /«. 


MO 


i it 


(2.5.1) 

(2.5*2) 


;  e  c  t  i 


on  axes  £  and  n  of  the  real  acceleration  of  the 
t->  the  fixed  system  of  coordinates  £nc,  and 


a  t  i  ve 


A-AU(0.n(*).C(');<l. 
/•-/•U(0.« 1(0.  CW  M 


(2.5*3) 

( 
l 
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are,  respectively,  projections  of  acceleration  of  the  force  of  the 
earth’s  gravity,  which  depend  in  evident  form  not  only  on  coordinate 
of  the  missile,  but  also  on  time  t.  For  an  explanation  of  the  last 
fact  let  us  note  that  even  at  the  point  with  constant  coordinates 
relative  to  the  nonrotating  system  the  acceleration  of  the  force 
of  gravity  is  changed  with  the  course  of  time.  Actually,  with 
respect  to  the  system  of  coordinates  ShC  the  position  of  the  earth 
is  continuously  changed.  At  the  same  time  the  field  of  gravity  of 
the  earth  does  not  have  radial  symmetry,  specifically,  as  a  result 
of  the  deviation  in  its  form  from  a  sphere.  If  the  last  fact  Is 
disregarded  and  the  earth  is  considered  a  sphere  with  radial  and 
symmetric  distribution  of  density,  then  projections  of  acceleration 
of  the  force  of  the  earth’s  gravity  will  be  represented  by  equations 


no  longer  dependent  in  evident  form  on  time.  Into  them  is  the 
value  of  acceleration  of  the  force  of  gravity  on  the  earth's  surface 
examined  as  a  sphere  of  radius  R,  and 

p-Zurtr+.iiwi’+icwi*  (2.5.5) 


is  the  distance  between  the  missile  and  center  of  the  earth.  Becaus< 
of  the  smallness  of  the  lateral  deviation  in  the  rocket  eft’  In  com¬ 
parison  with  the  radius  of  the  earth  R,  this  coordinate  in  equation 
(2.5.5)  can  be  cropped.  Correspondingly,  it  is  possible  to  drop 
coordinate  t. (t)  in  equations  (2.5*3),  considering  the  effect  of  this 
coordinate  on  the  magnitude  of  acceleration  of  the  force  of  gravity 
unimportant.  Then  relations  (2.5*1)  can  be  considered  as  a  combina¬ 
tion  of  two  differential  equations 


^Sr-M0  +  All(0.n(0;  <1. 
-«,(<)  +  /»  16(0. 1(0;  <! 


(2.5.6) 
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for  the  search  of  desired  functions  i,(t)  and  n (t)  on  the  given 
current  readings  a^(t)  and  of  two  newtonmeters . 

The  direct  solution  of  the  system  of  nonlinear  differential 
equations  (2.5.6)  aboard  the  missile  for  the  purpose  of  seeking  the 
current  values  of  coordinates  Z(t)  and  n (t)  it  is  possible  by  means 
of  the  use  of  a  rather  high-speed  computer.  In  the  subsequent  two 
chapters  the  approximation  solution  to  these  equations,  which  in 
principle  can  be  used  in  the  system  of  control  of  the  flight  range 
of  the  ballistic  missile  in  the  presence  of  onboard  of  its  newton¬ 
meters  or  integrators  of  accelerations  and  simplest  computers  is 
used . 

§  6.  Auxiliary  Relation  Connecting  Magnitudes  of 
Ballistic  Coefficients  in  the  Nonrotating 
System  of  Coordinates" 

At  different  conversions  of  the  initial  ballistic  equation 
(2.,l.l0)  one  relation  proves  to  be  useful,  and  It  connects  magnitudes 
of  the  ballistic  coefficients  •£.  35^  ^  and  £  in  the  nonrotating 

system  of  coordinates  with  projections  of  velocity  u*  and 

with  projections  of  acceleration  of  the  force  of  gravity  n*; 

a*),  f  (5*,  n*;  c*),  which  refer  to  the  instant  t  =  ot  of  the  termina¬ 
tion  of  the  powered-flight  section  of  the  calculated  motion  of  the 
missile. 

Let  us  introduce  functions  Z*(t)  and  n *(t),  which  are  current 
coordinates  of  the  missile  in  its  calculated  motion  on  the  powered- 
flight  section.  They  satisfy  the  totality  of  differential  equations 
(2.5.6),  if  In  them  the  current  values  of  projections  a^(t)  and  a^(t) 
of  the  apparent  acceleration  are  replaced  by  their  calculated  values 
a£(t)  and  a U t) .  Thus, 


«;<«)  +  /*  IV  <0.  u*(0;M. 

»«;<£)  +  /,  ino. 
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(2.6.1) 


Correspondingly,  let  us  introduce  functions  E,**(t)  and  r\**(t),  which 
are  coordinates  of  the  missile  in  the  free-flight  section  of  its 
calculated  motion.  Let  us  assume  that  in  the  beginning  of  motion 
of  the  missile  in  the  free-flight  section  the  resistance  of  the 
atmosphere  is  unimportant,  and,  consequently,  the  apparent  accelera¬ 
tion  is  absent.  In  this  case  one  should  consider  that  functions 
K**(t)  and  n **(t)  satisfy  the  totality  of  differential  equations, 
which  is  obtained  from  (2.6.1)  if  in  equations  of  the  latter  we  drop 
quantities  a*(tJ  and  a*(tj  and  replace  variables  £*(t)  and  n*(t) 
respectively  by  iA*(t)  and  We  obtain 


*»“<*>;  <J. 

n. 


(2.6.2) 


Let  us  designate  the  current  values  of  projections  on  the  axes  £  and 
n  of  the  velocity  of  the  missile  in  the  powered-flight  section  with 
its  calculated  motion  by  u£(t)  and  u*(t) ,  and  by  u£*(t)  and  u^*(t)  - 
the  corresponding  quantities  which  refer  to  the  free-flight  section. 
We  have  the  evident  equalities 


«i(0— (2.6.3) 

and,  similar  to  them 

(2.6.4) 

The  powered- flight  section  of  the  calculated  motion  of  the  missile  at 
the  instant  t  -  a*  will  pass  over  into  the  free-flight  section,  whence 
it  follows  that 


and  also 


!*(«*)  -r  (0-1%  »!•<«•) -n~  («*>-»,%  (2.6.5) 

“i  («*)  -  «r  (<*’)  -  B?.  «;  (a*)  -  u~  (a*) .  (2.6.6) 
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According  to  equations  (2.6.1)  and  equations  (2.6.5)  for  the  end  of 
the  calculated  powered-flight  section,1  we  obtain 


>  “  ae(9  )  +  /*(6\  «"). 

(2.6.7) 

-.;<*•)  +  /,<  e*.  n*:  o*). 


Similar  relations  can  be  ootained  on  the  basis  of  equations 
(2.6.2)  and  the  same  equations  (2.6.5)  for  the  beginning  of  the 
calculated  free-flight  section.  They  have  the  following  form; 


-  ft  (l\  V;  O.  -  /-  (6*.  v,*;  «•).  (2.6.8) 

* 

Let  us  examine  now  on  the  calculated  free-flight  section  of 
flight  of  the  missile  the  certain  instant  t  -  tJt  which  follows 
after  the  instant  of  the  switching  off  of  the  engine  t  =  c»*.  Correct 
tc  smallness  of  the  second  order  we  have  expansions  in  Taylor  series 
for  coordinates  of  the  missile  and  n^,  which  refer  to  this  instant, 
namely : 


-  r  (o +('»-  «*>  . 

*  -  «r  (*.)  -  n~  («•>+<<»-  «*)  . 


(2.6.9) 


or,  taking  account  equalities  (2.6.4),  (2.6.5)  and  (2.6.6), 

-«*)“;•  (2.6.10) 

Similarly  for  projections  u ^  and  of  the  calculated  velocity  of 
the  missile  at  the  mentioned  instant  t  =  wc  have 


'in  designations  of  the  type  dZ,{a)/dt  it  follows,  of  course,  to 
consider,  that,  from  the  beginning  an  operation  of  differentiation  of 
the  appropriate  function  is  produced,  and  then  the  value  shown  in 
parenthesis  ic  giver;  to  the  argument. 


(2.6.21) 


-i-«r  <  v  -  k  («•> + <n  -  . 

•;  -  K  <*»)  -  OT  +  <*,  -  «*)  • 


On  the  basis  of  equations  (2. 6.  it),  It  13  possible  here  to  produce 
the  following  replacement: 


K>_.Ay ,  ^ 


(2.6.12) 


Taking  into  account,  furthermore,  equality  (2.6.6),  and  also  rela¬ 
tions  (2.6.8),  we  obtain 


■i-»;+(<i-o/,(r.  v.  <j*) 


(2.6.13) 


Let  us  imagine  such  a  motion  of  the  missile  in  -he  powered- 
flight  section  different  from  the  calculated,  as  a  result  of  which 
at  the  Instant  t  *  1 2  it  proves  to  be  at  the  point  with  coordinate? 

C-  and  n.  ,  possessing  a  velocity  the  projections  of  which  or.  axes 

1  -*■  ii 

6  and  n  are  respectively  equal  to  quantities  u,  and  u  .  If  at  this 

s  n 

instant  t  =  the  engine  of  the  missile  Is  turned  off,  then  the 
following  free-flight  section  of  its  motion,  because  of  the  unique¬ 
ness  of  solution  to  the  problem  of  dynamics  completely  coincides 
with  the  calculated.  Actually,  both  motions  are  such  that  at  the 
same  instant  t  **  t  in  them  positions  and  velocities  of  the  missile 
relative  to  the  system  of  coordinates  £nc  coincide.  Furthermore, 
when  t  >  t1  they  are  subordinated  to  the  same  totality  of  differential 
equations,  which  refers  to  the  free-flight  section  of  its  motion. 

Thus,  at  such  a  powered-flight  section  the  rocket  will  not  nave  an 
error  in  the  range  of  its  flight.  Consequently,  parameters  of  the 
end  of  the  powered-flight  section  of  its  motion  £ ^ and 
should  satisfy  the  initial  ballistic  equation  (2.  lJ.10),  i.e.. 


+  <•«- »*0  ^ 0. 


( 2 . 6 . 1 4 ) 
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Replacing  here  quantities  and  by  their  expressions  according 

to  equations  (2.6.10),  and  and  in  accordance  with  equations 
(2.6.13),  we  obtain  after, reduction  to  the  common  factor  t1  -  o# 
equality 

*»■*[  +  +  Ml’.  n';  «*)*-  + 

+/«(F.  n*;  «*)^  +  ‘R‘“0,  (2.6.15) 

U  K  Mu  k 

which  connects  the  ballistic  coefficients  -jj-.  and  in  the 

nonrotating  system  of  coordinates  tor,  with  projections  of  accelera¬ 
tion  of  the  force  of  gravity  /» (E*.  V;  «*).  /»  (6*.  n‘;  ®*l  and  with  projec¬ 
tions  of  velocity  of  the  missile  u|  and  u*  at  the  calculated  instant 
cf  the  switching  off  of  the  engine. 

Let  us  exclude  from  equality  (2.6.15)  projections  of  accelera¬ 
tion  of  the  force  of  gravity  by  means  of  relations  (2.6.7)  and  use 
equations  (2.6.3)-  As  a  result  let  us  obtain  the  new  equality 

+[4!gfl_<lo]^+*-«.  (2.6.16) 

which  will  be  used  in  the  following  section  during  the  conversion  of 
the  expression  of  the  initial  ballistic  equation  (2.4. 10)  tc  the  form 
convenient  for  applications. 


§  7 .  Isochronal  Variations  of  Coordinates  and 
Projections  of  the  Velocity  of  the  Missile. 
Basic  Ballistic  Equation 


With  the  juxtaposition  of  the  real  and  calculated  motions  of 
the  missiles,  let  us  call  the  isochronal  variations  of  its  co¬ 
ordinates  the  differences 


dl(i)-E(0-F(«).  4n(0-n(0  —  n*(0-  (2.7.1) 
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Similarly,  let  us  introduce  isochronal  variations  of  projections  of 
velocity  of  the  missile 


K(') -MO—  «»(*>.  WO  «  MO— “5(0  (2.7.2) 

and  projections  of  the  apparent  acceleration 

fat  4(f).  *.,-<.,(0— «;<f>.  (2.7.3) 

We  will  consider  the  enumerated  isochronal  variations  small  quantities, 
the  squares  and  products  of  which  can  be  neglected.  Thereby,  it  is 
proposed  that  the  problem  of  thrust  of  the  engine  both  in  direction 
and  in  magnitude  is  produced  to  a  sufficient  degree  accurately. 


Functions  V(t),  n *(t)3  u}(t),  a{(t)  and  a^(t)  refer  t0 

the  powered-flight  section  of  the  calculated  motion  of  the  missile 

and,  consequently,  are  determined  at  values  of  the  argument  t  not 

exceeding  the  duration  of  this  section  cr*.  At  the  same  time  in 

equations  (2.7.1),  (2.7.2)  and  (2.7-3)  functions  Z(t),  nft), 

u  ( t ),  ar(t)  and  a  (t)  refer  to  the  powered-flight  section  of  the 

n  4  n 

real  motion  of  the  missile,  the  duration  o  of  which  can  be  both  less 
and  more  than  the  calculated  value  a*.  Therefore,  for  the  complete 
certainty  of  variations  £>£(£),  6r\(t ),  6a^(tJ  and 

<$ a  (t)  in  the  whole  time  interval  of  the  powered-flight  section  of 
the  real  motion  of  the  missile,  one  should  agree  upon  what  is  under¬ 
stood  by  functions  Z*(t),  n *(t),  u*(t)}  u*(t),  a*(t ),  a*(t)  with 
argument  t  somewhat  exceeding  quantity  o*  (within  limits  of  the 
allowed  variance  of  duration  of  the  powered-flight  section  for  the 


specific  type  of  missile).  If  one  assumes  that  the  thrust  of  the 
engine  with  termination  of  the  powered-flight  section  is  lowered 
gradually  or  by  stages,  then  for  these  functions  when  i  >  a*  It  is 
natural  to  take,  respectively,  functions  Z**(t)J  r\**(t),  and 


which  refer  to  the  calculated  free-flight  section  and  assume 
n 

functions  and  a*(£)  equal  to  zero.  In  this  case  it  follows  to 

expect  that  the  isochronal  variations  (2.7.1),  (2.7.2)  and  (2.7.3) 


will  be  small  at  the  short  time  interval  directly  following  the 


calculated  instant,  of  the  switching  off  of  the  engine. 


For  missiles  with  the  sudden  switching  off  of  the  engine,  on 
the  contrary,  it  is  more  preferable  to  consider  functions 
r\*(t) ,  u£(t)  and  u*(t)  when  t  >  a*  equal  to  their  calculated  values 
calculated  on  the  assumption  that  the-  engine  of  the  missile  at  the 
instant  t  =  o#  was  not  switched  off.  Here  the  isochronal  variations 
(2.7.1),  (2.7.2)  and  (2.7.3)  again  can  be  considered  small. 

Let  us  use  now  equations  (2.7.1)  and  (2.7.2)  for  Isochronal 
variations  of  coordinates  and  projections  of  the  velocity  in  order 
to  transform  the  initial  ballistic  equation  (2.4.10)  by  which,  as 
was  already  mentioned,  the  necessary  time  of  cessation  of  operation 
of  the  engine  of  the  missile  for  its  hitting  of  the  assigned  target 
is  determined. 


Let  us  present  in  the  beginning  equation-  (2.4.10)  in  the 
following  manner: 

l«0-  V  (01  %  +  ln<0  -  tf  (01^ + (0 -  «;(')! 

+[«,(»)- «;  (oi  ~  -  {it*  (o  -ri-jf + 

+ iv  (o  -  If + K<o  -  “»*i  £ + 

+i®;<o-»;i^+(*-a*)-!y.  (2.7.4) 


Into  the  left  side  of  the  last  equation  let  us  substitute  expressions 
for  isochronal  variations  of  coordinates  and  projections  of  the 
velocity  of  the  rocket,  according  to  equations  (2.6.5),  (2.7.1)  and 
(2.7.2).  In  its  right  side  let  us  expand  each  of  the  functions 

U*rt; ,  u‘(t;  in  Taylor  series  near  the  value  t  =  s*,  retain- 
Jrg  ir.  the  expansion  only  terms  of  the  first  order  of  smallness.  V/c 
obtain 


*<o-|+^(o£ +  «»!«,£, +»“•<'>£, - 

,  *V(r)  dt  ,  dt  ,  an 

'  37*  du^  '  '  dt 
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(2.7.5) 


The  right  side  of  the  last  equation  can  be  considerably 
simplified  if  we  take  into  account  equality  (!’. 6. 16)  of  the 
previous  section.  As  a  result  we  arrive  at  the  following  equa¬ 
tion: 

*1  («)  J  4-  +  6*M0  ^  (0  ^  - 

[<<«*)£  +•;<«*>  £].  (2.7.6) 

which,  subsequently,  we  Will  call  the  basic  ballistic  aquation. 

For  the  actual  use  of  equation  (2.7.6)  in  the  system  of  inertial 
control  of  the  flight  range  of  ballistic  missiles  its  left  side  should 
be  constructed  aboard  the  missile  at  current  readings  of  newtonmeters 
in  the  form  of  a  certain  electrical  or  mechanical  magnitude. 


CHAPTER  III 


APPROXIMATE  PRESENTATIONS  OF  THE  BALLISTIC  EQUATION- 
DETERMINING  THE  INSTANT  OF  SWITCHING 
OFF  OF  THE  ENGINE 


Differentia}.  Equations  for  Isochronal 
Variations  of  Coordinates  of  the 
Missile"  In  the  ftowered-Fllght 
Section  of  Its  Motion 


In  the  previous  chapter  there  was  obtained  the  so-called  basic 
ballistic-  equati  r.  12.7*6),  which  should  satisfy  parameters  of  the 
end  of  the  powered- flight  section  In  order  that  the  rocket  would 
not  have  an  error  In  the  range  of  its  flight.  For  the  construction 
■of  the  left  side  of  this  equation,  it  is  necessary  to  know  aboard 
the  missile  the  Isochronal  variations  of  coordinates  6Crtl  and  6n(t), 
and  also  projections  of  the  velocity  t)  and  The  latter 

are,  of  course,  time  derivatives  of  variations  of  coordinates,  i.e.. 


(3.1.1) 


The  coordinates  themselves  i(t)  and  nft)  of  the  real  motion  of 
the  missile  in  the  powered-flight  section  satisfy  the  totality  of 
di f  l ercnt  i  al  equation^.  ('“.5.6),  and  coordinates  of  the  calculated 
motion  -  totality  (2.6.1).  Let  us  form  the  differences,  respectively 
of  the  left  and  right  sides  of  these  equations  and  equate  them  to 
each  other.  Taking  into  account  equalities  (2.7.1)  and  (2.7.3),  we 
obtair. 


CD 


5^'!  -  8*t  (0  +  U\v  (I)  +  n  (0,  If  (*>  +  8n  (0;  ‘1  - 

'  -Mf(‘W(0;  M. 

-  80,  (0  +  Mr  co  +  8E  (0.  i*  (0  +  (0;  *1  - 
—MHO.  n*(0;  0- 


In  the  right  sides  of  the  last  equalities  we  can  use  the 
expansion  in  the  Taylor  series  of  the  function  according  to  increases 
in  two  of  its  arguments.  Retaining  only  terms  of  expansion  of  the 
first  order  relative  to  variations  &£(t)  and  Snft;  and  carrying  out 
simplifications,  we  will  arrive  at  the  following  totality  of  two 
linear  differential  equations: 


■  “*®i(0+-jjr4K0  +  -^*®n(0. 

-  8a, (l)  +  ^  81(0+  %-*l(0* 


The  desired  functions  of  these  equations  are  isochronal  varia¬ 
tions  of  coordinates  of  the  rocket  and  6n (t).  Initial  condi¬ 

tions  for  them  can  be  obtained  from  the  fact  that  at  the  instant  of 
the  launch,  i.e.,  when  t  «  0,  coordinates  and  velocities  of  the 
missile  are  the  same  both  in  real  and  in  calculated  motions.  Con¬ 
sequently,  according  to  equations  (2.7*1)  and  (2.7*2),  we  have 

8*(0)-0,  8«i(0)  =  0  (3.1.4) 

and  further 

«*<e»-^®-0.  (3.1.5) 


The  right  sides  of  equations  (3*1*3)  contain  variables  - 
isochronal  variations  of  projections  of  the  apparent  acceleration 
of  the  rocket  and  In  accordance  with  equations 

(2.7*3),  in  principle  they  can  be  obtained  aboard  the  missile  by 
means  of  a  continuous  formation  of  the  difference  between  the  real 
reading  of  the  appropriate  newtonmeter  and  its  current  calculated 
value . 


'*7 


Coefficients  of  equations  (3.1.3)1 


■rt  *jc 


•u 

■*r 


(3-1.6) 


should  be  considered  as  time  functions,  known  earlier  for  each 
specific  case  of  flight  of  the  missile.  In  order  to  obtain  these 
coefficients,  in  accordance  with  rules  of  expansion  in  Taylor  series, 
in  partial  derivatives  of  functions  /*  (|,  tj;  t)  and  /«(S.n:0  according  to 
variables  5  and  n,  it  is  necessary  to  produce  replacement  of  the 
latter  by  functions  £*(t)  and  Functions  5*(t)  and  n*(t)  are 

earlier  known  according  to  calculated  motions  of  the  missile  in  the 
powered-flight  section.  Thus,  if  aboard  the  missile  continuous 
(without  substantial  lag)  integration  of  equations  (3-1.3)  is  carried 
out  and  thus  magnitudes  of  functions  6u^(t)  and 

become  known,  then  the  left  side  of  the  basic  ballistic  equation 
(2.7.6)  can  be  constructed  by  means  of  multiplying  and  adding  devices. 
Similarly,  with  the  use  of  a  clock  mechanism  the  right  side  of  this 
equation  is  constructed. 


§  2 .  Approximate  Expressions  for  Isochronal 
Variations  of  Coordinates  and  Projections 
of  Velocity  of  the  Missile. 
Simplification  of  the 
Ballistic  Equation 


The  simplest  method  of  the  use  of  integrators  of  accelerations 
in  the  system  of  control  of  the  flight  range  of  ballistic  missiles 
is  based  on  the  simplification  of  differential  equations  (3.1.3). 

In  them  terms  containing  as  factors  the  isochronal  variations  them¬ 
selves  of  coordinates  6£(t)  and  are  dropped.  This  is  equivalent 

to  the  assumption  about  the  fact  that  the  effect  on  the  missile  of 
forces  of  gravity  with  its  actual  and  calculated  motions  can  be 
considered  practically  equal.  Such  an  assumption  is  admissible  only 


1  The  equality  in  (3.1.6)  follows  from  the  fact  that  projections 
cf  acceleration  or  the  force  of  gravity  are  partial  derivatives  of 
the  potential  of  gravity  according  to  appropriate  coordinates. 


with  sufficiently  small  deviations  in  the  real  motion  of  the  missile 
from  the  calculated,  i.e.,  at  increased  requirements  for  the  accuracy 
of  motion  cf  the  missile  according  to  the  assigned  program.  Equations 
(3.1.3)  are  replaced  in  this  case  by  the  approximate  equalities 

(3.2.1) 

whence  allowing  for  relations  (3.1.1)*  (1.5-^)  and  (3.1.5)*  there  are 
equations  for  variations  of  projections  of  the  velocity  of  the  missile 


*««>-  -*$!>. 
«■,(*)  -  *2jLSi>  - 


I 

(««,(»)*  =  6V,(0. 


(3-2.2) 


Quantities  &V^(t)  and  61’^ft),  which  enter  into  these  equations, 
can  Le  caileu  isochronal  variations  of  the  apparent  velocity  cf  the 
missile.  According  to  equations  (2.7-3)*  they  are  differences 

t "  (0  -  Vt,  (<)  ~  Vl  (I),  (0  -  v,  (()  -  v;  (0 ,  (3.2.3) 


in  which 

<  i 

*»<*>- r.(0-ja.(0*  (3.2.4 

are  projections  of  the  apparent  velocity  of  the  missile  during  its 
actual  motion  and 


i  i 

*  j *1(0^.  ^«(0  c  (3.2.5 

are  calculated  values  of  the  same  magnitudes. 

Integrating  in  turn  expressions  (3 .2.2)  with  respect  to  time, 
we  obtain  the  following  equations  for  isochronal  variations  of 
coordinates  of  the  missile: 


49 


Here 
path , 
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VI  I 

«<0  -  =  j  APa(l)A  =  W,(/), 

1 1  I 

*n(0  -  J J  te.(QA*  -  J  -  w.(0. 


(3.2.6) 


Here  6S^(t)  and  fis^ft)  are  variations  of  projections  of  the  apparent 
path,  i.e.,  differences 


where 


•*  (0  -  *10 — <  <0 .  *s,  (0  -  S.  (0  -  5*.  (0 . 

«  «■ 

*(<)«  jvt(o*.  j.(o-Jv,(o<» 


(3-2.7) 

(3-2.8) 


are  projections  of  the  apparent  path  of  the  missile  and 


*(o  -  j  v;(0J<.  <$*(0  -  j  v;(i)<fc 


(3.2.9) 


their  calculated  values. 


Let  us  substitute  now  expressions  (3*2.6)  and  (3-2.2),  respec¬ 
tively,  for  fiCft),  fin (t)a  6u^(t)  and  into  the  basic  ballistic 

equation  (2.7.6),  which  determines  the  instant  of  the  switching  off 
of  the  engine.  We  have 


^(0j  +  W*(0£  +  «ri(O|^  +  «f.(0j5j- 

— <*-o[«o£+«;(o£J. 


(3-2.10) 


Replacing  here  variations  of  projections  of  the  apparen*-  velocity 
and  apparent  path  by  their  expressions,  according  to  equations  (3.2.3) 
and  (3-2.7),  we  obtain 


l*  (0  -  Si  (0)  +  l st  (t)  -  S\  (01  + 

+ 1^(0 - ^(01  ^  + 1^,(0  -  v;(/; 
- (<-a‘)[-i(c.-)^  +  a;(0^-J. 


(3.2.1D 


i 


The  last  equation  can  be  considerably  simplified,  if  we  note 
that  correct  to  smallness  of  the  first  order  inclusively  the  follow¬ 
ing  expansions  in  Taylor  series  take  place: 

(3.2.1?) 

V\(t)  =V\(e)  +  (t-a  )**(')■ 

By  means  of  these  equalities  equation  (3*2.11)  is  reduced  to  the  form 


wo —^<01  ■£.+ -  j;  wi  £ + + 

+  +  (3*2.13) 


where  terms  dependent  on  time  in  evident  form  are  already  absent. 

Thus,  the  basic  ballistic  equation  (2.7.6)  can  be  substituted 
bv  the  approximate  equation  (3-2.10)  or  the  equation  (3.2.13) 
equivalent  to  it.  For  the  construction  of  the  left  side  of  equation 
(3.2.13),  the  presence  is  necessary  aboard  the  missile  of  two 
integrators  of  accelerations,  axes  of  which  during  the  whole  powered- 
flight  section  of  flight  should  retain  directions,  parallel  respec¬ 
tively  to  axes  F„  and  n  of  the  nonrotating  system  of  coordinates  CnC. 
Furthermore,  there  must  be  a  computer,  which  Includes  in  its  com¬ 
position  two  additional  elements  for  the  time  integration  of  current 
readings  of  the  very  integrators  of  accelerations  and  a  special 
element  for  the  reproduction  of  the  calculated  current  values  of 
projections  of  the  apparent  path  S£(t)  and  S^(t). 

§  3.  Construction  of  the  Ballistic  Equation  by  Means 
of  the  Use  of  Readings  of  Two  Integrators 
of  Accelerations  with  Special 
Orientation  of  Their  Axes 
of  Sensitivity 


There  can  be  considerable  interest  in  the  possibility  of  the 
reduction  in  the  number  of  elements  of  the  computer  system  of 
inertial  determination  of  the  flight  range  of  the  ballistic  missile 


by  means  of  the  proper  selection  of  orientation  of  axes  of  sensitivity 
of  integrators  of  accelerations.  It  turns  out  that  there  can  be 
produced  a  total  of  one  additional  integration  of  current  readings 
of  integrators  of  accelerations,  and  not  two,  as  in  the  method  of 
the  construction  of  the  ballistic  equation  given  in  the  previous 
paragraph.  Per  this  purpose,  as  will  be  shown  below,  axes  of 
sensitivity  of  integrators  of  accelerations  should  be  parallel  to 
certain  directions  fixed  in  each  specific  case  of  flight  of  the 
missile  (so-called  A-  and  u-directions ) . 

Another  method  of  simplification  of  the  system  of  inertial 
control  of  the  range  with  the  use  no  longer  of  an  integrator  of 
accelerations,  and  the  meter  of  the  apparent  accelerations  with  the 
changing  orientation  of  the  axis  of  sensitivity,  is  stated  in  the 
following  paragraph.  Subsequent  modification  of  this  method  is 
given  in  §  5  of  this  chapter. 

Tnus ,  let  us  position  the  axis  of  sensitivity  of  one  of  the 
Integrators  of  accelerations  in  the  plane  of  the  nonrotating 
system  of  coordinates  6ns  at  a  certain  constant  angle  A  to  the  axis 
£.  The  direction  Itself  of  the  axis  of  sensitivity  of  this  Integra¬ 
tor  will  be  called  the  A-direction  (Pig.  10). 

On  the  basis  of  equation  (1.6.5),  after  the  replacement  in  it 
of  letters  x,  y  and  v,  respectively,  by  £,  n  and  X  and  also  angle  a 
by  angle  X ,  we  obtain  equation 

— +  (3.3-1) 

whi ih  expresses  the  projection  of  the  apparent  velocity  of  the 
missile  on  the  mentioned  A-direction  by  its  projections  on  axes  6 
arni  ri.  It  is  obvious  that  in  accordance  with  formula  (1.4.7) 

i 

(3-3.2) 

where  a^(t)  —  component  of  apparent  acceleration  along  the  same  A- 
' i i reef  Ion. 
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Fig-  10. 


Thus,  quantity  is  the  current  reading  of  the  integrator  of 

accelerations  with  the  axis  of  sensitivity,  oriented  in  direction  X. 

Let  us  select  angle  X  so  that  there  would  take  place  equalities 

—  tfcwX,  j— ™Af*tafc.  (.  3  -  3  •  3 ) 

Here,  as  it  is  easy  to  see  (see  Fig.  10),  quantity 

is  for  each  specific  case  of  flight  of  the  missile,  an  earlier 
known  quantity. 

Using  equations  (3-3.1)  and  (3-3.3),  let  us  transform  the  sum 
of  the  last  two  corresponding  terms  of  the  left  side  of  equation 
(3.2.13).  We  have 

F*(I).£L  +  V,(l)  jiL  -  M  lKi(<)eojX  +  K,(0»ta  M  - 

(3.3.5) 

Tivus  ,  the  reduced  sum  correct  to  the  constant  factor  M  is 
determined  by  the  reading  of  the  integrator  of  accelerations,  the 
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axis  of  sensitivity  of  which  is  located  along  the  X-airection. 


Ac  a  result  of  relation  (3*3-5)*  we  have 

<«’>  R-j  +  V';  (3)  £  -  MVl  (a*).  (3.3.6) 


Similar  to  relation  (3- 3- 5),  we  oan  obtain  the  following 
equal! ties 

(0  -  J  (0  £  -  N  | St  (t)  co.  H  +  J,  (<)  sin  (!)  -  NS,  (t) 

<(*)$  +  •?«(')  ^  -  .V  t^(<)co.«  +  S'jt),\ntf-Ns;(i).  (3.3.7) 

In  them  angle  u  (Pig.  11)  is  determined  by  means  of  equations 

-yj-  .Voosh.  —  A/  »ln  p,  (3.3.8) 

and  the  constant  factor  N  -  by  equality 

*-/(£)'  +  (£)*  (3.3.9) 


In  right  sides  of  equations  (3-3.7)  S ft)  -  projection  of  the 

vector  of  the  apparent  path  of  the  missile  on  the  so-called  u- 

dlrection,  which  forms  the  constant  angle  p  with  axis  £  (Fig.  11) » 

and  S*(t)  -  current  calculated  value  of  this  projection.  By  analogy 
u 

with  equations  (1.6.7)  and  (1.4.7)* 

1  11 

ii.j.  1 0) 


Let  us  turn  now  again  to  equation  (3.2.13).  Taking  into  account 
equalities  (3.3.6),  (3.3.6)  and  (3-3. 7),  we  obtain 


mvi  (i) + n  1  (<)  -  s;  (oi  -  m  n  (8*). 


54 
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(3.3.11) 


The  right  side  of  the  last  equation  is  a  constant  determined 
for  each  specific  case  of  flight  of  the  missile.  This  magnitude 
correct  to  the  constant  factor  M  is  equal  to  the  calculated  value 
of  projection  on  the  X-direction  of  the  vector  of  apparent  velocity 
of  the  missile  at  the  calculated  Instant  of  the  switching  off  of  its 
engine  t  =  a*.  If  in  equation  (3* 3* 11)  we  introduce  another  designa¬ 
tion  1 


then  it  can  be  reduced  ,to  the  following  final  form: 

r»(0+#[<M0-J&(01 -»*(*•>.  (3'  3.13) 

For  the  construction  aboard  the  missile  of  the  left  side  of  the 
ballistic  equation  (3.3*13),  an  additional  integration  of  current 
readings  of  only  one  Integrator  of  accelerations  is  necessary  with 
the  axis  of  sensitivity  oriented  in  the  u-direction.  Furthermore, 
the  reproduction  of  the  calculated  values  5^(t)  of  this  repeated 
Integration  is  necessary. 


1  Let  us  note  that  coefficient  p,  similar  to  factor  1/t  in  equa¬ 
tion  (1.6.11),  has  a  dimensionality  opposite  to  time  (in  particular 
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§  n. 


Construction  of  the  Ballistic  Equation  wltn 
the  Help  of  a  Special  Meter  of 
Apparent  Accelerations  of 
the  Computer 

Let  us  examine  another  conversion  of  the  ballistic  equation 
(3.2.13),  as  a  result  of  which  the  possibility  of  solution  to  the 
problem  of  inertial  control  of  flight  range  of  the  ballistic  missile 
by  means  af  a  single  newtonmeter  is  explained.  Its  readings  in  the 
computer  should  be  integrated  after  preliminary  multiplication  by 
^  certain  assigned  time  function.  The  axis  of  sensitivity  of  the 
newtonmeter  in  turn  change  its  orientation  of  the  nonrotating  system 
of  coordinates  according  to  the  assigned  law. 

Let  us  turn,  first  of  all,  to  the  conversion  of  the  first  two 
corresponding  terms  of  the  ballistic  equations  (3.2.13).  According 
to  equations  (3.2.1*)  and  (3.2.8)  and  also  the  known  Cauchy  equation 
of  the  conversion  of  the  repeated  integral  into  a  simple  one  contain¬ 
ing/  the  ur p< t'  threshold  in  the  form  of  the  parameter  in  the  sub- 
integral  expression,  we  have 


.  i .  n  i 

St  <0  -  $  ?'  <0  *  -  ( \  (0  **  -  \  (1  -  r)  at  (r)  dr, 

•  ••  • 


(3.^.1) 


where  i  -  now  variable  of  integration.  Representing  here  the 
difference  t  -  t  In  the  form 

I_T  =  «_o,*)  -T  (0‘-  T).  (3.4.2) 


we  obtain  equality 


9  9 

st  (0  =  (<  -  0-)  (T)  dx  +  ^  (3-  - 1)  a.  ( t)  dr.  (3.4.3) 

•  • 

Returning  in  the  last  integral  of  the  right  side  of  this 
equality  to  the  initiai  variable  of  integration  t  ar.  i  taking  into 
account  equation  (3.2.**),  we  obtain 


( 3 .  k .  <0 


«%<»)  -  <f  -  a1)  K,(i)  +  <o*  -  0  *i  {<) *■ 


Let  us  note  further  that  correct  to  smallness  cf  the  first  order 
inclusively  we  can  assume  that 


•S»  <0  «  At  («)  +  <<- •*>  W) 


(3.*. 5) 


Equating  to  each  other  differences  of  the  left  and  right  sides 
of  the  last  two  equalities,  we  have 


I 

*(0— $(<)  ~  j(a*-i)«t(0<ft-Ai(9‘)  + 

+(<-«*)in(l)-V;<a*)l. 


v  3- 4.6> 


The  last'terrr.  of  the  right  side  of  equality  C  3 .  h .  6 )  has  a 
second  order  of  smallness  and  can  be  omitted.  Actually,  as  a 
res  i 5  \  :.f  t.*e  proximity  cf  the  real  anu  calculated  instants  r  *  c 

and  t  «  o*  of  the  switching  off  of  the  engine  quantity  V^ft)  differs 
Uttle  from  V^(o*),  and  this  latter  is  distinguished  from  f£(o>*)  by 
a  small  Isochronal  variation  6K|(a#).  Therefore  the  difference 
V^(t)  -  y«(o»)  is  a  3mal3  magnitude,  which  in  equation  (3- 4.6)  is 
multiplied  ir.  turn  by  another  3mall  difference  t  -  o*. 

Thus,  correct  to  smallness  of  the  second  order 


*(‘>^$(0“  f («■-<)< H(0*— Sl(o’) 


(3-^.7; 


In  perfect  analogy 


<)«•(»)* -  -£<»*) 

< 


(3.^.o) 


Substituting  differences  (3-^.7)  and  (3-^.8)  into  equation 
( 3.  ?.]  3)  and  replacing  in  it,  furthermore ,  quantities  V^(t)  and  vr^t' 
by  their  ix-precentations ,  according  to  equations  (3-2.10,  we  obtain 
aft>  r  i;.=r  11  fi  cations 


(3- **.9) 


[\[£i  +  -  o  ilt(0  +  [£  4-  «.•  - 1)  £1 «.  («))  *  * 


T > i «.  J'igh'.  aide  of  the  last  equation  for  the  specific  flight  of 
the  rocket  is  a  certain  constant,  which  we  designate  by  the  letter 
b.  in  accordance  with  equations  similar  to  (3*3*6)  and  (3*3*7),  this 
constant  can  also  be  represented  in  the  form 

C"  )-r'2j-<S|(a)  +  'j“‘y*(o),“ 

-  MVU*)  +  NSU*)-  (3* 4.  10) 


Conversion  of  the  integrand  expression  of  the  left  side  of 
equation  (3*4.9)  can  be  produced  by  the  same  methods  which  were 
incorporated  in  the  previous  section.  Let  us  present  in  this  equa¬ 
tion  alternating  factors  before  projections  of  the  apparent  accelera¬ 
nt:..  u  :  ar.d  j  (t )  in  '.he  following  manner; 

*  , 


jjL  +  (3*  —  1T(«)oo»k(<). 

jjL  -f  (a*  —  l)^j|  ■  K (<) »ln  x  (I). 


(3*4.11) 


[n  accordance  with  the  last  equalities. 


*«*«)- jr - 5r- 

£;+<-*-  *>n 


(3*4.12) 


The  second  equation  (3*4.12)  determines  <(t)  changing  with  time 
angte  between  the  so-called  x-direction  and  axis  £  (Fig.  1?). 


Pig.  12. 

Using  equations  (3.3-3)?  (3*3.8)  and  (3-3.12),  expression 
(3.4.12)  for  variable  K(t) ,  after  comparatively  simple  conversions, 
can  also  be  presented  in  the  form 


K(ty=M  y  i  +  ip  (3'  —  t)  cos  (fi  - 1)  +  (3-  -  <)*. 


(3.4.13) 


Let  us  produce  in  the  left  side  of  equation  (3. 4.9)  replacement 
of  coefficients  according  to  equalities  (3.4.10)  and  ,3-4.11). 
obtain  the  relation 


« 

\  K  (0  l«l  (0  COS K (/)  +  a,  (/)  sin  x  (/)|  dt  -=  C, 


(3.4.14) 


where  in  integrand  expression  the  sum 


8'(f)eosx(t)  -4-  a„(/)  sinx{/)  =  o,(f) 


(3.4.15) 


is  the  projection  of  the  apparent  acceleration  of  the  missile  towaru 
the  mentioned  variable  ic-direction  (Pig.  12),  and  constant  C  is 
expressed  by  equation  (3.4.10).  As  a  result  we  obtain  the  ballistic 
equation  iri  the  following  form: 


-  C. 


(3.4.16; 


•Ions tract  Lon  of  the  left  side  of  equati  on  ( 3. 4.  *6)  aboard  the 
missile  requires  the  oresenee  of  highly  accurate  fel icw-up  systems 
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located  on  a  gyroscopic  stabilizer  for  the  change  in  orientation  of 
the  axis  of  sensitivity  of  the  meter  of  acceleration  according  to 
the  assigned  program. 

Of  course,  approximation  methods  of  the  formation  of  the  left 
side  of  the  ballistic  equation  (3.4-16)  by  means  of  the  usual  integra¬ 
tors  of  accelerations,  for  example,  by  means  of  replacement  of  func¬ 
tions  K(t)  and  <(t)  by  their  certain  mean  values. 

§  5.  Inertial  Control  of  Range  of  the  Ballistic 
Missile  by  Means  of  Longitudinal  and 
Standard  Integrators 
of  Accelerations 

In  systems  of  control  of  ballistic  missiles  one  can  use  the 
so-called  standard  integrators  of  apparent  accelerations,  the  axes 
of  sensitivity  of  which  are  guided  in  plane  £n  at  the  assigned  angle 
to  the  axis  K  according  to  the  program  of  flight.  The  control  of 
fright  controls  is  produces  in  this  case  so  that  together  with  the 
fulfilling  of  the  assigned  program  of  the  charge  In  pitch  of  the 
missile  (i.e.,  change  In  angle  between  its  longitudinal  axis  and 
axis  £)  the  reading  on  the  standard  Integrator  would  be  reduced  to 
zero.  The  angle  between  the  axis  of  sensitivity  of  the  standard 
integrator  and  axis  of  the  missile  itself  is  selected  close  to  a 
straight  line. 

I'he  use  of  the  standard  integrator  of  the  apparent  accelerations 
in  the  system  of  control  of  motion  of  the  missile  leads  to  stabiliza¬ 
tion  in  the  assigned  direction  of  the  resulting  force  of  thrust  of 
the  engine  and  aerodynamic  forces  acting  on  the  rocket.  Thus,  the 
direction  of  the  vector  of  the  apparent  acceleration  of  the  missile 
( h-  -vie  ver ,  of  course,  not  its  magnitude)  Is  stabilized.  With  the 
known  approximation,  being  distracted  from  errors  of  the  system  of 
control,  it  can  be  considered  that  the  standard  Integrator  controls 
the  flight  of  the  missile  so  that  the  projection  of  the  apparent 
acceleration  on  the  axis  of  sensitivity  of  this  integrator  would  be 
equal  to  zero. 
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Let  us  present  In  the  ballistic  equation  (3.^-l6)  of  the 
previous  section  quantity  o^(t)  in  the  form  of  the  product 

*.(<)  —  «(0co*t(0  (3-5.1) 

Here  y(t)  is  the  angle  which  forms  the  vector  of  the  apparent 
acceleration  alt)  with  K-direction  inclined  at  angle  <(t)  to  axis 
5  (Fig.  13)- 


Function  <(t)  is  determined  by  the  second  equality  (3.^.1?), 
arid  the  direction  of  the  apparent  acceleration  is  changed  according 
to  the  assigned  law.  Because  of  this  the  magnitude  of  angle  y(t) 
should  be  considered  as  a  known  function  of  time. 

Producing  in  equation  (3.**.16)  the  replacement  of  quantity 
a^(t),  according  to  equation  (3-5-1)  we  will  arrive  at  the  follow¬ 
ing  modification  of  the  ballistic  equation: 

\  -C.  (3.5.?) 

in  which  the  alternating  coefficient 

^(0  —  (3.5.3) 
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should  be  considered  the  known  time  function  for  the  given  cal¬ 
culated  case  of  flight  of  the  missile. 

For  construction  of  current  values  of  the  left  side  of  equation 
(3.5.2)  aboard  the  missile  just  as  earlier  (see  S  4  of  this  chapter), 
integration  of  current  values  of  the  apparent  acceleration,  multiplied 
by  the  assigned  time  function  is  required.  However,  now  there  is  no 
need  in  special  high-precision  equipment  for  the  change  in  orientation 
of  the  axis  of  sensitivity  of  the  meter  of  the  apparent  acceleration 
of  the  missile.  Actually,  if  this  newtonmeter  is  rigidly  connected 
to  the  standard  integrator  so  that  axes  of  their  sensitivity  would 
be  perpendicular  to  each  other,  then  with  the  proper  accuracy  of  the 
control  of  flight  of  the  missile  the  apparent  acceleration  of  the 
latter  (more  accurately,  apparent  acceleration  of  location  of  these 
instruments )  will  be  wholly  directed  along  the  axis  of  sensitivity 
of  the  mentioned  newtonmeter. 

yiU.  some  less  of  accuracy  of  flight  of  the-  missile  the  newton- 
meter  can  be  fastened  directly  aboard  the  missile,  the  guiding  axis 
of  its  sensitivity  in  parallel  to  the  longitudinal  axis  of  the 
mis.-.ile.  The  latter  is  deflected  from  the  direction  of  the  apparent 
acceleration  of  the  missile,  i.e.,  from  the  direction  of  the  resulting 
force  of  thrust  of  the  engine  and  forces  of  aerodynamic  actions,  as 
a  rule,  at  small  angles.  Because  of  this,  approximately  to  lay 

'  «,(j)sk«(0.  (3*5.4) 

where  a„(t.'  -  projection  of  apparent  acceleration  of  the  missile  on 

•J 

its  longitudinal  axis. 

If,  in  accordance  to  the  last  approximate  equality,  we  produce 
the  appropriate  replacement  in  equation  (3*5.2)  and,  furthermore, 
substitute  function  Q(t)  by  its  certain  mean  value  (j ,  then  we  will 
obtain  the  approximate  equation 


$5  «.(<)*-??•<  0 


C 


(3*5-5) 
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for  the  determination  of  the  instant  of  switching  off  of  the  engine 
of  the  missile.  Here  expression 

I 

MO*  (3.5.6) 

• 

is  the  current  reading  of  the  integrator  of  acceleration  with  the 
axis  of  sensitivity  parallel  to  the  longitudinal  axis  of  the  missile.1 
Such  an  Integrator  is  called  longitudinal.  As  is  known,  it  was  used 
in  the  guidance  system  of  the  German  missile  V-2. 

Let  us  return  again  to  equation  (3*5.2)  and  examine  additionally 
certain  approximate  methods  of  construction  of  its  left  side  on  the 
basis  of  current  readings 


t 

(3-5.7) 

of  the  integrator  of  accelerations  with  the  axl.  of  sensitivity 
located  perpendicular  to  the  axis  of  the  standard  integrator.  Con¬ 
sidering  the  latter  equality  the  left  side  of  equation  (3.5*2';  is 
converted  to  the  form 


) 


Q(t)a(t)it-$Q(t)^l 


it 


(3.5.8) 


and  can  be  Integrated  by  parts.  As  a  result,  taking  into  considera¬ 
tion  that  k(0)  »  0,  we  obtain,  according  to  relations  (3*5.8)  and 
(3.5.2),  the  ballistic  equation  in  the  following  form; 

C(  0?(0-  5  Q'(t)  Fvw  c.  ( 3 . 5 . 9 ) 


‘Let  us  note  that  quantity  V^Ct)  in  equation  (3.5.6}  because 

of  the  variability  of  direction  of  the  longitudinal  axis  is  not  the 
^rejection  on  this  axis  of  the  apparent  velocity  of  the  missile  in 
the  nonrotating  system  of  coordinates  (for  more  detai)  on  this  see 
5  4  c  r  Chapter  I ) . 
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The  construction  aboard  the  missile  of  the  left  side  of  the  last 
equation  can  already  be  produced  by  means  of  only  one  Integrator  of 
acceleration. 

Let  us  note  in  conclusion  that  the  same  method  of  integration 
by  parts  can  be  used  for  constructing  the  left  side  of  the  ballistic 
equation  (3.4.16). 
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CHAPTER  IV 


INERTIAL  CONTROL  OP  THE  RANGE  OF  PLIGHT  OF  THE 
BALLISTIC  MISSILE  ALLOWING  POR  THE 
CHANGE  IN  THE  FORCE  OF  GRAVITY 


§  1 .  One  Method  of  the  Solution  of  Differential 
Equations  by  Which  Isochronal  Variations  of 
Coordinates  of  the  Missile  with  It~ 
Motion  In  the  Powered-Flight 
Section  Are  Determined 


In  the  previous  chapter  the  basic  ballistic  c-^ailon  (2. 7.6;  was 
transformed  to  such  a  form  at  which  its  left  side  could  be  constructed 
aboard  the  missile  with  the  help  of  integrators  of  accelerations.  In 
this  case  Isochronal  variations  of  coordinates  and  projections  of 
velocity  of  the  missile,  which  enter  into  the  composition  of  equa¬ 
tion  (2.7.6),  were  expressed  by  the  approximate  equations  (3.2.6) 
and  (3.2.2).  This  corresponded  to  the  neglect  in  differential  equa¬ 
tions  (3*1*3)  of  terms  containing  as  factors  the  isochronal  varia¬ 
tions  5^(t)  and  oni't;  themselves.  As  was  already  indicated  in  Chapter 
III,  the  mentioned  terms  of  equations  (3*1.3)  take  into  account  dur¬ 
ing  the  determination  of  current  coordinates  c-f  the  rocket  changes  in 
acceleration  of  the  force  of  gravity  because  of  the  noncoincidence  of 
its  actual  motion  with  the  calculated.  Error  appearing  from  such  a 
simplification  of  equations  (3- 1.3)  in  the  determining  of  the  range 
jf  fiight  of  the  missile  is  small  only  at  rather  small  deviations  in 
its  real  motion  from  the  calculated.  Otherwise  it  is  necessary  to 
solve  approximately  aboard  the  missile  differential  equations  (3.1.3) 
allowing  for  terms  reflecting  the  effect  of  the  change  in  accelera¬ 
tion  of  tr.r  force  of  gravity. 
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Below  for  the  solution  of  equations  (3.1.3)  &  similar  approxima¬ 
tion  method  is  used,  and  it  was  shown  to  the  author  by  corresponding 
member  of  the  Academy  of  Sciences  of  the  Ukrainian  SSR  Yu.  D.  Sokolov. ' 
By  means  of  this  method  the  problem  of  inertial  control  of  range  car; 
be  quite  accurately  solved  by  the  means  given  in  Chapter  III  with 
the  help  of  the  application  of  standard  integrators  of  accelerations 
and  the  simplest  computers. 

The  ides  of  the  method  consists  in  the  conversion  of  equations 
*3.1.3;  to  their  equivalent  integral-differential  and  integral  rela¬ 
tions  with  suosequent  replacement  under  the  sign  of  integral  of 
expressions  of  desired  functions  by  their  simplest  approximations, 
which  satisfy,  however,  initial  and  final  conditions  of  the  problem. 

In  accordarice  with  this,  let  us  integrate  the  right  and  left  sides 
of  equations  (3.1*3)  with  respect  to  time  from  the  Initial  instant 
t  -  0  up  to  the  unknown  until  the  instant  of  switching  off  of  the 
engine  t  *  o,  which  provides  the  hitting  of  the  missile  on  the  target. 

1  -king  Into  ac-  i:.  this  case  initial  conditions  (3.1.51  for 

isochronal  variations  of  projections  of  the  velocity  of  the  missile, 
we  obtain  two  integral-differential  relations 

-  fa* (3)  -  (j)  +  ^  «(<>*  +  \  Q  *1  (0 

J~ftn(0<&.  (A.  1.1) 

in  which,  similar  to  equalities  (3.2.2), 

j &»,(/)  A  (A.  1.2) 

the  isconronal  variations  of  projections  of  the  apparent  velocity 
of  the  missile  on  axes  C,  and  n  at  Instant  t  »  o. 


'See,  for  example,  Yu.  D.  Sokolov.  One  method  of  observed 
solved  linear  integral  differential  equations.  —  Doo.  AN  URSF;, 
1955,  No.  2. 
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If,  however,  the  right  and  left  sides  of  equatJons  (.3.1.3)  are 
integrated  with  respect  to  time  twice:  one  time  within  the  limits 
of  t  =  0  to  the  current  instant  t  and  the  second  -  also  from  t  *  0 
to  the  instant  of  switching  off  of  the  engine  t  -  o,  then,  taking 
account  initial  conditions  (3.1.4)  and  (3.1.5),  we  obtain  an  addi¬ 
tional  two  integral  relation: 


*•!<«)-  «,(«)+ 


(4.1.3) 


where  in  turn,  in  accordance  with  equalities  (3.?. 6), 


AS*  (a)  -  j  M j  (*)*•. 

M,(«)  -  j »%( 0*  -  [  j  *.(*) ** 


(4.. .;) 


isochronal  variations  of  projections  of  the  apparent  path  cf  the 
rocket,  which  refer  to  the  instant  t  *  o. 

In  accordance  with  the  mentioned  method,  one  should  substitute 

in  the  right  sides  of  relations  (4.1.1)  and  (4.1.3)  functions  65(t> 

and  6n (t)  by  (heir  approximate  representations,  which  turn  into  zero 

at  the  initial  Instant  of  time  and  respectively  into  3£(o)  and  6n(o' 

when  t  =  a.  Then  values  of  quantities  6£(o'',  5n(cr),  6ur<a)  and 

‘  > 

6^(0),  found  in  accordance  with  the  mentioned  relations,  prove  to 
be,  as  a  rule,  more  accurate  than  those  calculated  according  to 
equations  (3.2.2)  and  (3.2.6),  founded  upon  approximate  equalities 
(3.2.1). 

;  ?  is  possible  to  take  as  approximate  representations  6£(  •.)  and 
in  right  sides  of  relations  (4.1.1)  and  i4.].3),  for  example, 
th>  1  :■  .  inear  approximation 


n7 


(4.1.5) 


**(')«  A.l(/)-ftl|(3)|. 

where  o  -  instant  of  time  for  which  the  determination  of  desired 
functions  is  produced;  in  this  case  o  —  instant  of  termination  of 
the  powered-flight  section  of  flight  of  the  missile. 

The  best  results  should  be  expected  with  the  quadratic  approxima¬ 
tion  of  the  form 

A*  («)•£.  An(/)-A«i(a)-J-.  (4.1.6) 

The  reason  for  this  is  that  functions  (4.1.6),  unlike  functions 

(4.1.5) ,  satisfy  simultaneously  initial  conditions  (3*1.4)  and 

(3.1.5) ,  which  concern  both  variations  themselves  of  the  coordinates 
and  their  time  derivatives  at  the  instant  t  *  0.  Functions  (4.1.5) 
do  not  satisfy  initial  conditions  (3.1.5). 

Leading  to  an.  even  greater  accuracy  should  be  the  assignment 
of  functions  and  6n(t)  in  the  form  of  the  following  polynomials 

of  the  third  power 

«*«  -«(»)(*  4  ~  *  4)  +  Amt  (*)•«  (4  -  4) . 

4^(1)  *■  4tj(<j)  (3  4—  2  4)  +  Am,(o)-«  (4  ~  4)’  (4.1.7) 


These  polynomials  not  only  satisfy  initial  conditions  (3.1.4)  and 

(3.1.5),  but,  furthermore,  when  t  *  0  turn,  respectively,  into  55(o) 
and  5n(o)  and  their  derivatives  —  into  6u^(a)  and  61^(0). 

&  2 .  Approximate  Solution  of  Differential  Equations 
for  Variations  of  Coordinates  Using  the  Quadratic 
Approximation  and  Also  the  Approximation  in  the 
Form  of  Polynomials  of  the  Third  Power 

Let  us  examine  from  the  beginning  the  solution  of  differential 
equations  >'3.1.3)  which  corresponds  to  the  quadratic  approximation 
V‘i.1.6).  ."uostituting  into  right  sides  of  relations  (4.1.3) 


? 


expressions  for  6£(t)  and  5n (t) ,  according  to  equations  (4.1.6), 
we  obtain 


«(•)  -  Mk(«)  +*Jr*(S  + 


Let  us  introduce  here  designations 


(4.2.1) 


(4.2.2) 


where,  because  of  relations  (3.1.6), 

%*<•>-**<«»•  (4.2.3) 

As  a  result  let  us  arrive  at  the  two  algebraic  equations: 


ae<o) -  w»(«) +AttW«6(o)+*i,(9)4n(»).  , ^  2 

*1  («)  -  «•(*)  +  **  (ff)  H  («)  +*«  (a)  4«l  («) 

with  respect  to  the  desired  quantities  $£(o)  and  6n(o). 

Values  of  coefficients  4^(0),  h^(a)  »  4^(0)  and  4^(0)  with 
an  accuracy  sufficient  for  practice  can  be  taken  with  argument  o  = 

-  a*.  Actually,  the  real  duration  of  the  powered-flight  section  o 
only  by  a  small  magnitude  is  distinguished  from  the  calculated  o*. 
Isoonronal  variations  6£(o)  and  6n(o)  should  also  be  considered  as 
small  magnitudes.  Consequently,  the  error  which  occurs  from  the 
mentioned  replacement  in  coefficients  4^(0) ,  h^(a)  and  4^(0), 
which  stand  in  the  right  side  of  equations  (4.2.4),  of  the  value  of 
argument  o  by  o*,  has  the-  second  order  of  smallness. 
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Let  us  designate  coefficients  (*4.2.2)  when  o  ■  a*  by  the  same 
letters  but  without  indication  of  the  argument.  With  the  made 
simplifications  and  replacement  of  the  designation  of  argument  o  by 
t  equations  (*4.2. *4)  take  the  following  form: 


«l(t)  -«*(!)  +  Ml  (!) + (0. 
*  (0  -  w,  (<)  +  (I)  +  Mn  (I), 


(4.2.5) 


where  *. .  ^ ,  h^r.  ■  are  quantities  known  for  the  given  cal¬ 

culated  flight  of  the  missile,  which  are  small  in  comparison  with 
onity  (see  5  3  of  this  chapter) .  Because  of  this  the  totality  of 
equations  (*4.2. 5)  is  most  convenient  to  solve  by  means  of  iterations 
taking  6C(t)  and  6n(t)  in  the  right  side  of  equations  (*4.2.5)  the 
basic  expressions  represented  by  equations  (3.2.6),  namely 


&t(t)-Ui(0.  *n(0-W,(0. 


(*4.2.6) 


by  substituting  these  expressions  into  right  sides  of  equations 
(*4.2.5),  we  arrive  at  the  following  equations  for  current  isochronal 
variations  of  coordinates  ■:(( t)  and  6n(t): 


H(i>  -  (1  +  *u) (0  +  M?.  (0. 
«n  (t)  -  M*  (0  +  (t  +  M  (»). 


(*4.2.7) 


it  is  easy  tc  be  convinced  thftt  subsequent  iterations  of  equa¬ 
tions  (4j,?.5)  lead  accurately  to  the  same  equations  for  and 

C'n(t),  if  only  in  the  calculations  we  drop  terms  with  squares  and 


products 


,  and  h  . 

nt  nn 


To  search  for  variations  of  projections  of  the  velocity  of  the 
missile  6u^(t)  and  let  us  substitute  into  right  sides  of 

relations  (*4.1.1)  expressions  for  6£(t)  and  tfnCt),  represented  by 
equations  of  the  same  quadratic  approximation  (*4.1.6).  We  obtain 
for  the  d'v,  i  reo  magnitudes  these  equations 


+  «(<») -srj 

*  at  *  at  (4.2.8) 

*M«)  =  «M®)  +  flSW-jr  j-J-Aft  +  *1 

Let  us  introduce  designations 

(4.2.9) 

*<*(a>  =--  *«W  =  "IrjjlRf’***- 

Then  equations  (4.2.8)  are  presented  in  the  form 


«■»<«)  -  »»(«)■+  ftt(«)*6(«) + ft,(«)»n(o). 

a«,(o) = «%(<»)+  r<(a)«l(a)  +  (o)  *n  f«). 


Similar  to  the  previous  in  equations  (4.2.10),  instead  of 
coefficients  g^(o),  g^Co)  *  g^to)  and  gnn(°)>  is  possible  to 
take  their  values  at  the  calculated  instant  of  the  termination  of  the 
powered- flight  section  of  flight  of  the  missile  a*.  We  designate 
these  values,  respectively,  by  g^,  =  g^  and  g^.  Let  us  sub¬ 

stitute  further  into  equations  (4.2.10)  expressions  (4.2.5)  for 
variations  of  coordinates  65(a)  and  6n(a)  and  disregard  the  products 
of  coefficients  g  and  h  with  arbitrary  indices.  Ultimately  we 
arrive  at  equations  for  the  desired  Isochronal  variations  of  velocity 
of  the  missile  6u^(o)  and  6u^(o)  with* respect  to  immobile  system  of 
coordinates  £n£.  Changing  in  them  the  designation  of  argument  o  by 
i,  we  obtain 


*«,  (|)  .  «F,  (I)  +  g*fiSt  <t)  +  (/).  .4.2.11) 

Equations  (4.^.7)  and  (4.2.11)  are  the  approximate  solution  of  the 
totality  of  differential  equations  (3.1*3)  in  a  f’or m  similar  for  the 
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application  in  systems  of  inertial  control  of  flight  range  of 
ballistic  missiles. 

Somewhat  more  bulky  are  equations  for  isochronal  variations  of 
coordinates  and  projections  of  velocity  of  the  missile,  if  we  take 
polynomials  (4.1.7)  for  the  initial  approximations  in  equations 
(4.1.1)  and  (4.1.3).  As  a  result,  similar  to  the  case  of  the 
quadratic  approximation,  it  is  possible  to  arrive  at  the  following 
equations  for  variations  of  coordinates: 


« (0  -  (1  +  3A*  -  2fa)  iSt  (I)  +  (3*,- 2/j,)  6S,  ( t )  + 

+  a*  (it*- h*)  W*  (0. 

*1(0  -  (»*-  2/a)  +  (1  +  3A„_  2/„)  6S,  (t)  + 

+  «•</*-*«)  MY<0  +  «*(/„-  (4.2.12) 

and  for  variations  of  projections  of  the  velocity,  correspondingly 

+  II  +  a*  ('a — gtt)\  4K(0  +  °*  (k*  —  ft*)  W*  (0. 

*».«  -  {3ra-2^)«»(0  +  (3y„~2/„)  85,(0  + 

+  <J*(‘\*-r*)«\(')  +  li  +  a*  «',(*).  (4.2.13) 

In  equations  (4.2.12;  and  (4.2.13),  besides  designations  already 
introduced  in  this  section,  quantities  L  in£,  *rr|  anci  «?’££;> 

Jnn  are  values  of  Integrals 

«« («)  *=  3  "s* 

•  #/  (4.2.14) 

it*  (<0  =  /<(«)  “  •gr^7^-‘,rf' 

fn(v)  ~  ~Sf i"  (*)  *  *3r  \ '  "*f 

II  •• 


(*4.2.15) 


which  they  take  with  argument  o  equal  tj  i.e., 
instant  of  switching  off  of  the  engine. 


at  the  calculated 


$  3.  Simplification  of  Equations  for  Isochronal 
Variations  of  doordlnates  and  Projections  of 
Velocity  of  the  Missile  at  Extensions  of 
the  Powered-Plight  Small  in  Comparison 
with  the  fladius  of*  the  Earth 

Equations  (4. 2.7 )  and  (4.2.10)  and  also  (if. 2. 12)  and  (4.2.13), 
obtained  in  the  previous  section  for  variations  of  coordinates  and 
projections  of  velocity  of  the-missile ,  can  without  great  harm  for 
their  accuracy  be  represented  in  simpler  form,  which  dees  not  require 
preliminary  calculation  of  integrals  of  the  type  (4.2.2),  (4.2.9) » 
(4.2.14)  and  (4.2.15).  Actually,  partial  derivatives  of  the  accelera 

tions  of  the  force  of  the  earth’s  gravity  xp  »  taking 

part  in  these  integrals,  are  functions  of  moving  coordinates  of  the 
calculated  motion  of  the  missile  and,  consequently,  finally  -  time. 
However,  with  a  small  extension  of  the  powered- flight  section  .  •:  the 
trajectory  of  the  missile,  in  comparison  with  the  radius  of  the 
eirth,  these  functions  can  be  substituted  by  their  values,  which 
refer  to  the  point  of  the  missile  launch,  and  the  derivatives  them¬ 
selves  can  be  calculated  on  the  assumption  that  the  earth  is  a  sphere 
with  a  radial  distribution  of  density. 

The  last  assumption,  if  one  considers  equations  (2.5*4),  leads 
to  the  following  expressions  for  the  mentioned  derivatives  of 
accelerations  of  the  force  of  gravity: 

cn.3.1) 

Here  C,  n  and  c  -  coordinates  of  the  arbitrary  point  of  space 
and  p  -  its  distance  to  the  center  of  the  earth. 


For  the  point  of  thd  missile  launch  (see  S  1  of  Chapter  II)  we 

have 

8-t-0  and  (4.3.2) 

and,  consequently,  at  this  point,  according  to  equations  (4.3.1) , 


Tr”-TT“~,,*>  *S‘"TT‘"0, 

=  (“.3.3) 


where  there  is  introduced  designation 


-in- 


(4.3.4) 


In  the  theory  of  gyroscopes  and  inertial  navigation  quantity 
v,  determined  by  equation  (4.3-4),  is  known  by  the  name  of  Schuler 

frequency  (v  :  C.  00125  l/sK 


Using  the  approximate  expressions  (4.3-3)  for  derivatives 
~5C  ’  ~S*  ”  IF  and  ir'  equalities  (4.2.2),  (4.2.9),  (4.2.14)  and 

(4.2.15),  we  arrive  at  the  following  simple  equations: 


*«(«)“• - •  4.  (a)  “  (v)  “•  0.  k„  (a)  —  — j—  . 

fU  (a)  ™ —  y  •  A*  (•)  “  (*)  “  0.  (4.3.5) 

and  also 

Ai  (a)  “  ~  •  h*  (°)  “  t*  (3)  “  0.  /*»  (v)  ™  -Jy  - 

*«(°)  *=  —  -y ,  •»*  («)  (o)  “  -y  ■  (4.3.6) 

Having  assumed  in  these  equations  o  *  o*,  one  can  present  expressions 
(4.2.7)  and  (4.2.11)  for  isochronal  variations  of  coordinates  and 
projections  of  velocity  of  the  missile  with  the  initial  quadratic 
approximation  (4.1.6)  in  the  following  final  form: 


«uo*=[i--^T-]*^(o, 

*«t  <0  -  4Ft  (0  -  ^  4^(0.  (4.3.7) 

te,(<)  =  *K,(l)  + 45,(0- 

Similarly,  for  the  case  of  the.  initial  approximation  of  functions 

and  6n(i)  in  the  form  of  polynomials  (4.1.7),  on  the  basis  of 
equations  (4.2.12)  and  (4.2.13),  we  obtain 

4$(0~[1-?^4*<0  +  ^4Kl(0, 

*1(0 -[i  b-v*(o*4k,(o, 

«*<*>-[*+  ^g^]^t(0-T*-a^(0. 

4-,(0  -  [l-  ^^]4F,(0  +  vV45,(/).  (4.3.8) 

Equations  (4.3.7)  or  (4.3.8),  which  approximately  consider  the 
change  in  action  on  the  missile  of  the  force  of  the  earth's  gravity 
' because  of  the  motion  of  the  missile  not  according  to  the  assigned 
calculated  law),  should  be  more  accurate  than  similar  equations 
(3*2.2)  and  (3.2.6)  given  in  Chapter  III. 

§  4 .  Construction  of  Ballistic  Equations  in  Which 
the  Effect  of  Changes  in  Acceleration  of  the 
Force  of  the  Earth's  Gravity  Is 
Approximately  Taken  Account 

In  the  previous  paragraph  equations  (4.3./  and  (4.3.8) 

were  obtained  for*  the  determination  of  variations  of  cooi’di nate -  arid 
projections  of  the  velocity  of  the  missile  with  a  ippr< :  xlmate 
calculation  of  the  effect  of  the  change  in  the  force  of  the  earth's 
gravity  with  motion  according  to  the  law  which  is  somewhat  distin¬ 
guished  from  the  calculated.  The  mentioned  equations  allow  more 
accurately,  in  comparison  with  expressions  '3-3*13'  and  '3*4.Jb> 


obtained  in  Chapter  III,  constructing  the  left  side  of  the  basic 
^ ^ ® t 1 ®  equation  (2.7. o)  aboard  the  missile,  using  readings  of 
the  came  integrators  of  accelerations.  For  this  we  will  use  in  the 
beginning  equations  (4.3.7)  and  substitute  into  equation  (2.7.6) 
expressions  for  6 i(t),  6n(t),  <5 u^(t )  and  6u^(t)  determined  by  them; 
we  obtain 


(1. 4.1) 


Further,  similar  to  that  which  was  done  in  §  2  of  Chapter  III, 
let  us  present  variations  of  projections  of  velocity  6v^(t)  and 
5F  (t)  in  the  form  of  differences  (3.2.3)  and  use  equations  (3.2.12). 
In  this  cast-  just  as  in  the  conversion  of  equation  (3-2.10),  terms 
standing  in  the  right  side  of  equality  (4.4.1)  are  reduced  with 
oimilar  terms  of  its  left  side.  Taking  into  account  still  equations 
i.w.  7/,  we  arrive  at  the  equation 


in  which 


i4(f)-<(oi  Ai + 1  j,(o-  £<om;+  + 

+  -  K  (O  .  (4.4.2 ) 

S  91  +(<ff  91  k*s*  ,V 

ai~7T - UT7i~-r  HZ' 

(4.4.3) 

A’  «  .  91  ,  2vV  91 

•  *  jh|  "*  6  9i\  **  J  9u%  ’ 


Equation  (4.4.2)  differs  from  equation  (3.2.13)  only  by 

$ 

■ot-f ficlents  of  isochronal  variations  of  projections  of  the  apparent 
path  and  Such  coefficients  in  equation  (3-2.13)  are,  I 

respectively,  partial  derivatives  —  and  and  in  equation  (4.4.2)  -  \ 


I 


are  quantities  and  A ^  distinguished  little  from  them.  Therefore, 
subsequent  conversions  of  equat  on  (4.4,2)  can  be  produced  exactly 
as  in  5  3  of  the  previous  chapter.  Then  instead  of  equation  (3.3.13), 
we  will  arrive  at  the  following: 

n  (*)  +  ✓  I<V  (0  -  (1)1  -  (4.4.4) 

The  left  sides  of  equations  (4.4.2)  and  (4.4.4)  are  constructed 
aboard  the  missile  with  the  help  of  the  same  technical  means  as  for 
equations  (3.2.13)  and  (3-3-13),  i.e.,  two  integrators  of  accelera¬ 
tions  and  the  appropriate  computer,  which  includes  elements  of 
repeated  integration.  In  this  case  the  axis  of  sensitivity  of  one 
of  the  integrators  is  at  the  same  constant  angle  X  to  axis  £  as  .1  r. 
the  construction  aboard  the  missile  of  the  left  side  of  equation 
(3.3.13).  As  regards  the  axis  of  sensitivity  of  the  other  Integrator, 
the  latter  should  be  inclined  toward  the  axis  f.  at  an  angle  of  u', 
which  is  little  distinguished  from  the  corresponding  angle  p.  Similar 
to  equality  (3.3.8)  the  magnitude  of  angle  o'  is  found  from  relations 

A^mm  N'coan’.  A\mm  V'llnp',  f  I4  _  U  .  c  ) 

where,  of  course, 

N'-y{Ait  +  (A‘,t  (4.4.6) 


In  turn  the  coefficient  p'  is  distinguished  little  from  coefficient 
y  of  equation  (3.3-13)  and  is  determined  by  equation 

(4.4.7) 

Here  quantity  M  is  the  same  as  that  in  equation  (3.3.4). 
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similar  to  the  form  of  equation  (3. 4.16).  The  use  of  equation 

(4.4.8)  allows  solving  the  problem  about  inertial  control  of  flight 

range  of  ballistic  missiles  just  as  in  S  4  of  Chapter  III  by  means 

of  one  meter  of  acceleration  with  an  alternating  direction  of  the 

axis  of  sens  i  ti  v ty  .  In  equation  (4.4.8)  the  variable  coefficient 

K'lti  is  formed  by  the  same  equations  (3*4.11)  as  coefficient  K(t) 

in  equation  (3-4.16)  but  with  the  replacement  of  angle  <(t)  by  angle 

<  ’  f  t )  and  derivatives  ■jj-  and  -|^  ,  respectively,  by  quantities 

and  A  '  .. 
n  ” 

The  use  in  the  transformation  of  the  basic  ballistic  equation 
(2.7.6)  of  more  complex  equations  (4.3.8)  for  coordinates  and  pro¬ 
jections  of  velocity  of  the  missile  leads,  naturally,  to  more  bulky 
results  in  -comparison  with  equations  given  in  the  beginning  of  this 
section.  Instead  of  equation  (4.4.1),  in  this  case  we  arrive  at 
the  more  complex  equation 

+«£.{[(  K+’,’’cj  + 

+^1*gi«+[t+aS£]«1+ 

+w4--S**+['-JTE]£)- 


In  turn  equation  (4.4.2)  is  replaced  by  the  following: 


4  +  U«(0-^(0l  A\  +  »'|(I)  Si  + 

+  r,(i)  B\ -  Vlit) (Bl-  £■) -  v\(‘){s\-Q  - 

(4.4.10) 


In  It  these  deslgnatior.t  are  introduced 


i  v  w  <t( 

JH 


7« 


(4.4.11) 


». 


•'(a*)*  9* 
15“*? 


+[*- 


I**’ 

v*(3,),'i  ai_ 

8  J  i%,  • 


Similarly,  equation  (4.  *1.4),  when  using  equations  (4.3.8),  is 
replaced  by  the  following: 

»V  <0  +  P"  IV  <0  -  K-  (01  =  V{.  (a)  + 

+  [(«;  -  a{  (o')  +  („•)] .  (4.4.12) 

Here  X"  and  p"  are  directions  of  axes  of  sensitivity  of  two  integra¬ 
tors  of  accele'  ations ,  which,  respectively,  differ  little  from  X- 
and  p-directions ,  introduced  in  i  3  of  Chapter  111.  Angles  X"  and 
p",  which  these  directions  from  with  axis  £,  are  determined  by 
relation 


Jtf'ooaX*, 

Coefficient  p"  in  equation  (4.4.12)  is  expressed  by  equation 


P‘ 


r 

IT 


Tv* 

(B\f  +  (B'J  ’ 


( 11 .  4  ,  x  4  ) 


The  construction  of  the  right  side  of  the  ballistic  equation  in 
the  form  (4.4.12)  aboard  the  missile  is  complicated,  as  compared  to 
the  case  of  the  equation  of  the  form  (4.4.4),  by  the  necessity  of 
introduction  into  the  computer  of  an  additional  term  linearly 
lenenJent  on  time.  Finally,  subsequent  conversions  of  equation 
(4.4.12),  specifically,  to  the  form  similar  to  the  ballistic  equa¬ 
tion  i  4. ,4. 8)  and  also  (3.5.2)  are  allowed. 

in  conclusion  let  us  note  that  the  selection  for  the  specific 
sv-'em  t.f  control  of  the  flight  range  of  one  of  the  ballistic 
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equations  given  in  this  and  previous  chapters  should  be  produced  on 
the  basis  of  estimates  of  quantities  of  the  so-called  systematic 
errors  inherent  to  the  system  of  inertial  control  of  the  given  type 
of  rocket.  Systematic  errors  are  understood  usually  as  those  errors 
in  the  determining  of  flight  range  of  the  missile,  which  appear 
exclusively  because  of  simplifications  made  in  the  derivation  of 
the  appropriate  ballistic  equation.  Because  of  this,  even  with  the 
accurate  switching  off  of  the  engine  of  the  rocket  at  the  instant 
Indicated  by  the  ballistic  equation  and  with  the  accurate  operation 
of  all  remaining  instr aments  and  devices  of  the  system  of  Inertial 
control,  the  actual  flight  range  of  the  missile  can  be  somewhat  dis¬ 
tinguished  from  that  assigned  according  to  the  calculation.  As 
regards  to  the  technical  difficulties  of  the  construction  of  systems 
of  inertial  control  of  flight  range  of  ballistic  missiles,  they  are 
included  in  the  manufacture  of  meters  and  integrators  of  the  apparent 
accelerations  (i.e.,  newtonmeters  and  impulse  meters)  with  extremely 
small  instrumental  errors  and  gyroscopic  instruments  with  minimum 
servicing  of  their  stabilized  axes  relative  to  directions  on  fixed 
star.  The  same  refers  to  the  immediate  switching  off  of  the  engine 
on  the  signal  of  the  achievement  by  the  ballistic  function  of  a  value 
■corresponding  to  the  assigned  flight  range  of  the  missile. 
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CHAPTER 


CONTROL  OF  LATERAL  MOTION  OK  THE  BALLISTIC  MISSILE 


§  1 .  Lateral  Deviation  of  the  Missile  from  the  Target, 
Expressed  In  the  Form  of  a  Function  of  Changes  In 
Parameters  of  the  End  of  the  Powered-Flight 
Section  in  the  Starting  System 
of  Coordinates 

In  the  solution  of  the  problem  about  the  elimination  of  lateral 
deviation  In  the  missile  from  the  assigned  target,  there  appear 
additional  alfficulties  connected  with  the  fact  that  cf  the  duration 
of  the  free-flight  section  of  real  motion  of  the  missile,  as  a 
rule,  is  somewhat  distinguished  from  its  calculated  value.  Tn  view 
of  this,  control  according  to  the  earlier  assigned  law  of  the  motion 
of  the  missile  in  a  direction  perpendicular  tc  the  programmed  plane, 
i.fe.,  plane  xy  of  the  starting  system  of  coordinates  xyz  or  (in 
another  possible  variant)  to  plane  £n  of  the  nonrotating  system 
(hC,  in  general  does  not  provide  the  absence  of  lateral  deviation 
in  the  missile  from  the  target.  Really,  because  of  the  laws  of 
mechanics  the  missile  on  the  greater  part  of  the  free-flight  section 
of  Its  flight  is  moved  with  minute  deviations  from  a  certain  plane, 
which  does  not  change  Its  orientation  relative  to  directions  at 
fixed  stars  (the  mentioned  deviations  are  connected  basically  wit* 
the  nonsphericity  of  the  earth).  Therefore,  if  the  duration  of  the 
free-flight  section  is  not  equal  to  its  calculated  value,  then 
because  of  the  rotation  of  the  earth  lateral  deviation  in  the 
missile  from  the  target1  continuously  appears.  However,  with 

'An  exception  is  the  case  of  flight  of  a  missile  in  the  plane  of 
the  equator  and  from  one  pole  on  another. 
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small  distinction  in  the  actual  motion  of  the  missile  from  cal¬ 
culated  this  deviation  should  also  be  small.  Correspondingly  the 
control  of  motion  of  the  missile  in  the  powered-flight  section  in 
a  lateral  direction  can  be  reduced  only  to  the  requirement  of  the 
fulfillment  of  an  earlier  assigned  (calculated)  law  of  motion  in 
the  projection  on  the  z  axis  (or  in  another  variant  -  on  axis  c). 
However,  for  a  more  accurate  elimination  of  lateral  deviation  in  the 
missile  from  the  target  it  is  necessary  to  take  into  account  how 
at  current  instant  of  time  its  real  motion  in  a  longitudinal 
iircctior  is  distinguished  from  the  calculated.  Por  this  Durpose 
one  can  use  current  readings  of  instruments  of  the  control  of 
flight  range  of  the  missile. 

It  is  obvious  that  an  especially  accurate  control  of  lateral 
motion  of  the  missile  is  necessary  only  during  a  short  interval  of 
time  when  cutoff  of  its  engine  can  occur  in  flight  at  the  assigned 
range.  At  the  remaining  time  of  the  powered-flight  section  the 
accuracy  of  control  of  lateral  motion  can  be  somewhat  less.  The 
duration  of  the  mentioned  interval  of  allowed  cutoffs  of  the  engine 
depend  on  the  quality  of  control  of  motion  of  the  missile  according 
to  the  assigned  program. 

In  the  starting  system  of  coordinates  xya  (see  5  1  of  the 
second  chapter.  Pig.  9)  quantity  b  -  of  .lateral  deviation  in  the 
missile  from  the  target,  just  as  the  range'  l",'  is  a  function  of 
coordinates  x,  y,  z  and  projections  v  ,  v  .  v  of  its  velocity 

xya 

relative  to  this  system  of  coordinates  at  the  instant  of  termination 
of  the  powered-flight  section  of  the  flight.  Thus, 

*-*(*,  p,.*,).  (5.1.1) 

The  astronomical  time,  which  corresponds  to  instant  of  termina¬ 
tion  of  the  powered-flight  section,  does  not  play  any  role  here  as 
in  equation  (2.1.2),  unless  one  does  not  take  into  account  the 
negligible  effect  of  such  factors  as  the  mutual  arrangement  of  the 
earth  and  moon  or  the  earth  and  sun. 
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With  the  replacement  of  arguments  of  function  (5.1.1)  by  their 

calculated  values  x *,  y* ,  z* ,  y*  v*  and  v*  it  should  turn  into 

x  u  z 

zero,  since  the  lateral  deviation  in  this  case  is  absent.  Thus, 

j  9 

=  o  (5.1.?) 

and,  consequently, 

A&  =  b(r,  y,  z.  p„  p„  p,)-6(*\  y\  z\  pj,  p*)  *  b,  (5.1.3) 

i.e.,  the  error  in  the  lateral  deviation  Lb  is  equal  to  the  lateral 
deviation  b . 

In  accordance  with  the  Taylor  series  for  the  function  of  many 
variables,  we  have,  correct  to  smalls  of  the  first  order  inclusively 
the  following  expansion: 


b  -  (*-,*)  £  +(y-  y*) £  +  (s  -  -)  £  + 

+{*«  -  <>  £  +  (\  -  v\)  %  +  («Y -  O  S:  • 


(5.1.^ 


Here  derivatives  ,  -fr ,  4r~  .  -7—  11  '  -|r  'xr’r  *aken  at  cal- 

culated  values  of  their  arguments  and,  consequently ,  for  each 
calculated  case  of  flight  of  the  missile  are  earlier  known  quantitie 
They  can  be  determined  in  a  way  similar  to  the  ball!  --1 1-.  inefficient 
giver:  in  §  1  of  the  second  chapter. 

5  ? .  Lateral  Deviation  as  a  Function  of  i grumcler. 
of  the  End  of  the  Powered-Flight  .Section  in 
the  Nonrotating  System  of  Coordinates 

In  examining  the  problems  of  inertial  control  '  ■.  feral 

motion  of  the  ballistic  missile,  the  starting  system  ■'  .'.vrd  i  iirv 

is  Inconvenient.  Therefore,  let  us  Introduce  the  same  nor.-:  *  r  '  :.g 

system  of  coordinates  as  that  in  the  examining  of  the  tht.r;- 

'nertial  control  of  range  (see  5  2  of  the  second  chapter).  Current 

coordinates  x(t),  y(t),  z(t)  and  projections  v  (t),  v  (t)  and  v„(t) 

x  y  z 

.T  trc  velocity  of  the  missile  relative  to  the  starting  system  of 


coordinates  xyz  are  expressed  in  this  case  by  current  coordinates 
ZU),  n(t),  t,(t)  and  projections  u^t),  u^(t)  and  u  (t)  of  the 
velocity  of  the  missile  in  the  nonrotating  system  of  coordinates 
SnC  by  equations  (2.2.5)  and  (2.3.5).  Entering  into  the  mentioned 
ecuations  also  in  evident  form  (through  angle  <p  «  ut )  is  time  t,  which 
passed  at  the  instant  of  launch  of  the  missile.  In  view  of  this, 
in  accordance  with  equation  (5.1.1),  the  magnitude  of  lateral 
deviation  in  the  missile  from  the  target  can  be  represented  in  the 
form  of  a  function1 

«c:a)  (5.2.1) 

of  coordinates  5,  n,  b  and  projections  u n,  u of  the  velocity 

of  the  missile  relative  to  the  system  of  coordinates  at  the 
instant  t  =  a  of  termination  of  the  powered-flight  section  and  also 
its  very  duration  o. 


Calculated  values  of  coordinates  x»,  y*,  3*  and  projections 
J’y*  ’*  ~  ^  velocity  of  the  missile  relative  to  the  starting 
system  xyz  correspond  to  certain  calculated  values  of  coordinates 
n«,  and  projections  u*,  u»,  u »  of  velocity  of  the  missile 
in  the  nonrotating  system  CnC.  It  Is  obvious  that  function  (5.2,1), 
which  represents  the  lateral  deviation  of  the  missile  from  the 
target,  will  turn  into  zero  at  these  values  of  their  arguments  if, 
furthermore,  we  assume  the  duration  of  the  powered-flight  section 
is  also  equal  to  its  calculated  value  a*.  Thus, 

«  *1  •  C  >  wj;  tf*)  *0.  (5.2.2) 


Expanding  function  (5.2.1)  In  Taylor  series  and  being  limited 
in  the  latter  only  by  terms  of  the  first  order  of  smallness,  we 
obtain,  taking  Into  account  equality  (5-2.2),  the  equation 


b  -  <*■ -  n  £ + oi.-  »,•)£ + (;- c-) »  + 

(5.2.3) 


‘Here,  Just  uj  In  §  ^  of  the  second  chat  tern  ’  mat  ’ ono 
^2.^.1)  are  -v:  cep  ted  .  ' 
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in  which  partial  derivatives  -5^.^,  jg-.  and  -jj  are 

constant  coefficients,  since  their  arguments  £,  n,  ?,  un, 
and  o  are  taken  equal  to  their  calculated  values,  l.e.,  respectively 
£*,  n#,  ?•,  ujf,  u*,  u*  and  a*.  Magnitudes  of  these  -nef f 1 r l ents 
depend  finally  on  the  selection  of  a  certain  calcv.'.ated  case  of 
flight  of  the  missile.  In  a  way  similar  to  the  appropriate  partial 
derivative  of  function  Z,  which  represents  the  flight  range  of  the 
missile,  they  are  expressed  by  equations  of  the  type  ‘  Z  .  .  5 )  and 

(2. *4.6)  in  terms  of  partial  derivatives  ■£,  J7 .  .  ■£-, 

and  also  quantities  x* ,  y*,  a*,  u»  u* ,  ujj. 

Equations  (5.2.3),  correct  to  squares  and  products  of  differences 
5—6*.  «|  —  n\  C  —  6*.  “t  —  **t.  «»  —  “*,  nj  —  mJ,  «  —  a*  ,  determines  the  lateral 
deviation  in  the  missile  from  the  target  when  magnitudes  of  coordi¬ 
nates  and  projections  of  velocity  of  the  missile  at  the  end  of  the 
powered-flight  section  are  not  equal  to  their  calculated  values. 


magn 

deri 


With  real  motion  of  the  missile  neither  the  instant  of  ^witching 
off  of  the  engine  o  nor  its  coordinates  £,  n  and  projection  of 
velocity  u^,  u^,  which  correspond  tc  this  instant  of  time,  are 
known  earlier.  Therefore,  for  elimination  with  the  abovp-mentl oned 
accuracy  of  lateral  deviation  in  the  missile  from  the  target,  a 
continuous  change  Is  necessary  in  quantities  Ci t)  and  u  ft:  charac- 
terizing  it.-  motion  in  a  lateral  direction  during  a  certain  time 
interval  of  the  end  of  the  calculated  powered-flight  section.  The 
purpose  of  such  a  control  is  the  continuous  reduction  to  zero  of 
funct ion 


makl. 


mlssl. 


p(o- 16(0  -  n  £  +  iu<o~  vi  £  +  (C(o-  ;*(  £  + 
+  K<0~  «;ii£  +  i\<o  -  «;i  £  + 


+  i*«<o-u;i  £+<<-*•)£ 


<  ^ .  ;■ .  u 


ccnsi 
and  u 


i  1.  1 1  ■«.  wh.'ie  .  r.terval  of  possible  instants  of  the  switching  off  f 
ttic  engine  for  a  missile  of  the  giver,  system.  Actually,  function 


f  unc  1. 1 
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6 (t),  according  to  equations  f?.1*.].)  and  (5.2.31.  turns  into  the 
magnitude  of  lateral  deviation  in  the  missile  b  if  we  assume  in 
it  argument  t  equal  to  the  actual  value  of  the  Instant  of  the 
switching  off  of  the  engine  a. 

The  current  value  of  function  Q(t),  which  we  will  call  the 
function  of  lateral  deviation ,  can  be  used  as  the  controlling 
signal  in  the  system  of  control  of  lateral  motion  of  the  missile. 
This  system  should  operate  in  such  a  way  that  function  6(t'  would 
turn  into  zero  because  of  the  corresponding  change  in  current 
magnitudes  of  the  coordinate  of  the  missile  c(t)  and  its  time 
derivative  u^(t). 

§  3 .  Construction  of  Corresponding  Terms  of  the  Pune  cion 
of  Lateral  Deviation  Dependent  Only  on 
the  Basic  Motion  of  the  Missile 


As  was  shown  in  the  previous  section,  for  the  inertial  control 
of  lateral  motion  it  is  necessary  to  reproduce  the  current  value 
of  function  8  (i)  aboard  the  missile.  The  totality  <~f  terns  <  5 .  2  .  "  '> 
making  up  this  function  is  decomposed  into  two  group?  6^  t)  and 
6  :'t>.  The  first  of  them  is  the  sum 

Mo  -  iKo  -  ri  £  4-  wo  -  htHJ  + 

(o  — ;i  £ + i-,(o  -  «;i  £  +«-  <o£.  (5.3.1 , 


the  magnitude  of  which  is  determined  by  the  basic  motion  of  the 
missile,  i.e.,  by  the  motion  of  its  projection  on  plane  ;.n  of  the 
nonrotating  system  of  coordinates  £nc.  The  second  group 

M4-1C<4  — +  !  5.  " 

consists  cf  two  terms  containing,  respectively,  quantise  C'f' 
and  u  ( t  . 

Some  method?  if  the  construction  if  the  current  value  f 
function  6. it  >  aboard  the  mlssl.*-'  are  examined  in  the  >'<  l  owing 

bb 


section.  As  regards  the  stun  (5*3. 1),  then  it,  according  to  its 
structure,  is  completely  similar  to  equation  (2.  *1.9)  for  the  initial 
ballistic  function  e(t).  Therefore,  the  approximate  construction 
of  the  sum  (5.3.1)  can  be  produced  aboard  the  missile  by  the  same 
methods  as  were  used  for  the  construction  of  function  e(t).  Let 
us  examine,  for  example,  the  construction  of  this  sum  aboard  the 
missile  by  means  of  using  readings  of  two  integrators  of  accelera¬ 
tions,  the  axes  of  sensitivity  of  which  are  oriented  respectively 
parallel  to  axes  £  and  n  of  the  nonrotating  system  of  coordinates 
£n£.  Similar  to  equation  (2.7.*0,  the  sum  (5.3.1)  can  be  represented 
in  the  form 


Pi  (<)  -  (g  (0  -  g*  (i)j  £  +  [n  (/)  -n‘(0l£  + 

+ Km  - 1 »;  (oi  ^ + («,  to  - <01  £ + 

•  +[s*(o-n-£+ + i*i<o -■&£'+ 

+  1*»*(0  —  + <*"■*)  "S'-  (5.3-3) 


With  subsequent  conversion  of  the  sum  B-^(i)  we  will  use  equations 
(2.7.1)  and  (2.7.2)  for  isochronal  variations  of  coordinates  and 
projections  of  the  velocity  of  the  missile  and  also  equations  of  the 
expansion  in  Taylor  series  (being  limited  only  to  terms  of  the 
first  order  of  smallness) 


r  to -g*-g-  (o -g-(o-(<- a*)  &P- , 

n*  (*) — n*-n*(0— n*  (« )=■<<- o')  , 


(5. 3- *) 


and  in  exactly  the  same  manner 

(0 -t/j -«•(/)■  u;  (o*>  -  ( t  -  o')  - , 

.  *  .  *. 

«;(0-«;(a*)-(<-3’)~5^-  (5-3.5) 

Taking  into  account,  furthermore,  that  in  accordance  with 
(2.6.3) 


(5.3.6) 


_  -n‘(r)  <*«>’> 

.  *  *  ~ E  & “V“  • 


we  transform  the  expression  for  sum  $^(t)  to  the  form 


P.  (0  =  **  (0  #  +  *1 <0  £  +  «"t  w  £  +  fcr, «  ~  + 

+*-*tw+*spv 

.  **•<«*)  »  .  •PtiVa*)  3*  ,  an 

+  “E*- 5^  + -7?r- 5T,  +  a?  J  • 


(5.3.7) 


Just  as  in  the  third  chapter,  let  us  disregard  in  the  beginning 
the  effect  of  the  change  in  the  acceleration  of  the  force  of  gravity 
during  motion  of  the  rocket  not  according  to  the  calculated  law. 

In  accordance  with  equations  (3.2.6),  (3.2.7),  (3-2.2)  and  (3.2.3), 
we  obtain  at  small  values  of  the  difference  t  -  a*  relations 


H<0  -  4*<<)  ■=  -  MO-  + 

+  (<-«>*';  <«*>1.  (5.3.8) 

d>)  (0  -  W,  (0  -  5,  (/)  -  S\  (0  -  S,  (0  -  [£>  («*)  + 

+  (<  —  «*)  K  (a’)l» 

and  also 

«M0  -  «\<0  =  Ft(0-  Fj(0  =  Ft(l)-  |K(a  )  + 

+  (<-o')«i(3')l. 

flu,  <t)  -  IV,  (/)  -  V,  (t)  -  K  (0  =  V*  (0  — iK  (s’)  + 

+  (*  —  »*)««(<**)).  (5.3.9) 

Used  here  again  are  equations  of  the  expansion  in  Taylor  series, 
but  for  projections  s{  (t).  Si  (0  of  the  apparent  path  of  the  missile 
arid  projoctionr.  !'•(/).  1"Jj  (t)  of  Its  app^^ent  velocity  in  the  calculated 
motion  equations  (3-2.5)  and  (3-2.9)  are  taken  into  account. 

By  means  of  equalities  (5.3*8)  and  (5-3-9)  the  sum  (5.3.7) 
is  reduced  to  the  form 


fc(«>  -  ■%.<<)  £  +  <5.(0 £  +  V K <i>£  +  K, </)  ± _ 

~  [•* <**>  *& + i*)  £ + vi  (3‘>  ^ + ri  (*’>  /y + 

+  <'-’*>  {[^- *<*•>]£  +  l*P- ^*>]£-f 
+ psp  _  **•>]  -  +  pp»  -  <o]  A + 5} . 


(5.3.10) 


In  the  last  equation  we  can  produce  further  simplifications,  since, 
according  to  equality  (2.6.1),  we  have 

~ *1  (0  “  Alt  (0.  *»  W;  *1. 

V«>;  0  (5.3.11) 

as  a  consequence,  as  a  result  of  the  integration  with  respect  to 
time  within  limits  of  J  ■  0  to  t  ■  c*,  allowing  for  equations 

(3.2.5) 


— jP’  -  vt  <«■)  -  ‘  A  It*  <0.  *f  (i),  1 1  <u  + 

« 

^  («*)  -  j  A  It* «).  n*  <»).  'J  dt  +  , 


Equations  (5. 3 .11)  and  (5.3.12)  can  be  used  In  the  subsequent 
conversion  of  expression  (5.3.10)  for  the  sum  6-^(t).  An  a  result 
let  us  obtain  the  following  final  expression  for  this  sum: 

fc  -  *(*)  £  +  s%{ o£  +  rt</>£  +  »mo£- 
-.C  +  (i-o’)0.  '  (5.3.13) 

Introduced  here  are  designations 
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D  =  £  *' /i  If  (0,  *r (0;  <I*  +  5^M*,(0,V («): »l *H- 

•  ft 

+  /l  If  («•).  V  (3*).  yj  £  +  /,  If  (3*>.  <1*  (3’);  3*1  £-  + 


(5.3.15) 


It  is  obvious  that  C  and  D  are  constant  quantities,  which  can  be 
determined  earlier  for  each  specific  case  of  flight  of  the  missile. 

For  the  construction  of  expression  (5.3.13)  aboard  the  missile, 
one  can  use  the  same  integrators  of  accelerations  and  elements  of 
repeated  Integration  as  in  the  solution  of  the  problem  on  inertial 
control  of  flight  range  by  the  method  given  in  I  2  of  the  third 
chapter.  Furthermore,  an  additional  computer  and  clock  mechanism 
are  necessary  here. 

If,  however,  the  control  of  flight  range  of  the  missile  is 
produced  with  the  help  of  integrators  whose  axes  of  sensitivity 
are  parallel  to  X -  and  u-directions,  then  for  the  formation  of 
expression  (5.3.13)  there  will  be  required,  furthermore,  an  auxiliary 
device  of  the  type  of  converter  of  coordinates,  which  continuously 
solves  the  system  of  equations  (see  §  3  of  the  thlra  chapter) 


Kill) co. i  +  F,(/)sln  \  -  K,(/), 
Vi(J)eos|»  +  F,(l)siu|i- V,(«) 


(5.3.16) 


relative  to  quantities  V^(t)  and  ^(t)  according  to  data  of  current 
readings  v.(t)  and  k  (t)  above-mentioned  integrators. 

A  p 

§  H .  Construction  of  Main  Terms  of  the  Function 
of  Lateral  Deviation 


?he  second  group  of  corresponding  terms  of  equation  (5-2.1*) 
for  function  Bit),  as  was  indicated  in  the  previous  section,  is 


the  si:.: 


two  terms 


ICIO-n-J  and 


(5.H.1) 
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They  can  be  called  the  main  terms  of  the  function  of  lateral 
deviation  B(t),  since  they  play  the  main  role  in  the  control  of 
lateral  motion  of  the  missile. 

For  the  construction  of  current  values  of  the  sum  (5.4.1) 
aboard  the  missile,  one  can  use  readings  V^(t)  of  the  integrator 
of  accelerations,  the  axis  of  sensitivity  of  which  is  oriented  in 
parallel  to  axis  C  of  the  nonrotating  system  of  coordinates 
(so-called  lateral  integrator). 

Let  us  note  that,  according  to  equation  (1.4.1)  (S  4  of  Chapter  I) 
there  occurs  the  relation 

-*(<)  + A.  (5.4.?) 

where  a^(t)  -  projection  of  apparent  acceleration  of  the  missile 
on  axis  c,  and  /  -  projection  on  the  same  axis  of  acceleration  of 

the  force  of  the  earth’s  gravity.  Th£  latter,  similar  to  equations 
(2.5.4),  with  sufficient  accuracy  can  be  represented  in  the  form 

(5.4.3) 

Thus,  relation  (5.4.2)  can  be  considered  the  differential  equation 
in  the  desired  function  C(t). 

At  the  instant  t  =  0  the  axes  of  the  starting  system  of  coordi¬ 
nates  xyz,  which  is  connected  to  the  earth,  respectively  coincide 
with  axes  of  the  nonrotating  system  CnC,  and  the  center  of  masses 
of  the  missile,  on  assumption,  is  found  at  their  common  origin. 
Furthermore,  at  this  instant  the  missile  does  not  have  velocity  with 
respect  to  the  st*~l 1 nm  system  of  nates  xun.  Corseauently , 

initial  conditions  of  equation  (5.4.2)  are  such: 

1(0)- 0,  ^p.-«c(0  )-u».  (5.4.4) 

where  -  velocity  of  the  beginning  of  the  starting  system  of 
coordinates  in  the  nonrotating  system  at  the  instant  when 


t 


t  -  0,  i.e.,  at  the  beginning  of  the  powered-flight  section  of 
its  flight. 

It  is  possible  to  assume  that  the  acceleration  of  the  force 
of  gravity  proves  to  have  an  insignificant  effect  on  the  rate  of 
change  in  function  l(t).  Therefore,  as  a  first  approximation  it 
is  possible  not  to  consider  the  projection  /  in  equation  (5.4.2) 
and,  respectively,  consider  that 

(5.4.5) 

Hence,  taking  account  the  second  initial  condition  (5.4.4),  we  have 

(/)  +  *.  (5.4.6) 

Here 

< 

r«(0-J«c(0*  (5.4.7) 

is  the  current  value  of  the  projection  on  axis  c  of  the  apparent 
velocity  of  the  missile  relative  to  the  nonrotating  system  of 
coordinates  This  Is  the  reading  of  the  integrator  of  accelera¬ 

tions  with  the  axis  of  sensitivity  parallel  to  the  axis  £.  Secondary 
integration  of  both  parts  of  equality  (5.4.5),  allowing  also  for 
the  first  of  the  initial  conditions  (5-4.4),  leads  to  the  equation 

C<0 -*(<>  +  •*».  (5.4.8) 

in  which  S  (t)  -  projection  on  axis  C  of  the  vector  of  the  apparent 
path  of  the  missile  in  the  nonrotating  system  of  coordinates  £nC. 

In  accordance  with  the  last  equation,  for  the  construction 
aboard  the  missile  of  current  values  $( t ),  an  additional  Integration 
rf  readings  of  the  aforementioned  integrator  of  accelerations  Is 
necessary.  Actually 

I 

& 


(5.4.5) 


On  the  basis  of  equations  (5.4.6)  and  (5.4.8)  the  expression 
for  the  sum  (5.3.?)  is  reduced  to  the  form 

MO  -  IMO  +  «v  -  n  $  +  ,r;(0  + 1-  •  (5.4.10) 

t 

This  sum,  Just  as  sum  8^(t),  obtained  in  §  3  of  this  chapter, 
can  be  constructed  directly  aboard  the  missile. 

Thus,  under  the  assumptions  made  about  the  unimportance  of 
quantity  /  ,  it  is  possible  by  means  of  integrators  of  accelerations 
and  computers  to  obtain  current  values  of  the  function  of  lateral 
deviation 


PW-MO  +  MO.  ‘  (5.  ,4. 11) 

§  5.  Calculation  of  the  Acceleration  of  the  Force 
of  Gravity  in  the  Construction  of  the  Function 
of  Lateral  Deviation 

Methods  given  in  the  fourth  chapter,  allow,  as  will  be  shown 
below,  approximately  integrating  equation  (5.4.2)  for  function 
5(0  taking  into  account  the  projection  of  the  acceleration  of  the 
force  of  the  earth's  gravity  /  . 

Let  us  note,  first  of  all,  that  with  an  accuracy  sufficient  for 
practice  at  not  too  great  an  extension  of  the  powered-flight  section 
of  the  flight  of  the  missile  in  equation  (5.4.3),  the  distance  p  from 
the  missile  to  the  center  of  the  earth  can  be  substituted  by  the 
radius  of  the  earth  R.  After  this  let  us  arrive  at  equation 

<0.  (5.5.D 

in  which  quantity  v  (Schuler  frequency)  is  determined  by  equation 
(4.3.4).  Finally  equation  (5.4.2)  can  be  replaced  by  the  following: 

v*C(/).  (5.5.2) 
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Let  us  integrate  wltii  respect  to  time  the  right  and  left  sides  of 
this  equation  within  limits  of  0  to  t  taking  into  account  the  second 
initial  condition  (5.^.4).  Equating  the  obtained  results,  we  arrive 
at  integral  differential  relation 

Vt(0  —  +»£•  (5.5.3) 

Let  us  produce  the  sam^  operation  of  integration  again  but  taking 
into  account  the  first  initial  condition  (5.^.**).  As  a  result  we 
obtain  the  integral  relation 

i  t 

C<0**  +  «#,  (5.5.JO  ; 

which  is  the  equation  for  the  determination  of  lateral  shift  cf  the 
missile  ift)  according  to  the  assigned  function  S^(t).  The  expression 
represented  by  equation  (S.^.B)  can  be  considered  the  first  approxi¬ 
mation  to  the  solution  of  this  equation.  To  obtain  the  following 
approximation,  in  which  already  the  effect  of  acceleration  of  the 
force  of  gravity  will  be  taken  into  account,  it  is  possible  to  use 
the  method  given  in  the  fourth  chapter.  First  of  all,  let  us  assume 
in  relation  (5.5.10  t  »  o.  Then  it  takes  the  form 

•  t 

t<«)  -  ■$<(< j)-  v*  j  j  {(<)**  +  (5.5.5) 

Let  us  further  replace  in  the  right  side  of  equation  ( 5  •  5 .  *0  function 
Z(tJ  by  its  approximation  in  the  form  of  a  polynomial  of  the  second 
power 

+  (5.5.6) 

the  coefficients  of  which  are  selected  so  that  initial  conditions 
(5.^. 4)  would  be  satisfied,  and,  furthermore,  in  order  that  the 
polynomial  would  take  value  c(o)  at  the  instant  t  *  o  of  termination 
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of  the  powered-flight  section  of  real  mtiou  of  the  missile.1 

if  then  we  produce  a  similar  replacement  also  in  relation 
(5.5.3)  and  accomplish  the  necessary  operations  of  integration, 
then  we  will  arrive  at  equality 


*»t  (a)  «=  Vc  (a)  —  v*u«  —  v»t  (a)  -i  +  «£, 
C(a)  — 5t(a)— v*t  (a)-^--f  u»a, 


(5-5.7) 


which  can  be  examined  as  the  equation  for  determining  the  value  of 
the  function  itself  c(a)  and  its  time  derivative  u^(o)  at  the 
instant  t  »  o. 


The  solution  of  the  second  equation  15.5.7)  can  be  obtained 
with  sufficient  accuracy  i f  in  its  right  side  quantity  5(0)  is 
substituted  by  its  first  approximation,  according  to  equation 
(5.4.8),  i.e..  If  we  assume 


t(o)-^(«)  +  «i»a.  (5.5.8) 

As  a  result  we  obtain  equation 

C<«) -  (l - +  (<- ( 5 .  ‘ .  '•  • 

We  produce  the  same  replacement  in  the  first  equation  of  (3- 5.'?). 
Finally  we  arrive  at  the  following  approximate  equation  for  the 
determination  of  the  pro,1ecti?©n  of  the  actual  velocity  of  the 
missile  on  ax  if  ■' 

“c(<*>=  Kc(a)— (a)  +  (l—  (5.5.10 

Considering  the  argument  a  to  be  close  to  0*,  it  is  ^ 00s Abie 
without  great  loss  of  accuracy  in  equations  (5.5.9)  and  '5.5.10) 


1 7h 'Te  is  a  possibility,  similar  to  §  1  of  the  fourth  chapter,  if 
nr.  approximation  by  means  of  a  polynomial  of  the  third  rower  ir 
which  t he  derivative  when  t  =  c  still  turns  into  u^(o). 


I 


to  substitute  argument  a  in  terms  containing  the  small  factor  v2, 
by  the  rated  value  of  duration  of  the  powered-f light  section  o*. 
Substituting  in  other  places  of  these  equations  argument  o  for  the 
current  time  t,  we  obtain  the  following  final  expressions  for  c(t) 
and  u  At) 


(5.5. 

liquations  (5.5-11)  lead  to  a  presentation  of  the  sum  B2(t),  more 
accurate  in  comparison  with  expression  (5. 4.10),  namely: 


M'> .  Ji_  -  Sp-H  -  C-)$  + 

+{n«)-^j,(o+[i-  spjs-  "-:}£■■■ 


15.5.U) 


The  construction  of  function  8 ?(t)  aboard  the  missile  can  be 
produced  by  means  of  the  integrator  of  acceleration  with  axis  of 
sensitivity  parallel  to  axis  C  (and  also,  of  course,  the  integrating 

devi ce ) . 


»  «  * 

I r.  this  book,  as  was  already  Indicated  in  the  introduction,  an 
account  of  problems  of  the  theory  of  inertial  /octroi  c”  *  hr  flight 
f  ballistic  missiles  was  produced  under  the  a -sumption  of  the  ideal 
tabilizaticn  of  axes  of  sensitivity  of  '  ,-wtc  •- r.eters  and  integrators 
of  accelerations  and  completely  accurate  operation  of  the  latter, 
accurate  indication  by  them  of  current  values  of  projections 
cf  the  vector-  jf  apparent  acceleration  or  the  missile  and  integrals 
of  t:  .i  vaui-  .  . 

An  estimate  of  the  effect  of  instrumental  errors  of  gyroscopic 
devices  and  integrators  of  accelerations  on  the  magnitude  of  possible 

fissile  with  flight  at  ‘  he  issighcd  t  -.  ..r- 


do-.  I  at  tor. 


t  0nd  also 


an  analysis  of  causes  of  the  appearance  of  these  errors  and  develop¬ 
ment  of  measures  on  their  elimination  represent  an  independent 
interest  and  fall  outside  the  framework  of  this  book. 

Let  us  note,  finally,  that  problems  of  inertial  control  of 
space  missiles  can  be  solved  by  the  same  methods  as  the  problem 
examined  above  about  the  flight  of  ballistic  missiles  within  limits 
of  the  earth. 


APPENDIX 

DERIVATION  OF  EQUATIONS  FOR  COSINES  OF  ANGLES  BETWEEN 
AXES  OF  TWO  SYSTEMS  OF  COORDINATES  TURNED  RELATIVE 
TO  ONE  ANOTHER  AT  A  FINITE  ANGLE 

In  manuals  on  analytical  geometry  and  according  to  theoretical 
mechanics,  equations  for  the  cosines  of  angles  between  axes  of  two 
systems  of  coordinates  arbitrarily  turned  relative  to  each  other  are 
derived  usually  with  the  help  of  several  additional  geometric  con¬ 
structions.  The  derivation  of  these  equations  is  also  possible 
cy  means  of  the  frequent  use  of  tables  of  cosines  of  angles  between 
axes  of  the  basic  and  auxiliary  system  of  coordinates,  each  of 
which  is  turned  relative  to  the  previous  around  one  of  the  coordi¬ 
nate  axes1  common  with  it.  Specifically,  the  main  axis  of  rotation 
is  simultaneously  the  axis  of  coordinates  of  two  auxiliary  systems 
turned  one  relative  to  the  other  at  the  same  angle  as  the  main 
axes  . 

Given  below  is  a  purely  analytical  derivation  of  the  mentioned 
equations,  founded  exclusively  on  the  simolest  theorems  of 
analytical  geometry,  and  some  of  these  theorem.-  are  used  in 
vectorial  form. 

I.ct  ui  assume  that  the  system  of  coordinates  is  turned  at 

an  arbitrary  angle  $  around  a  certain  axis  d,  which  passes  through 

'See,  for  example,  Whittaker.  Analytical  Kinetics.  M.-L.,  1937 
and  A.  Yu.  Ishlinskiy.  Mechanics  of  Gyroscopic  Systems.  I zd-vo 
AN  RASP  ,  1 9 (-.  3  . 
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the  origin  of  this  system.1  The  final  position  of  the  system  is 
designated  by  xyz.  If  Z,  m  and  n  are  osines  of  angles  which  form 
the  mentioned  axis  of  the  final  turn  d  with  axes  of  the  system  of 
coordinates  CnC,  then  it  is  obvious  that  the  same  magnitudes  are 
respectively  cosines  of  angles  between  the  axis  d  and  axes  of  the 
system  xyz.  Specifically,  we  have 

J  —  cos  a,  (1) 

where  a  -  angle  which  forms  the  axis  of  rotation  d  with  the  axis  £ 
of  the  system  of  coordinates  Ch?  and  simultaneously  the  same  axis 
of  rotation  d  with  the  x  axis  of  the  system  xyz  (see  Pig.  I1*). 


Let  us  designate  the  direction  cosines  of  the  x  axis  relative 
to  the  system  of  coordinates  Ch?  respectively  by  a ,  b  and  o.  They 
are  in  this  case  the  desired  quantities.  Besides  the  apparent 
equality 

«*  +  **  +  e*«=i  (21 


' In  §§  2  and  3  of  the  second  chapter  of  this  book  such  an  axis 
va.i  the  axis  of  rotation  of  the  earth,  and  angle  i  was  the  product 
if  the  angular  velocity  of  the  earth  U  for  time  t,  which  passes  from 
that  instant  when  the  axes  of  systems  of  coordinates  and  xyz, 

rvTOv  ot  ivoly  ,  coincided  . 


v 


( 


9 


the  direction  cosines  a,  b  and  a  satisfy  also  the  relation 

la  +  mb  +  ne  =  /,  ( 3 ) 

which,  in  accordance  with  the  known  equation  of  the  analytical 
geometry,  expresses  the  magnitude  of  the  cosine  o£  angle  a  between 
the  x  axis  and  axis  of  rotation  d. 

For  the  determination  of  the  three  desired  quantities  a,  b  and  <2, 
one  equation  is  necessary.  Specifically,  it  can  express  the  fact 
that  plane  xd  is  turned  relative  to  plane  £d  about  the  axis  of 
rotation  d  at  the  assigned  angle  Thus  the  measure  of  the  dihedral 
angle  between  the  mentioned  planes  also  should  be  equivalent  to 
angle  4>. 


Since  the  dihedral  angle  is  measured  by  the  angle  between  two 
perpendiculars  to  the  planes  forming  it,  then  we  will  introduce 
two  vectors  h  and  Tc ,  the  first  of  which  is  perpendicular  to  the 
piano  £<?,  and  the  other  -  plane  xd.  According  to  the  property  of 
the  scalar  pr<_  :uct  of  two  vectors,  we  have 

hkcoav- Mt  +  Mt  +  A-A.  (If) 

where  h^a  h^,  h ^  -  projections  on  the  axis  £,  n ,  t  of  vector  Ji 
and,  respectively,  k^,  k ^  -  projections  on  these  axes  of  vector 

k. 


Let  us  take  as  vector  h  che  vectoriai  product  of  the  single 
—1 n  _n 

vector  £'  ,  located  on  the  £  axis,  and  single  vector  d  ,  which 

possesses  the  direction  of  the  axis  of  rotati or  d.  Projections  of 

the  first  of  them  on  axes  £,  r)  and  c,  of  course,  are  equal, 

respect.!  vely  ,  to  numbers 

0,  0,  ( 5 ) 


and  of  the  second 


l,  m,  ». 


According  1  '  !  .•  ■■1'  ih-  composition  of  the  vector  ial  product, 

obtain 


and,  consequently. 


I  —  frx 


t» 

t  0  0 

I  M.  It 


(7) 


\-m.  (8) 

In  turn  for  vector  k  we  take  the  vectorial  product 


i 


l<  i  t  • 

f  .  MM 


(9) 


Here  the  second  line  of  the  determinant  is  made  up  of  projections 
of  the  single  vector  on  the  axes  £,  n  and  c,  equal,  of  course, 
to  the  direction  cosines  of  the  *  axis  in  the  system  of  coordinates 
5ns . 


Expanding  the  mentioned  determinant,  we  arrive  at  the  following 
equality : 


l^mm  nb —  me,  h%  •»  te — no,  kfmma  —  lb. 


(10) 


Moduli  of  both  vectorial  products  (7)  and  (9)  are  identical 
and  equal  to  the  sine  of  angle*a  between  axes  C  and  d  or,  which  Is 
the  same,  between  axes  *  and  d.  Thus,  taking  into  account  still 
another  equation  (1),  we  have 

M  m,  h  mm  »!•  <t  —  1  —  &  (11) 


By  means  of  equalities  (8),  (10)  and  (11)  relation  <  *0  can  be 
presented  In  the  form 

(1—  /*)  cos  f  ™  —  t(mb  +  nc)  -t  (**-fm*)  a.  (jj) 


It  is  the  deficient  third  equation,  besides  equations  (2)  and  (3), 
for  the  search  of  the  three  unknown  cosines  a,  b  and  a,  which 
determine  the  direction  of  the  x  axii'  with  respect  to  axes  of  the 
system  of  coordinates  SnC  . 


t 
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According  to  equation  (3),  we  have 


«*  +  iw-;<i -«).  0.3) 

Furthermore , 

!*•  (1H) 

Taking  account  the  last  two  equalities  in  equation  (12),  we  obtain 
after  reducing  similar  terms  to  the  equation  for  one  of  the  desired 
quantities 

*  »  eaay  +  ?  (1  —  mi  f).  (15) 

To  determine  quantity  b ,  we  exclude  from  relation  (2)  cosines 
a  and  e  by  means  of  equations  (12)  and  (15).  As  a  result  we  arrive 
at  the  quadratic  equation 

—  l*)eo«  ♦]•  +  *» +  ^.  (/—[!•  + 

+  (1  —  fcw)*  — I,  (l6) 

which  after  simplifications  with  the  help  of  relation  (1*0  is 
reduced  to  the  form 

*■—24  lm  (t—  cos?)  (I— cos  f)*— a*  «in*  f>™0.  (17) 

Of  the  two  roots  of  this  equation 


»,  ■ta(I  —  coaf)  +• »  ataf. 

*•  —  lm  (1  —  COB  —  n  cia  f 


e  should  discuss  the  first.  Actually,  assuming  specifically. 


}  k  a  o  0,  *  “■  J.  T”y  i 


(19) 


we  out-aln  accordingly 


-  1. 


(  20) 


However,  taking  into  account  the  accepted  designations  (see  Fig.  Ik), 
we  have  in  this  case 

b  -  eon  *n  “  -f  1  (21) 

and,  consequently,  the  second  root  of  equations  (18)  should  be 
dropped.  Thus,  in  general 

b  -  6»(i  —  e«9)  +»»«lnf.  (22) 

Substituting  now  expressions  (15)  and  (22)  for  a  and  b  into  equation 
(3),  we  will  arrive  at  the  following  equation  for  the  determination 
of  quantity  q,  namely: 


e  =  In  (1  —  (Mf)  —  m  sfnf.  (23) 

Similarly  cosines  of  angles  between  the  y  axis  and  axes 
r. ,  c  and  further  between  the  a  axis  and  the  same  axes  £,  n,  £  are 
determined.  Finally,  the  table  of  cosines  of  the  angles  between 
axes  of  the  system  of  coordinates  xyz  and  can  be  represented 

iri  the  following  form: 


I 

•  (I  —  eosf)/*-f 

+  «»? 


V  (1 —  eon  ?}/*»— 

*  (1  —  eos  f)  In  + 

+  M«ia? 


*1 

(1  —  cos  9)  m/q 
+  Jiaiaf 
(1  —  cosf)m*  + 
+  COST 

(1  —  cos  %>)  mn  — 
— -/slnf 


c 

(1  —  cosf)  ni  — 

—  m  sin  f 

(i  —  cos  V)  *»*  + 
+  /  sin  * 

(F  —  eos<p)n*  -f- 
+  cos<p 


2k) 


”*•  is  table  was  used  in  the  second  chapter  of  this  i  :'ok . 
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TtMfSTer  — — — 

This  monograph  discusses  the  mathematical  principles  of  certaii! 
allowed  variants  of  inertial  guidance  of  the  flight  of  ballis+i 
missl  les.i.e. ,  control  without  the  use  of  any  extern*-.!  in  format 
(radio  signals,  radiation  of  stars  and  others).  It  is  assumed 
that  the  providing  of  the  assigned  flight  range  of  the  missile 
is  produced  as  a  result  of  the  well-timed  switching-cj.  f  of  itr. 
engine  by  a  signal  entering  from  the  computer.  Fed  into  the  in- 
nut  of  the  latter  are  readings  of  sensitive  elements  of  the 
system  of  inertial  guidance,  and  these  measure  the  apparent 
acceleration  of  the  missile  or  integrals  from  the  apparent  tame 
acceleration.  Control  should  be  such  that  the  deviation  of  the 
actual  motion  of  the  missile  from  the  rated  does  not  r"iva  an 
effect  on  the  range  of  its  flight.  To  accomplish  thi  .  •  .oasurlr . 
instruments  are  placed  on  the  missile  in  a  quite  definite  manner 
and,  in  particular,  the  direction  of  their  axes  of  st^oltivity 
in  a  number  of  cases  is  stabilized  by  means  of  gyrosc  ;p:i- • .  Thi^ 
monograph  is  intended  for  specialists  in  uhe  field  of  ^hc  theory 
of  control  processes.  It  can  be  useful  in  the  investigation  of 
r:  ’w  problems  of  this  discipline  and  also  as  a  training  manual. 
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